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ABSTRACT 


The  influence  of  polarization  filtering  on  maximum 
detection,  recognition,  and  identification  ranges  of  a 
generic  second  generation  FLIR  sensor  is  examined  with  a 
computational  model.  The  scenario  studied  represents  a 
second  generation  FLIR  sensor  mounted  on  an  aircraft  in 
level  flight  at  3  00m  approaching  a  ship  target.  The  target 
ship  radiant  signature  is  modeled  with  an  advanced  infrared 
signature  prediction  program,  MuSES  (Multi- Service  Electro- 
Optic  Signature) .  A  weather  file  representative  of 
Midlatitude  Summer  at  sea  conditions  was  utilized. 
Polarized  sea  background  and  path  radiance  calculations  are 
performed  with  a  polarized  version  of  the  SEARAD  Radiance 
and  Propagation  Code.  Results  showed  that  there  is  an 
improvement  in  maximum  range  of  the  sensor  for  detection, 
recognition,  and  identification  tasks  when  a  horizontal 
filter  is  included,  provided  that  the  target  does  not  have 
a  negative  degree  of  polarization.  For  detection  task  the 
improvements  were  found  to  be  33.48%,  35.65%,  and  39.78-6 
when  the  target  has  0%,  +2%,  and  +8%  degree  of  polarization 
respectively .  A  better  modeling  of  Apparent  Temperature 
Difference  (ATD)  calculation  is  also  developed.  To  improve 
the  model  use  of  polarized  target  model  is  recommended. 
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X.  INTRODUCTION 


Infrared  systems  have  been  a  powerful  technique  to 
detect  objects  in  a  remote  and  passive  way  for  some 
decades .  The  absence  of  an  active  source  to  perform 
detection  makes  this  mechanism  widely  suited  for  military 
application.  Due  to  the  large  number  of  potential 
applications,  intensive  scientific  and  technological 
efforts  have  been  made  to  improve  the  performance  of  these 
systems . 

Thermal  imaging  systems  respond  to  a  difference  in 
radiant  exitance  between  a  target  and  its  background.  They 
are  mostly  sensitive  enough  to  detect  weak  targets,  but 
background  radiation  hampers  the  process,  making  detection 
difficult.  Proper  modeling  of  background  and  target 
characteristics  is  considered  as  a  key  factor  to  develop  a 
successful  target  detection  technique.  One  such  promising 
model  is  use  of  properly  polarized  filters  in  sensors  for 
background  radiance  suppression. 

The  phenomenon  of  polarization  in  the  infrared  has 
important  potential  applications  for  cueing  and  detecting 
targets  in  forward  looking  imagery  by  improving  the  target- 
to-background  contrast.  Previous  experiments  and  analysis 
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done  at  the  Naval  Postgraduate  School  have  shown  that  sea 
surface  radiance  has  a  vertically  polarized  characteristic, 
particularly  at  viewing  angles  within  a  few  degrees  of  the 
horizon.  [Ref.  1,  16].  As  a  consequence  of  these  studies 
the  concept  of  improvement  in  target-to-clutter  background 
contrast  ratio  by  implementing  polarization  filters  gained 
interest.  However  having  a  polarized  filter  in  the  sensor 
also  effects  the  incoming  radiation  from  a  target,  which 
may  have  apparent  polarized  features. 

The  main  objectives  of  this  project  are  to  determine 
the  effect  of  degree  of  polarization  on  maximum  range 
performance  analysis  using  an  advanced  signature  prediction 
model,  and  also  to  verify  previous  project  results  and 
models  which  have  studied  target-background  contrast 
improvement  using  horizontal  polarizers  on  sensors.  In 
Chapter  II  some  fundamentals  of  infrared  theory  are 
presented.  Chapter  III  is  dedicated  to  background- target 
radiation  theory  and  polarization  theory.  In  this  Chapter 
beginning  from  the  basic  radiation  theories,  all  the  known 
effects,  such  as  environmental  radiation  components,  are 
studied.  The  polarization  effect  on  background  and  target 
signatures  is  also  studied  in  detail. 


In  Chapter  IV  descriptions  of  the  scenario  and  the 
model  are  given  and  using  the  MuSES  (Multi-Service  Electro- 
Optic  Signature) [Ref.  25]  infrared  signature  prediction 
program,  target  signatures  for  the  appropriate  conditions 
are  computed. 

In  Chapter  V,  computations  of  Minimum  Resolvable 
Temperature  Difference  (MRTD)  of  a  second-generation 
sensor,  Apparent  Temperature  Difference  (ATD)  between 
target  and  background,  and  sea  background  radiance  are 
made.  Following  definition  of  the  degree  of  polarization, 
the  effect  of  polarization  is  included  into  the  overall 
computations.  Having  these  results,  maximum  range  analysis 
for  a  second-generation  FLIR  (Forward  Looking  Infrared 
System)  sensor  with  a  horizontally  polarized  filter  is 
performed  by  plotting  the  MRTD  and  ATD  values  on  the  same 
graph.  Finally  conclusions  and  recommendations  are  provided 
in  Chapter  VI . 
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II.  INFRARED  TECHNOLOGY  FUNDAMENTALS 


This  Chapter  summarizes  the  basic  infrared  theory  that 
includes  some  important  points  necessary  for  full 
understanding  of  this  work. 

A.  ELECTROMAGNETIC  SPECTRUM 


The  propagation  of  electromagnetic  energy  can  be 
described  in  terms  of  spatial  and  temporal  variations  in 
electric  and  magnetic  fields.  Typical  examples  of 
electromagnetic  waves  are  light  rays,  radio  waves,  radar 
beams  and  infrared  radiation.  The  optical  radiation  covers 
the  electromagnetic  spectrum  from  the  Ultraviolet  (UV) 
range  throughout  the  visible  and  up  to  the  Infrared  (IR) 
portion  of  the  spectrum.  Figure  2.1  shows  the  regions  of 
the  spectrum. 
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Figure  2.1  The  Electromagnetic  Spectrum.  From  Ref . [3] 
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The  infrared  portion  of  the  spectrum  lies  between  the 
visible  (shorter  wavelength)  and  the  radio  frequency 
(longer  wavelength)  regions.  This  region  corresponds  to  the 
range  of  wavelengths  from  approximately  0.7pm  to  1000pm.  The 
near  infrared  (NIR)  is  characterized  by  wavelengths  from 
0.7pm  to  0.9pm,  and  the  short-wave  infrared  (SWIR)  is 
characterized  by  wavelengths  from  0.9pm  to  3.0pm.  Since  the 
short-wave  band  includes  mostly  reflected  radiation  it  can 
be  imaged  with  night  vision  devices  that  collect  and 
amplify  the  reflected  light  waves.  The  mid-wave  (MWIR)  is 
characterized  by  wavelengths  from  3.0pm  to  6.0pm  and  the 
long-wave  infrared  (LWIR)  characterized  by  wavelengths  from 
6.0pm  to  15.0pm.  The  long-wave  band  is  called  the  "thermal 
radiation  region"  because  thermal  emission  is  greater  than 
the  reflection  in  these  wavelengths .  Due  to  this 
characteristic  of  the  mid-wave  and  long-wave  bands,  thermal 
imaging  systems  such  as  FLXR  (Forward  Looking  Infrared 
System),  IRLS  (Infrared  Line  Scan),  and  IRST  (Infrared 
Search  and  Track)  work  in  these  bands . 

The  factors  such  as  the  form  of  the  radiation,  methods 
of  detection,  and  the  atmospheric  transmission  mostly 
determine  these  sub-bands .  The  region  from  5pm  to  8pm  is 
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also  a  part  of  the  infrared  spectrum,  but  it  is  not  used 

due  to  the  high  absorption  of  the  atmosphere  in  this  area. 

The  infrared  radiation  complies  with  the  laws  of 
electromagnetism,  as  stated  by  Maxwell's  Equations. 

B.  INFRARED  DEFINITIONS  AND  UNITS 

Since  infrared  science  makes  use  of  many  specific 

definitions  and  uses  varied  units,  it  is  necessary  to 
introduce  a  terminology  for  the  sake  of  better 

understanding.  Table  2.1  shows  some  fundamental  radiometric 
quantities  and  units  associated  with  infrared  radiation.  In 
this  table: 

•  Radiant  exitance  (emittance)  or  flux  density,  M,  is 
the  radiant  flux  density  leaving  the  target  surface 
per  unit  area. 

•  Irradiance  or  radiant  flux  surface  density,  E,  is 
the  radiant  power  incident  on  a  unit  area  of 
surface. 

•  Radiant  intensity,  I,  is  the  radiant  power  exiting 
a  point  source  along  a  given  direction  within  a 
unit  solid  angle. 

•  Radiance,  or  radiant  intensity  surface  density  in 
a  given  direction,  L,  is  the  radiant  power  per  unit 
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solid  angle  per  unit  area  leaving  from  or  incident 


on  an  area  projected  perpendicular  to  the  direction 
of  radiant  energy  flow. 

•  In  table  U  represents  radiant  energy,  S  represents 
surface  area,  and  £1  (sr)  represents  solid  angle. 


Symbol 

Quantity 

Units (metric) 

Description 

U 

Radiant 

Energy 

Joules 

P 

Radiant 

Power 

Watts 

9u/9t 

M 

Radiant 

Exitance 

Watt/cm^ 

9p/9s 

E 

Radiant 

Incidance 

Watt /cm2 

9p/9s 

I 

Radiant 

Intensity 

Watt/sr 

9p/9£1 

L 

Radiance 

Watt/cm2  sr 

92p/  (9s9Q) 

Table  2.1  Radiometric  Quantities. 


C.  LAWS  OF  THERMAL  RADIATION 
1 .  Blackbody  Radiation 

When  we  consider  the  total  optical  power  that  is 
incident  on  a  object,  we  can  distinguish  three  distinct 
outcomes .  Some  fraction  of  this  total  radiant  energy  is 
absorbed,  some  is  reflected  and  some  is  transmitted  through 
the  object.  Therefore  the  ratios  of  each  of  these  to  the 
incident  power  must  add  up  to  unity.  [Ref.  4] 


a+p+t= 1 


(2.5) 


where : 

a  (absorptivity)  =  Pabsorbed/Pincident 
P  (ref  lectivity)  =  Preflected/Pincident 
X  ( transmissivity )  =  Ptransmitted/P  incident 

An  ideal  blackbody  can  be  defined  as  an  ideal  emitter 
and  an  ideal  absorber  i.e.  absorptivity (a) =emissivity (£)  =1 . 
Thus  the  blackbody  has  a  reflectivity  (p)  and  transmissivity 
(x)  both  of  zero.  The  blackbody  radiator  is  an  idealized 
source  of  radiant  energy  having  defined  properties.  It  is 
perfectly  diffuse,  radiates  at  all  wavelengths,  and  at  all 
temperatures  its  spectral  radiant  exitance  is  the  maximum 
possible  for  any  actual  thermal  source  at  the  same 
temperature.  The  characteristics  of  many  real  sources  of 
radiation  approach  the  ideal  quite  closely;  many  others  are 
conveniently  described  by  corrections  to  the  blackbody 
curve. [Ref.  2]  The  spectral  radiant  exitance  of  a  blackbody 
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is  presented  in  Figure  2.2,  as  a.  function  of  wavelength  and 


source  temperature. 


Figure  2.2  Spectral  Radiant  Exitance  of  Blackbodies  at  Various 
Temperatures.  Dotted  Curve  Represents  the  Displacement  Law. 
From  Ref. [23] . 


2.  Planck 's  Law 

The  equation  that  represents  a  blackbody's  spectral 
exitance  as  a  function  of  its  temperature  is  known  as 
Planck's  Radiation  Law  and  has  the  following  form: 


Irihc 


X5{ehc'm  -l) 


xlO 


-6 


(2.6) 
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where : 


M\,b  =  the  blackbody  spectral  radiant  exitance  at 
wavelength  X  (Watt/m2  |im) 

c  =  The  vacuum  velocity  of  light  (2.9979250  x  108  m/s) 
h  =  Planck's  constant  (6.626196  x  10  34  Joule/s) 
k  =  Boltzmann's  constant  (1.380622  x  10~23  Joule/K) 

T  =  Absolute  temperature  of  the  blackbody  (K) 

X  =  Wavelength  (m) 

In  infrared  theory  the  radiant  exitance  is  directly 
proportional  to  the  emissivity.  The  " emissivity" ,  £,  of  a 

surface  is  the  ratio  of  radiant  exitance  emitted  from  a 
surface  to  that  emitted  from  a  blackbody  (perfect  emitter) 
at  the  same  temperature.  The  radiant  exitance  of  a  non 
blackbody  then  can  be  defined  as 

Mx  (2.7) 

where  £*.  is  the  spectral  emissivity  which  is  always  less 
than  1  for  non-black  bodies.  [Ref.  2]  A  material  which  has 
an  emissivity  that  is  independent  of  wavelength  is  called  a 
gray  body.  Table  2.2  shows  emissivities  of  some  common 
materials . 
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gray  body.  Table  2.2  shows  emissivities  of  some  common 
materials . 


Material 

£ 

Aluminum  Foil 

0.09 

Copper , polished 

0.05 

Copper , oxidized 

0.78 

Carbon 

0.95 

Paint, oil 

0.94 

Concrete 

0.92 

Sand 

0.90 

Table  2.2  Emissivities  (8)  of  Common  Materials.  From 
Ref. [4]  . 

3 .  Wien  Displacement  and  Stef  an-Boltzmann  Laws 

Figure  2 . 2  shows  blackbody  curves  in  terms  of  the 
radiant  exitance  as  a  function  of  wavelength  for  a  number 
of  different  temperatures,  each  one  of  higher  temperature 
completely  above  the  others.  The  wavelength  of  the  peaks, 
i.e.  the  wavelengths  where  the  maximum  radiation  occurs  for 
a  given  temperature,  can  be  found  from  the  eguation  known 
as  the  Wien  Displacement  Law, 

_  2898  (2.8) 

'Knax  rp 

where : 

A,„ax  =  wavelength  where  the  peak  of  radiation  occurs  (Jim)  . 

The  peak  wavelength  in  micrometers  is  approximately 
3000  divided  by  the  temperature  in  Kelvin. 
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The  Planck  curve  can  be  integrated  to  obtain  the 


expression  for  total  radiant  exitance  at  all  wavelengths . 


Wt=<fT‘ 


(2.9) 


where : 

Wb  =  Total  radiant  exitance  of  a  blackbody 

—  8 

O  =  The  Stefan-Boltzmann  constant  (5.6691  X  10 
Watts/m2K4  ) 

T  =  The  temperature  in  K 

For  use  with  non-black  body  sources  this  law  is  modified  by 
the  inclusion  of  an  "effective  emissivity" ,  £,  giving  the 

form 


Wh  =  eoT4 

This  equation  is  known  as 


(2.10) 

the  Stefan-Boltzmann  Law. 


D.  ATMOSPHERIC  PROPAGATION  OF  INFRARED  RADIATION 

The  earth's  atmosphere  consists  of  different  molecules 
and  suspended  particles.  When  infrared  radiation  passes 
through  the  atmosphere  these  molecules  and  suspended 
particles  both  absorb  and  scatter  radiation.  The  degree  and 
type  of  the  attenuation  depend  on  the  properties  of  the 
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atmosphere  and  the  distance  between  the  source  of  the 
radiation  and  the  receiver. 

The  gases  and  particles  suspended  in  the  path  may 
themselves  radiate.  This  path  radiance  and  its  fluctuations 
may  "blur"  the  scene  and  reduce  the  contrast.  Scattering 
agents  in  the  path  may  also  deflect  or  scatter  radiation 
from  other  sources  in  such  a  way  that  the  source  radiance 

appears  to  be  increased.  [Ref.  7] 

The  two  main  causes  of  atmospheric  attenuation  of 
infrared  radiation,  however,  are  absorption  by  molecules  in 
the  atmosphere  and  scattering  due  to  particles  suspended  in 
the  path.  Beer's  law  quantifies  the  total  attenuation  at 
any  wavelength  by  the  atmosphere  as: 


where 

•  I0  is  the  initial  radiant  intensity 

•  lz  is  the  radiant  intensity  of  radiation  after 
propagation  distance  z  in  the  atmosphere 

•  (I  is  the  extinction  coefficient  of  atmosphere,  which 
is  defined  as  the  sum  of  absorption  and  scattering 
coefficients 

•  z  is  propagation  distance 
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1.  Absorption 


Absorption  of  radiation  within  the  atmosphere  can 
occur  in  electronic,  vibrational  or  rotational  transitions. 
Transitions  between  electronic  states  in  an  atom  generally 
take  place  in  the  visible  or  ultraviolet  region  of  the 
spectrum.  Infrared  molecular  absorption  is  due  to 
transitions  between  vibrational  and/or  rotational  energy 
states  of  a  gas  molecule.  [Ref.  5] 

Since  the  atmosphere  contains  many  different  particles 
and  behaves  in  a  dynamic  way  it  is  not  very  easy  to  predict 
the  extinction  coefficient  and  transmittance.  An  accurate 
atmospheric  transmission  calculation  requires  that  all 
molecular,  aerosol,  and  precipitation  effects  be  taken  into 
account.  When  exact  results  are  required  a  detailed  model 
must  be  used.  In  infrared  technology  computer  models  such 
as  the  LOWTRAN,  MODTRAN,  FASCODE  and  the  SEARAD  Radiance 
Models  are  used  to  compute  the  extinction  coefficient  and 
transmittance . 

In  this  project  the  latest  version  of  M0DTRAN4  and  the 
SEARAD  Radiance  Model  are  used  for  transmittance,  path 
radiance,  and  background  radiance  calculations.  An  example 
of  transmittance  values  obtained  for  a  1-km  path  length  is 
shown  in  Figure  2.3. 
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wavelength  -  micrometer 


Figure  2.3  Typical  Atmospheric  Transmission  for  a  1-km 
Length  and  Midlatitude  Stammer  Conditions.  After 
Ref. [24] . 

2 .  Scattering 

Scattering  is  a  process  that  changes  the  direction  of 
motion  of  individual  photons,  causing  divergence  of  the 
energy  and  decrease  in  the  forward  radiance. 

a.  Rayleigh  Scattering 

The  decrease  in  forward  spectral  radiance,  due 
to  scattering,  depends  on  the  wavelength  of  the  incident 
radiation  and  the  diameter  of  the  particles.  When  the 
particles  are  small  compared  with  the  wavelength  of  the 
radiation,  the  process  is  known  as  Rayleigh  scattering.  In 
this  case  the  radiance  loss  due  to  scattering  decreases 
rapidly  with  increase  in  wavelength. 
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In  the  infrared  spectrum,  scattering  is  greatest  in 
the  near  infrared  region  adjoining  the  visible  spectrum, 
and  decreases  rapidly  for  longer  wavelengths.  [Ref.  6] 

b.  Mie  Scattering 

For  larger  particle  sizes  greater  than  about  one 
tenth  of  the  wavelength,  the  Rayleigh  approximation  is  no 
longer  valid,  and  another  scattering  theory,  Mie  theory, 
applies.  [Ref.  2]  These  classes  of  particles  generally 
cause  haze,  rain,  and  fog  conditions.  Since  smoke  and  light 
mist  particles  are  small  compared  to  the  infrared 
wavelengths,  infrared  radiation  can  propagate  further 
through  smoke  and  mist  than  visible  radiation.  However, 
fog,  and  rain  particles  are  larger  and  as  a  result  they 
scatter  infrared  and  visible  radiation  to  a  similar  degree. 

In  the  atmosphere,  transmission  of  infrared  radiation 
is  attenuated  by  both  absorption  and  scattering,  which 
occur  together.  Absorption  may  be  much  greater  than 
scattering,  or  vice  versa,  depending  upon  the  nature  of  the 
atmosphere,  the  particle  size,  and  the  wavelength. 
Furthermore,  the  combined  sum  of  the  percentages  of 
absorption,  scattering,  and  transmission  always  totals  100. 
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E.  INFRARED  POLARIZATION  PHYSICS 

Polarization  has  recently  received  increased  interest 
in  discrimination  of  targets  from  background,  since  it 
provides  more  information  than  pure  intensity  imaging.  When 
the  intensity  contrast  between  target  and  background  is  not 
distinguishable,  the  polarization  contrast  may  offer  new 
dimensions  in  discrimination.  [Ref.  17]  In  this  section  the 
polarization  physics  is  summarized  while  the  effects  of 
polarization  phenomena  in  thermal  imaging  sensors' 
performance  and  target  and  background  thermal  signatures 
are  studied  in  detail  in  Chapter  III. 

For  a  plane  wave  propagating  in  the  z-direction,  the 
electric  field  components  in  the  plane  perpendicular  to  the 
direction  of  propagation  may  be  represented  by. 

Ex(z,t)  =  Eoxcos(wt  +  Sx)  (2.12) 

Ey  ( z,t)  =  Eoy  cosiwt  +  Sy)  ( "  ‘  X  ) 

Where  Ex  =  x-component  amplitude 
Ey  =  y- component  amplitude 
8X  =  phase  constant  of  x  component 
8y  =  phase  constant  of  y-component 

The  polarization  of  electromagnetic  radiation  is 

defined  by  the  behavior  of  the  electric  vector  over  time 
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intervals  that  are  long  compared  to  the  response  time  of 
the  detector  to  be  used. 

When  the  end  point  of  the  electric  vector  at  each 
point  in  the  plane  of  incidence  moves  periodically  around 
an  ellipse  with  increasing  time,  the  radiation  is  called 
"elliptically  polarized".  The  ellipse  reduces  in  special 
cases  to  a  circle  (circularly  polarized)  or  a  straight  line 

(linearly  polarized).  [Ref.  12] 

It  is  known  that  the  polarization  behavior  could  be 
completely  represented  in  terms  of  four  measurable 
quantities  known  as  the  Stokes  polarization  parameters.  The 
Stokes  polarization  parameters  for  a  plane  wave  are: 

[Ref.  2] 


I  =  E2  ox  +  E2  oy 

(2.14) 

■*i  __o 

(2.15) 

M  =  E2ox  -  E  oy 

C  =  2 EmE  oy  cos  8 

(2.16) 

S  =  2 EoxE  oy  sin  6 

(2.17) 

where 

1=  the  total  intensity  of  the  light 
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M=  the  degree  of  polarization  in  the  horizontal  or  vertical 
direction 

Eox=  the  component  of  the  wave  electric  field  normal  to  the 
plane  of  propagation 

Eoy=  the  component  of  the  wave  electric  field  parallel  to 
the  plane  of  propagation 

C=  the  degree  of  polarization  in  the  ±45°  direction 

S=  the  degree  of  polarization  of  right  or  left  circular 

polarization  contained  within  the  beam 

5=  the  phase  difference  between  the  s  and  p  components  of 
the  electric  field 

Unpolarized  light  has  M=C=S=0 .  To  describe  completely 
the  polarization  of  light,  all  four  of  Stokes  parameters 
must  be  used. 

An  important  quantity  that  collectively  describes  the 
parameters  is  the  degree  of  polarization  POL.  It  can  be 
calculated  by: 


POL  = 


>+("»> 


(2.18) 


where : 

<NV>  =  mean  apparent  vertically  polarized  radiance. 
<Nh>  =  mean  apparent  horizontally  polarized  radiance. 
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The  observed  degree  of  polarization  depends  on  several 
parameters  such  as  the  refractive  index  of  the  surface,  the 
viewing  angle,  the  surface  temperature,  the  radiating 
environment,  surface  roughness  and  so  on. 

Polarization  of  light  can  also  be  generated  by 
reflection  and  refraction.  The  polarization  of  the  wave  can 
change  while  passing  from  one  medium  to  another.  This 
phenomenon  including  the  reflection  and  refraction 
characteristics  of  light  is  governed  by  the  Fresnel 
Equations  and  by  Snell's  Law.  For  nonmagnetic  media  the 
Fresnel  equations  can  be  written 


sin(6>  -  9t) 

(2.19) 

Wndichr  -  sin(6,  +0J 

tan(0;  -6t) 

W‘el  tan(6>W) 

(2.20) 

2  •  sin(^)  •  cos($) 

Wndic.hr-  sin(0.+0) 

(2.21) 

2  •  sin(#r)  •  cos($) 

(2.22) 

sii&'+ej-cos&'-e,) 

where : 

0i  =  incident  angle. 

0t  =  angle  of  transmission, 
r  =  reflection  coefficient. 
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t  =  transmission  coefficient. 

"Perpendicular"  and  "parallel"  refer  to  the  component 
fields  perpendicular  and  parallel  to  the  plane  of 
incidence . 

For  plane  waves  incident  at  a  boundary  between  two 
media  of  different  complex  refractive  indices  nx  and  n2,  the 
angle  of  reflection  equals  the  angle  of  incidence.  The 
angle  of  transmission  (refraction)  0t  as  a  function  of  the 
incidence  angle  0i  is  given  by  Snell's  law: 

ttj  sin(#,)  =  n2  sin(0, )  (2.23) 
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III.  BACKGROUND  -  TARGET  RADIANCE  AND  POLARIZATION  EFFECTS 


A.  BACKGROUND  RADIANCE 

1 .  General 

In  an  environment  where  there  are  some  optical  targets 
present,  the  surrounding  elements  that  emit,  reflect  or 
scatter  random  radiation  in  the  same  spectral  region  as  the 
target  can  be  defined  as  "background" .  There  are  numerous 
background  infrared  radiation  sources  occurring  in  nature 
that  can  be  divided  into  three  groups;  terrestrial, 
celestial  and  atmospheric  sources.  Terrestrial  background 
refers  to  radiance  seen  by  infrared  systems  working  within 


the  atmosphere 

of  the  earth. 

This  type  of  background 

includes 

the 

radiation 

that 

several 

natural  sources 

occurring 

on  the  surface 

of 

the  earth 

are  emitting  or 

reflecting  (oceans,  clouds,  vegetation,  etc.).  It  also 
includes  sunlight  scattered  by  particles  in  the  atmosphere 
and  atmospheric  molecular  emission  (sky  radiance) . 
Celestial  refers  to  that  background  seen  by  a  system 
working  in  space  and  consists  of  radiation  from  the  stars 
and  the  sun,  reflectance  from  planets,  and  space  debris. 
Dust  particles,  water  droplets,  and  molecules  of  the 
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various  gases  occurring  in  the  earth's  atmosphere  behave 
like  background  sources  of  radiation  with  respect  to  an 
optical  target.  They  scatter,  reflect,  absorb  and  re-emit 
primary  radiation  illuminating  them  from  both  celestial  and 
terrestrial  sources. 

Background  radiation  is  strongly  affected  by 
environmental  conditions  and  the  nature  of  the  materials . 
Changes  in  emissivity  and  reflectance  of  the  various 
background  surfaces  under  varied  conditions  are  considered 
primary  parameters  that  effect  the  background  radiation  in 
terms  of  magnitude  and  polarization  characteristics. 

Since  the  target  is  always  surrounded  by  some  sort  of 
background  it  is  obvious  that  radiation  reaching  the 
detector  from  this  background  is  confused  with  the 
radiation  emitted  or  reflected  from  the  target.  [Ref.  5] 

Background  radiance  measurements  are  generally  made 
with  thermal  imagers  and  radiometers  under  various 
meteorological  conditions.  Technical  properties  of  imagers 
and  radiometers  sometimes  make  the  comparison  of  these 
devices  very  difficult  in  terms  of  thermal  resolution.  The 
geographical  locations  of  measurements  also  may  change . 
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All  of  the  above  make  the  background  phenomena 
complicated  and  the  modeling  of  a  universal  background 
radiance  function  very  difficult. 

Of  particular  interest  in  this  project  are  the 
background  radiances  of  the  sea  surface  and  the  sky.  In  the 
rest  of  the  Chapter,  the  main  characteristics  and  spectral 
radiance  properties  of  these  backgrounds,  especially  for 
the  spectral  region  8pm  to  14pm  in  which  we  are  most 
interested,  are  reviewed. 

2 .  Sea  Background 

The  sea  surface  infrared  radiance  Lsea  is  the  sum  of  the 
radiance  caused  by  its  inherent  blackbody  thermal  emission 
Lsst  from  below  the  surface  and  the  radiance  due  to 
reflected  incident  sky  and  solar  radiation  from  the 
surface.  This  can  be  written  in  the  following  form. 

4ea  =  d'sst  +  PLsky  +  P\m  (3.1) 

where 

£  =  the  sea  surface  emissivity 

p  =  the  sea  surface  reflectivity 

p'  =  reflectivity  integrated  over  the  angle  subtended  by 
the  sun . 

All  of  these  are  averaged  over  the  wave  slope  distribution. 
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Equation  3 . 1  must  then  be  corrected  for  atmospheric 
transmittance  and  path  radiance  over  the  range  "R"  so  that 
the  apparent  radiance  received  Lapp  can  be  written: 


Lapp  ( R ,  Z,  d,  <j),  £,  p,  t)  =  Lsea ( T ,  Z,  e,  </>,  £,  p)  *  t(R,  Z,  e,  (t>,  t )  +  Latm ( R ,  Z,0,<p,t )  (3.2) 


where 

R  =  range  to  the  sensor 
0  =  azimuthal  angle 
<])  =  elevation  angle 
X  =  wavelength 
t  =  time 

Lsea  =  total  surface  or  zero  range  background  radiance 
Latm  =  radiance  emitted  along  the  path  to  the  sensor 
x  =  transmittance  through  the  atmospheric  path  to  the 
sensor . 

[Ref.  1] 

Major  factors  that  determine  the  sea  background 
character  can  be  summarized  as;  Infrared  optical  properties 
of  sea  water,  geometry  of  the  sea  surface  and  wave  slope 
distribution,  and  finally  the  sea  water  temperature 
distribution  and  properties  of  the  bottom  material. 


26 


Each  radiance  component  except  path  radiance  is  a 
function  of  the  spectral  emissivity  or  reflectivity  (£x,px) 
of  sea  water  and  the  orientation  of  the  sea  wave  facet 
relative  to  the  reflected  sources  and  the  observer.  Since 
£x(6)  +px(9)  =1/  both  reflectivity  and  emissivity  will  depend 
on  the  polarization  of  the  incident  or  emitted  radiation. 

The  sea  radiance  may  be  modeled  as  a  blackbody  at  the 
sea  surface  temperature.  The  reflected  radiance  component 
is  integrated  over  the  sky  zenith  and  azimuth  angles. 
[Ref. 16]  Direct  reflected  solar  radiation  is  significant 
only  in  very  restricted  orientations  close  to  the  solar 
direction  and  dominates  at  short  wavelengths. 

To  compute  the  radiance  from  an  ocean  surface,  a 
FORTRAN  computer  code,  the  SEARAD  Radiance  Model  developed 
at  NCCOSC  can  be  used.  The  SEARAD  Radiance  Code  [Ref.  10] 
is  a  modification  of  MODTRAN  valid  for  a  spectral  range 
extending  from  the  visible  to  the  far  infrared.  The  latest 
SEARAD  code  takes  into  account  the  polarization 
characteristics  of  sea  surface  emission  and  reflected  sky 
radiance  components  in  calculations .  Several  variously 
modified  versions  exist  such  as  a  user-friendly,  multi¬ 
range  version  that  runs  in  MATLAB  and  a  self-contained  DOS- 
compatible  version  that  runs  in  the  DOS  environment.  The 
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SEARAD  Radiance  model  computes  four  contributions  (thermal 
emission,  path  radiance,  reflected  sky  radiance,  and 
reflected  solar  radiance)  to  sea  radiance.  It  is  based  on 
the  Cox-Munk  (Cox  and  Munk,  1954,1956)  [Ref.  26] 
statistical  model  for  wind-driven  capillary  wave  facets. 
The  assumption  for  the  model  is  that  the  strength  of 
interaction  between  an  optical  ray  and  a  capillary  wave 
facet  is  proportional  to  the  facet  area  projected  normal  to 
the  ray . 

SEARAD  also  allows  computation  of  the  atmospheric 
transmittance  as  in  normal  MODTRAN.  It  is  used  in  this 
project  many  times  for  sea  background  radiance  calculations 
and  for  part  of  the  Apparent  Temperature  Difference 
calculations.  Sample  SEARAD  input  and  output  data  files  are 
presented  in  Appendix  A. 

3 .  Sky  Background 
a .  General 

The  sky  background  radiance  in  the  infrared  is 
caused  by  radiation  emitted  from  the  molecules  in  the 
atmosphere  and  by  sunlight  scattered  from  the  molecules, 
particles,  and  water  droplets  in  the  atmosphere.  The 
scattering  factor  is  especially  influenced  by  the  cloud 
distribution.  It  is  important  in  the  short  wavelength 


range,  close  to  the  visible,  and  is  present  only  in 
day time .  The  thermal  emission  happens  at  wavelengths  longer 


than  4  ^m  and  is  present  day  and  night. 

An  idealized  representation  of  the  sky  spectral 
radiation  is  shown  in  Figure  3.1  [Ref.  4] 


Figure  3.1  Contribution  From  Scattering  and  Atmospheric 
Emission  To  Background  Radiation.  From  Ref. [4]. 
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b.  Spectral  Radiance  of  Clear  Sky 

Spectral  radiant  intensity  and  spectral 
distribution  of  the  clear  sky  in  the  infrared  region  depend 
on  angle  of  elevation,  water  vapor  and  carbon  dioxide 
contents  of  the  atmosphere,  the  temperature  of  the 
atmosphere,  and  the  altitude  above  the  earth's  surface. 
They  are  also  affected  by  clouds  and  haze. 

The  spectral  distribution  of  energy  radiated  from 
clear  day  and  night  skies  is  shown,  not  to  scale,  in  Figure 
3.2  which  illustrate  the  primary  differences  between  the 
day  and  night  sky  background  radiation  normal  to  the 
earth's  surface.  [Ref.  4] 


Figure  3.2  The  Spectral  Energy  Distribution  of 
Background  Radiation  from  Clear  Day  and  Night 
Skies.  From  Ref. [4] . 

The  elevation  angle  determines  the  length  of  the 
atmospheric  path  and  hence  the  emissivity,  while  the 
temperature  determines  the  blackbody  emission 
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characteristics.  The  elevation  angle  dependence  of  the 
spectral  radiance  is  shown  in  Figure  3.3. 

The  effect  of  the  ambient  temperature  is  very  strong 
for  this  kind  of  radiation,  as  shown  in  Figure  3.4.  This 
figure  shows  the  zenith  sky  spectral  radiance  for  two 
different  ambient  temperatures. 


Figure  3.3  The  Elevation  Angle  Dependence  of 
Sky  Spectral  Radiance.  From  Ref. [4] . 


Figure  3.4  Sky  Spectral  Radiance  Variation 
with  Ambient  Temperature.  From  Ref. [4]. 
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c.  Spectral  Radiance  of  an  Overcast  Sky  and 
Clouds 

The  infrared  radiation  from  clouds,  which  results 
as  clouds  emit  heat  and  scatter  sunlight,  is  an  important 
consideration  in  computations  and  produces  considerable 
variation  in  sky  background. 

Overcast  skies  have  been  found  to  be  very  good 

blackbody  radiators  (8=1),  having  a  temperature  at  low 
altitudes  within  one  or  two  °C  from  the  ambient  surface 
temperature.  Infrared  radiation  incident  on  a  cloud  water 
droplet  is  highly  absorbed  and  partially  reflected  and 
transmitted.  However  due  to  multiple  scattering  from  all 
neighboring  droplets  in  the  cloud  the  transmittance  of  the 
cloud  goes  essentially  to  zero,  and  thus  the  emissivity 
approaches  unity.  We  can  say  that  clouds  in  the  infrared 
region  are  good  absorbers  and  consequently  good  emitters, 
causing  overcast  skies  to  behave  as  good  blackbody 
radiators . 

4 .  Land  Background 

The  radiant  energy  emitted  and  reflected  by  the  ground 
is  determined  by  the  emissivity,  reflectivity,  and 
temperature  of  the  ground  and  the  modification  of  this 
radiation  by  the  intervening  atmosphere. 
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The  same  general  principles  that  govern  the  sky 
background  determine  the  amount  and  variation  of  land 
background  radiation.  The  radiation  from  the  terrain  is 
the  sum  of  the  reflected  sky  and  solar  radiation  plus  a 
percentage  of  the  radiation  emitted  by  a  blackbody  at  the 
ground  temperature.  When  the  observed  terrain  is 
sufficiently  distant,  the  radiation  is  modified  by  the 
atmosphere  at  the  wavelengths  where  it  absorbs . 

Below  4  | I'm  most  of  the  radiation  is  from  scattered  or 
diffusely  reflected  sunlight.  This  radiation  varies  greatly 
with  the  sun's  position,  the  amount  of  cloud  cover,  and  the 
nature  of  the  terrain. 

Above  4  Jim,  measurements  of  the  radiance  of  a  terrain 
show  that  it  closely  follows  the  appropriate  blackbody 
curve  for  the  temperature  of  the  ground.  A  few  minor 
deviations  from  the  blackbody  curve  may  occur  in  regions 
where  there  is  a  strong  atmospheric  absorption  band.  The 
diurnal  variation  of  the  radiance  from  various  types  of 
terrain  follows  the  temperature  changes  of  the  emitter. 
Differences  in  the  radiances  of  trees,  rocks,  grasses, 
cities,  and  other  types  of  terrain  are  easily  observed. 
These  result  from  differences  both  in  their  emissivities 
and  in  their  temperatures  under  natural  conditions.  The 
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radiance  of  the  terrain  is  higher  in  sunlit  areas  than  in 


areas  which  do  not  receive  sunlight  because  of  clouds . 
[Ref.  8] 

Figure  3 . 5  shows  the  surface  radiance  and  temperature 
of  a  land  background  at  two  different  times. 


Figure  3.5  The  Surface  Radiance  and  Temperature  of  a 
Land  Background  Compared  at  Morning  (0855)  versus  Pre-dawn 
(0135) .  From  Ref. [9] . 

Since  computation  of  radiation  from  backgrounds  is  one 
of  the  important  steps  in  overall  performance  analysis  of 
thermal  imaging  systems,  the  advanced  models  must  be  used 
to  make  these  computations.  In  the  sea  background  case  the 
SEARAD  Radiance  Model  is  considered  adequate  for  many 
applications . 

In  principle  there  are  three  different  methods 
available  to  model  land  surface  backgrounds.  These  are;  a 


pure  statistical  approach,  models  based  on  physical  laws, 
and  a  combination  of  both  methods .  Models  contain 
parameters  such  as  soil  type,  soil  moisture  content, 
vegetation  cover  and  type.  The  physical  properties  of  the 
medium  (soil  water  diffusivity,  hydraulic  conductivity, 
thermal  diffusivity  and  thermal  convectivity)  are  both 
temperature  and  moisture  dependent  and  should  be  considered 
in  the  models .  One  thing  common  to  all  the  different  models 
is  the  fact  that  they  only  give  mean  conditions .  No 
information  on  variances  and  the  bandwidth  of  values  found 
in  reality  can  be  obtained.  To  get  this  information  new 
methods  have  to  be  developed.  [Ref.  9] 

B .  TARGET  RADIANCE 

An  infrared  imaging  system  can  detect  a  target  when  the 
target  has  a  radiation  signature  that  can  be  distinguished 
from  all  the  background  radiation  within  the  field  of  view. 
If  the  target  is  an  intense  source  such  as  a  jet-engine  or 
rocket  exhaust,  this  can  be  achieved  more  easily.  On  the 
other  hand  if  the  target  emits  nearly  the  same  radiance  as 
the  background,  such  as  a  man  or  tank  on  the  ground,  than 
it  requires  more  work. 
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Infrared  signatures  of  objects  are  functions  of 
surface  temperature,  emissivity,  reflectance,  and  the 
surrounding  environment .  At  zero  range  the  radiance 
consists  of  an  emitted  component  and  several  reflected 
components . 

The  emitted  portion  of  the  signature  is  blackbody 
radiation,  reduced  by  the  directional  emissivity  of  the 
target  surface.  The  reflected  radiance  comes  from  other 
portions  of  the  target  and  the  surrounding  environment ( e . g . 
atmosphere,  sun,  moon,  and  other  sources) ,  all  acting 
through  the  reflectance  of  the  target  coating.  [Ref.  13]  . 
The  overall  contributions  to  the  signature  from  a  target 
are  shown  in  Figure  3.6. 


Figure  3 . 6  Infrared  Signature  Composition  of  a  Target . 

In  this  section,  equations  for  predicting  the  radiance 
at  a  target  observation  point,  using  the  fundamentals  of 
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radiometry,  will  be  discussed.  In  theory,  using  these 
expressions  one  can  calculate  the  radiance  from  any  given 
point  on  a  target . 

1.  Emitted  Radiance 

The  emitted  component  of  target  radiance  is  similar  to 
blackbody  radiation.  The  radiance  for  any  portion  of  the 
surface  is  given  by  Equation  3.3. 

LemiI(9r,  <pr)  =  -\  e(X,  9r,  <pr)Mhha,  T)dl  (3.3) 

71  \ 

The  term  e  (k,  <t>r,0r)  is  the  spectral  directional 

emissivity  at  the  target  intercept  point,  and  Mbb(^,T)  is 

the  spectral  radiant  emittance  at  target  temperature  T, 

whose  value  is  given  by  the  Planck  function.  The  angles  0r,i, 
and  <|>r,i  refer  to  zenith  and  azimuth  angles  in  reflected  and 
incident  directions  respectively.  The  integration  limits  Ai 

and  X.2  define  the  waveband  of  interest. 

If  the  emissivity  does  not  vary  with  wavelength 

Equation  3 . 3  becomes 

4 (3-4) 

n  i 
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2 .  Reflected  Radiance 

The  reflected  radiance  due  to  any  source  of  incident 
flux  may  be  computed  as  shown  for  any  desired  point  on  the 
target : 

^2 

Lref  <t>r)=  \\  ftt,  0r*  <t>r’  ftWft,  ft)  COS  OdQdX  (3.5) 

Xx  2k 

where  f  (X,  <t>r , ©r , <!>i , )  is  the  spectral  bi-directional 

reflectance  distribution  function  (BRDF)  ,  and  L  (X,  4>i,0i)  is  a 
specific  component (e.g. ,  target,  atmosphere,  sun,  etc.)  of 
the  input  radiance.  The  integration  is  performed  over  a 
hemisphere  defined  by  the  surface  element  under 
consideration,  and  over  the  wavelength  limits  defining  the 
waveband  of  interest.  [Ref.  13] 

Natural  sources  of  radiation  affect  target  signatures 
by  way  of  reflected  solar  radiation  and  earthshine  and  to  a 
much  lesser  extent  by  lunar  sources.  As  explained 
previously  the  reflected  solar  radiation  is  an  important 
factor  in  sea  background  signature  calculations  as  well  as 
target  signature  calculations.  Therefore,  in  this  project 
reflected  solar  contributions  are  involved  in  the  target 
signature  computation  model  and  in  the  background  radiation 
computation  model .  In  this  consideration  the  time  of  the 
day  when  the  measurements  are  made  becomes  important . 
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Reflected  earthshine  is  a  minor  component  of  the 
total  contrast  signature  for  nadir  viewing.  However,  as  the 
view  extends  toward  the  horizon  and  the  aspect  becomes  more 
side-on,  the  component  increases.  Since  in  this  project  the 
background  is  sea  background,  this  component  is  included  in 
the  "reflected  radiance  from  sea  surface"  part  of  the 
total  radiance. 

The  total  radiance  from  the  target  is  then 
computed  as : 


L^dr,t)  =  Lemit{6r,<Pr)  +  Lr/\dr^  (3.6) 


where  four  terms  refer  to  the  emitted,  reflected 
atmospheric,  reflected  direct  solar,  and  reflected  target 
radiances . 

The  radiant  intensity  is  found  by  integrating  the 
target  radiance  over  the  projected  target  area.  [Ref.  13] 

I(er,<pr)=  |  Llot  (0r ,  <pr )  cos  6rdA  (3'7) 

rarg  et 


3 .  Target  Types 

In  terms  of  the  radiances,  targets  can  be  considered 
in  three  different  groups:  ground  targets,  airborne  targets 
and  ship  targets . 
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Ground  targets  show  differences  in  emitted  radiation 
due  to  their  surface  emissivities  and  temperatures.  While 
some  targets  are  heated  with  respect  to  the  background  by 
internal  sources  such  as  the  power  train  in  a  tank,  some  of 
them  are  heated  by  natural  sources.  An  automobile  or  tank 
is  hotter  than  its  surroundings  for  several  hours  after 
use.  Personnel  are  usually  warmer  than  their  surroundings 
because  of  the  biological  processes  that  release  heat 
within  their  bodies.  Objects  that  do  not  release  heat 
internally  can  also  sometimes  be  detected  by  natural 
temperature  differences  between  the  object  and  the 
surroundings . 

Airborne  targets  such  as  airplanes  and  rockets  can  be 
detected  by  the  emission  from  the  heated  skin  and  from  the 
hot  motor  parts  and  also  by  the  radiation  from  the  plume  of 
hot  gases  from  the  motor.  [Ref.  8] 

A  ship  signature  is  made  up  of  two  main  components: 
Internally  generated,  and  externally  generated.  The  primary 
internal  IR  source  results  from  the  main  machinery  onboard 
any  vessel,  in  particular  drive  engines  and  electrical 
generators.  The  primary  external  sources  of  radiation  are; 
the  sun,  sky  radiance,  and  sea  radiance.  [Ref.  14] 
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In  the  case  of  a  ground  vehicle  or  ship,  the 
interaction  term  may  include  reflections  of  the  object  off 
its  environment  as  well  as  the  reverse.  For  these  cases, 
such  reflections  can  be  an  important  part  of  the  object's 
total  signature.  [Ref.  11] 

The  target  radiance  is  further  affected  by 
propagation  through  the  atmosphere.  To  determine  exactly 
the  radiance  originating  from  a  target  and  propagating  to 
the  sensor,  one  must  consider  the  atmospheric  transmission 
as  a  function  of  wavelength.  This  phenomenon  is  studied  in 
detail  in  Chapter  II. 

4.  Geometry  Effect 

In  the  general  case,  targets  are  complex  surfaces  and 
can  be  treated  by  dividing  the  surface  into  flat  facets  of 
constant  total  bidirectional  reflectance  and  then  summing 
over  all  facets  comprising  the  surface.  The  "faceted 
modeling"  technique  is  a  very  popular  method  for  describing 
the  geometry  of  targets  for  use  in  infrared  signature 
computations .  Many  infrared  simulation  models  in  current 
use,  including  the  one  which  is  used  in  this  project,  use 
faceted  representations  of  targets.  Figure  3.7  illustrates 
two  sample  targets  that  are  approximated  in  this  way.  Each 
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facet  is  assigned  thermal  and  optical  properties  and 
thickness  for  calculations  later. 

The  measured  value  of  radiance  is  affected  by  the 
location  of  the  observer  relative  to  the  target.  In  some 
cases  the  radiance  changes  significantly  as  the  aspect 
viewed  changes . 

Because  of  the  many  parameters  involved  in  computation 
of  radiance  from  a  target,  it  becomes  very  difficult  to 
establish  a  simple,  single,  rule  of  thumb  that  provides  an 
accurate  radiance  value  for  all  wavelengths  or  aspects. 

There  are  many  complex  Infrared  Signature  Simulation 
Software  Codes  in  the  electro-optic  technology  literature 
to  predict  the  infrared  radiance  as  well  as  the  zero  range 
and  apparent  temperatures  of  various  targets  for  different 
circumstances .  These  complex  models  are  mostly  based  on 
realistic  assumptions  about  the  optical  properties  and 
thermal  conditions  of  the  target  surface.  They  take  into 
account  quantitatively  the  potentially  significant  target 
radiation  resulting  from  reflected  earth,  sea,  skyshine, 
and  sunshine,  and  the  transmissivity  of  the  atmosphere 
between  the  target  and  the  observer.  Aspect -dependent 
signature  components,  such  as  variation  in  target  cross- 
sectional  area  and  weather  conditions  are  also  included  in 
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these  models.  Since  each  of  these  factors,  in  particular 
circumstances,  can  significantly  impact  the  infrared 
signature  of  an  object,  using  these  models  becomes  crucial 
for  obtaining  the  most  accurate  radiance  magnitudes.  Some 
of  these  models  including  short  descriptions  are  listed  in 
Appendix  B  from  Ref.  13. 


Figure  3.7  Typical  Samples  of  Faceted  Target  Models  Used  in 
Signature  Prediction  Programs.  After  Ref . [25] . 

The  computation  of  radiance  from  a  specific  target  in 
the  most  accurate  way  available  constitutes  one  of  the  main 
features  of  this  thesis.  In  the  following  Chapters,  the 
MuSES  (Multi-Service  Electro-Optic  Signature)  model  is  used 
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for  this  purpose.  In  Chapter  IV  detailed  parameters  and  the 
output  of  this  model  are  explained. 

C.  POLARIZATION  EFFECTS  ON  RADIANCE  &  TARGET 
DISCRIMINATION 

The  presence  of  interfering  background  radiation 
during  the  detection  process  is  a  very  common  problem  to 
most  infrared  sensors.  This  problem  degrades  the  sensors' 
target  detection  capability  and  limits  the  performance.  The 
use  of  polarized  infrared  energy  and  certain  selective 
polarization  properties  of  the  target  and  background  help 
to  distinguish  a  target  in  complex  natural  backgrounds  by 
improving  the  target  radiance  to  background  radiance  ratio . 
The  improvement  is  based  on  the  selective  rejecting  of  some 
part  of  the  polarized  radiation  from  backgrounds  or  target 
with  a  filter  at  the  entrance  of  the  sensor. 

While  researchers  have  recognized  the  potential 
usefulness  of  polarized  image  data  for  target  detection  and 
recognition,  wide-spread  use  of  polarization  in  the 
infrared  has  not  occurred.  This  is  probably  because 
polarization-sensitive  imaging  methods  have  been  slow  to 
develop,  and  require  detailed  studies. 

Polarization  sensitive  thermal  imaging  sensors  have 
many  potential  military  and  commercial  applications.  In 


particular,  stealth  ships  and  aircraft,  because  of  their 
large  planar  faceted  surfaces  are  expected  to  emit  highly 
polarized  thermal  radiation  which  reveals  the  shape  of  the 
object.  Therefore  polarization  characteristics  of  ships  are 
being  studied  for  use  in  detecting  ships  against  ocean  and 
sky  backgrounds.  [Ref.  15] 

Polarization  occurs  in  the  emitted  and  reflected 
infrared  radiation  of  targets  and  backgrounds .  Man-made 
objects  tend  to  be  glossier  and  more  regularly  shaped  than 
backgrounds  and  have  unnaturally  smooth  surfaces,  leading 
to  radiation  with  greater  polarization.  Although  natural 
backgrounds  such  as  grass,  trees,  dirt,  and  sand  emit  and 
reflect  radiation  that  is  less  polarized,  the  large  sea 
surface  near  grazing  incidence  shows  more  polarization  than 
the  more  detailed  faceted  man-made  ship  signature. 
Regarding  the  main  theme  of  this  thesis  some  important 
polarization  sources  are  summarized  below. 

1.  Sea  Emission  and  Reflection 

Background  infrared  radiation  from  the  sea  consists  of 
two  components;  the  partially  vertically  polarized 
emissions  from  the  sea  surface  and  the  partially 
horizontally  polarized  reflected  radiation  from  the  sky, 
sun,  moon,  etc.  Each  radiance  component  is  a  function  of 
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the  spectral  emissivity  of  seawater  and  the  orientation  of 
the  sea  wave  facet  relative  to  the  reflected  sources  and 
the  observer. 

Sea  surface  polarization  characteristics  have  complex 
mathematical  models.  Reference  20  presents  detailed 
descriptions  of  some  models  frequently  used  for  analysis  of 
the  polarization  characteristics  of  the  sea  surface. 
Previous  experiments  and  analyses  [Ref.  1  and  Ref.  16]  have 
shown  that  the  electric  vector  of  sea  surface  radiance 
appears  horizontally  polarized  within  the  sun  glint  regions 
and  vertically  polarized  otherwise.  Figure  3.8  shows  the 
predicted  degree  of  polarization  for  MWIR  and  LWIR  in  and 
near  the  sun  glint  region  as  a  function  of  solar  zenith 
angle  for  an  observer  looking  in  the  vertical  plane  of  the 
sun.  The  negative  values  of  polarization  outside  of  the  sun 
glint  regions  indicate  vertical  polarization.  In  these 
previous  experiments  measured  LWIR  surface  radiance  shows  a 
degree  of  polarization  with  vertical  orientation  around  7 
to  30%. 

In  this  project  the  degree  of  polarization  of  a  sea 
background  radiance  for  the  specified  geometry  and 
conditions  is  predicted  by  the  polarized  version  of  the 
SEARAD  Radiance  Model . 
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Figure  3.8  Theoretical  Degree  of  Polarization  as  Function 
of  Azimuth  Angle  and  Solar  Zenith  Angle,  for  MWIR  and  LWIR. 
From  Ref . [1] . 


2.  Sky  Emission  and  Reflection 

Atmospheric  path  radiance  will  not  contribute  to  a 
change  in  apparent  target  to  background  radiance  contrast 
ratio  seen  with  a  polarized  sensor.  The  only  absorption 
mechanism  for  path  radiance  will  be  the  transmissivity  of 
the  polarized  filter.  In  this  project,  atmospheric  path 
radiance  is  included  into  the  total  target  and  total 
background  radiance  values  separately,  and  in  each  case 
treated  as  having  0%  degree  of  polarization. 
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Due  to  the  high  concentration  of  water  in  the 
atmosphere,  the  sky  is  generally  an  excellent  absorber. 
Since  it  is  an  excellent  absorber,  it  is  a  poor  reflector. 

3 .  Target  Emission  and  Reflection 

As  sea  surface  radiance  shows  a  net  vertical 
polarization  characteristic,  primary  contrast  improvement 
can  be  obtained  by  placing  a  horizontal  filter  in  a  sensor. 
In  cases  where  the  polarization  due  to  emission  from  the 
target  has  a  perpendicular  direction  compared  to  the 
background  polarization,  the  use  of  a  polarized  filter  will 
cause  an  additional  contrast  improvement.  In  previous  work 
[Ref.l]  it  was  concluded  that  polarization  features  were 
not  prominent  in  the  images  of  ship  targets,  in  either  the 
MWIR  or  LWIR  bands,  with  no  observations  greater  than  ±5%  in 
degree  of  polarization.  A  parallel  conclusion  is  stated  in 
Ref.  19  also,  where  it  is  pointed  out  that  the  degree  of 
polarization  from  target  surfaces  in  the  IR  is  typically  5 
to  8  percent  under  appropriate  viewing  conditions.  Although 
5  to  8  percent  degree  of  polarization  seems  to  be  a  small 
value,  in  suitable  conditions  it  can  create  a  significant 
temperature  difference  and  consequently  a  significant 
signal  for  detection. 
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Painted  surfaces  such  as  found  on  all  types  of 
vehicles  and  ships  display  emission  polarization 
particularly  at  angles  close  to  the  Brewster  angle.  Figure 
3.9  shows  21  measurements  taken  on  various  paint  samples 
ranging  from  very  smooth  texture  through  sand-paint 
mixtures .  The  amount  of  polarization  depends  on  the  index 
of  refraction  of  the  paint  and  any  degradation  due  to  the 
surface  roughness . 
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Figure  3.9  Measured  Polarized  Signature  Components  in 
the  7.5-12(jm  band  for  Selected  Paint  Samples.  After 
Ref. [12] . 


The  index  of  refraction  of  the  surface  material  of  the 
target  determines  the  magnitude  of  the  polarized  components 
of  total  radiance  and  the  net  polarization  will  vary 
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according  to  the  balance  between  thermal  emission  and 
reflection  in  the  total  radiance. 

In  this  project  the  performance  of  a  polarized  second 
generation  sensor  was  studied  for  a  target  degree  of 
polarization  limited  to  0  to  8  percent  in  both  directions. 
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IV.  MODEL  DEFINITION  AND  TARGET  SIGNATURE  COMPUTATION 

A.  MODELING  APPROACH 

1 .  Background 

Maximum  range  analysis  for  thermal  imaging  systems  with 
and  without  polarizer  was  studied  in  previous  projects,  by 
Yu  (Ref.  27,  1997),  Lagaras  (Ref. 18,  1999),  and  by 

Guimaraes  (Ref.  23,  1999).  While  Yu  used  a  scenario 

selected  from  the  database  of  the  EOMET95  measurements 
series  in  Monterey  Bay  off  Moss  Landing,  the  analyses 
carried  out  by  Lagaras  and  Guimaraes  were  based  on  the 
EOPACE  data  acquired  in  April  1996  in  San  Diego  Bay.  In 
Ref.  27,  maximum  ranges  were  calculated  without 

polarization  filtering.  In  Ref.  18,  the  calculations  were 
performed  for  an  unpolarized  target  considering  the 
polarization  filter  effect,  since  no  target  polarization 
was  observed  in  the  MAPTIP  or  EOPACE  data.  The  effect  of  a 
polarization  filter  on  the  Minimum  Resolvable  Temperature 
Difference  of  a  sensor  is  studied  in  Ref.  23.  In  all  of  the 
above  target-to-background  models  the  common  inadequacy  is 
lack  of  a  valid  target  signature  computation  model.  In  the 
models  used,  by  Guimaraes  and  Lagaras  the  ship  targets  are 
assumed  to  be  unpolarized  and  the  sea  background  degree  of 
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polarization  is  averaged  within  the  required  bandwidth.  In 
Ref. 18  it  is  observed  that  when  a  horizontally  polarized 
filter  is  used,  at  some  distance  the  apparent  temperature 
difference  passes  through  zero  and  goes  negative,  whereas 
it  is  expected  that  at  greater  ranges  the  apparent 
temperatures  of  the  target  and  background  should  approach 
each  other  asymptotically,  so  that  the  apparent  temperature 
difference  should  approach  zero.  This  inconstancy  in 
apparent  temperature  difference  found  in  Ref.  18  is 
interpreted  as  probably  resulting  from  the  use  of  different 
incompatible  procedures  for  target  and  background  apparent 
temperature  calculations.  In  this  project  an  algorithm  is 
developed  for  apparent  temperature  difference  calculation 
to  overcome  the  anomaly. 

In  this  thesis  a  model  that  contains  a  commercial 
advanced  target  signature  prediction  program,  MuSES  (Multi- 
Service  Electro-optic  Signature) [Ref .  25],  is  developed  to 
predict  the  maximum  range  of  a  sensor  for  specific  tasks. 
The  model  includes  a  ship  target  and  sea  background,  each 
having  a  degree  of  polarization.  The  scenario  and  the  model 
are  defined  in  the  following  sections. 

2 .  Scenario 

For  the  model  a  single  scenario  is  designed.  The  target 
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in  the  scenario  is  a  modified  version  of  a  ship  target 
contained  in  the  MuSES  program.  The  ship  target  in  the 
program  is  a  faceted  model  of  a  210  foot  ship  (dimensions: 
64.8LX7.14WX7.93H  (m)  )  which  is  similar  to  the  Cyclone 

class  ships  (170  foot)  .  The  geometry  and  ship  dimensions 
are  detailed  in  Appendix  C. 

The  background  in  the  scenario  is  a  typical  sea 
background  with  an  average  polarization  value  -13-6  in  the 
specified  bandwidth.  The  overall  analysis  is  performed  in 

the  8-12|im  region. 

The  sensor  simulated  has  the  generic  characteristics 
of  a  Second  Generation  FLIR  sensor.  The  characteristics  are 
given  in  Appendix  F.  For  the  scenario  it  is  assumed  that 
the  hypothetical  sensor  is  mounted  on  an  aircraft  which 
approaches  the  target  at  300m  elevation.  Figure  4.1  shows  a 
schematic  representation  of  the  scenario. 


◄ - 

A  Ship  Target  on  Sea 
Background 


Sensor  Mounted  on  Aircraft 


300  m 


Figure  4.1  Schematic  Representation  of  the  Scenario. 
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Sea  radiance  and  atmospheric  propagation  are  modeled 
with  the  SEARAD  propagation  code.  For  the  atmospheric  data 
required  as  input  to  SEARAD  and  to  the  target  model  in 
MuSES  one  of  the  weather  data  files  available  in  the  MuSES 
program  database  is  used.  The  Natural  Environment  (weather) 
option  in  the  program  uses  models  or  measurements  of  the 
weather  conditions  and  the  solar  position  to  calculate  the 
cooling  or  heating  that  occurs  from  the  influence  of 
weather.  Currently  there  are  four  different  weather  files 
available  in  the  program.  For  the  target  signature 
calculations  in  this  project  the  file  dated  Jul  20,1984  is 
used.  The  parameters  are  as  follows: 

Global  Positioning: 

Latitude  (deg.):  47.2,  Longitude  (deg.):  88.5,  Time 

Zone  (hrs . ) :  4,  Elevation  (m) :  332 

Date  of  the  Measurements:  07/19/1984 

Available  Solution  Times :  From  07/19/1984  06:20  AM  to 
07/20/1984  08:59  AM 

Selected  Solution  Time:  07/19/1984  18:00  PM 

This  weather  file  is  selected  as  being  compatible  with 
a  standard  MODTRAN  atmosphere  (Midlatitude  Summer) . 
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3 .  Model  Definition 

The  model  developed  in  this  project  contains  a  defined 
scenario  and  a  series  of  calculations  that  are  carried  out 
in  this  and  in  following  Chapters.  It  can  be  summarized  in 
following  four  steps: 

•  Target,  Background,  and  Path  Radiance  Calculations 

•  Minimum  Resolvable  Temperature  Difference  (MRTD) 
Prediction 

•  Apparent  Temperature  Difference  (ATD)  Calculation 

•  Maximum  Range  Predictions 

The  steps,  more  detailed,  are  shown  in  Figure  4.2.  As 
opposed  to  the  previous  models,  in  this  model  an  advanced 
target  signature  prediction  program  is  used  and  both  target 
and  sea  background  polarization  effects  are  included. 
Target  and  background  computations  are  considered  to  be 
compatible.  The  program  description  and  computations  are 
presented  in  the  following  section. 

B.  TARGET  SIGNATURE  COMPUTATION 

The  performance  of  many  infrared  seekers  and  other 
sensor  systems  operating  in  the  far  infrared  spectral 
region  of  8|im  -12  |im  is  critically  dependent  on  the 
magnitude  of  thermal  signatures  (i.e.,  thermal  contrasts, 
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Polarization 

Filter 

Included 


Apparent  Background 
Radiance  Calculation 


Specific 
Task  and 
Target 
Description 


Target  Degree  of  Pol. 
Polarized  Ap.Rad.Calc. 


Background  Degree  of  Pol.Cal. 
Polarized  App.  Background  Rad.Cal. 


MRTD  (R) 
Calculation 


Conversion  of  Radiance  into  Temp. 
Planck's  Formula 
Look-up  Table 

Tjargcb  T Background 


ATD  Calculation 

ATD(R)—  T Tar  get)  ~  T Background 


Maximum  Range 
Calculations 


Figure  4.2  Flow  Chart  of  The  Model  Used  in  This  Project. 


or  equivalently,  differences  in  thermal  radiance  values 
between  various  scene  components) . [Ref.  11] 

To  detect  a  target  in  a  background,  a  contrast  in 
radiance  alone  is  not  enough,  but  is,  however  essential. 
Since  signature  description  is  a  primary  part  of  infrared 
systems  analysis,  signature  calculation  models  play  an 
important  role  in  performance  prediction  of  infrared 
systems . 

These  models  should  include  and  reproduce  the  effects 
of  the  most  important  radiative  signature  processes. 
Because  infrared  systems  typically  predict  contrast  between 
a  target  and  its  background,  and  because  natural  sources  of 
target  infrared  radiation  are  frequently  important, 
environmental  as  well  as  target  radiative  phenomena  must  be 
included  in  predictive  analysis  tools . 

To  achieve  reasonably  accurate  estimates  of  target 
signatures  it  is  most  important  to  understand  the  complex 
phenomenology  that  determines  the  levels  of  radiation.  The 
phenomenological  processes  are  particularly  dependent  on 
the  scenario  considered. [Ref .  11] 

In  the  rest  of  this  Chapter  the  MuSES  (Multi-Service 
Electro-optic  Signature)  infrared  signature  prediction 
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program  is  used  in  order  to  obtain  accurate  radiance 
predictions  of  ship  targets. 

C.  MUSES  INFRARED  SIGNATURE  PREDICTION  PROGRAM 

MuSES  [Ref.  25]  is  a  new  commercial  "next  generation" 
infrared  signature  prediction  program  developed  by 
ThermoAnalytics  Inc.  which  lets  the  user  do  complete 
thermal  modeling  and  infrared  signature  prediction  within 
an  integrated  easy-to-use  interface.  This  software  is  the 
replacement  for  PRISM  (The  Physically  Reasonable  Infrared 
Signature  Model),  the  previous  standard  US  Army  infrared 
prediction  code  and  the  TCMII  model  which  are  both  used  for 
ground  vehicles  and  other  targets  of  interest.  With  MuSES 
one  can  create  one's  own  target  model  or  use  one  of  the 
standard  models  in  the  program.  It  has  also  a  capability  of 
editing  and  manipulating  the  existing  model  geometry  e.g., 
to  scale,  translate,  rotate,  copy  and  delete.  This  model 
includes  the  detailed  environmental  radiation  effects  into 
the  radiance  and  temperature  computations .  Modifying  the 
characteristics  of  target  and  background,  such  as 
emissivity,  conductivity,  density  and  specific  heat. 


various  scenarios  can  be  analyzed. 


D.  COMPUTATIONS  WITH  MUSES 


In  this  project  a  ship  thermal  model  is  used  which  is 
modified  from  the  original  MuSES  sea  vehicle  thermal  model. 
The  target  dimensions,  the  other  important  characteristics 
of  the  target,  and  the  geometry  of  the  calculation  are 
presented  in  Appendix  C .  The  target  model  is  used  with  the 
weather  file  (included  in  MuSES)  which  is  shown  in  Appendix 
D.  In  the  model  scenario  chosen  the  target  is  heading  in 
the  direction  of  north  so  the  simulated  measurement  made  in 
the  late  afternoon  was  not  affected  by  sun  glint.  The 
result  of  these  calculations  as  well  as  the  time  varying 
temperature  and  radiance  measurements  are  presented  in 
Appendix  E.  These  results  are  then  used  in  overall 
performance  analysis  of  the  thermal  imaging  sensor  in 
Chapter  V. 
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V.  MRTD-ATD  COMPUTATIONS  AND  ANALYSIS 


A.  MINIMUM  RESOLVABLE /DETECTABLE  TEMPERATURE  DIFFERENCE 

The  overall  performance  of  a  thermal  imaging  system  is 
governed  by  both  its  thermal  sensitivity  and  its  spatial 
resolution,  and  is  generally  defined  in  terms  of  its 
Minimum  Resolvable  Temperature  Difference  (MRTD)  and 
Minimum  Detectable  Temperature  Difference  (MDTD) . 

MRTD  has  become  a  standardized  test  methodology  for 
thermal  imaging  systems.  It  measures  the  temperature 
difference  at  which  an  observer  is  just  able  to  discern  a 
four  bar  (7:1  aspect  ratio)  chart  against  a  background 
thermal  scene.  MRTD  gives  the  detectable  temperature 
difference  as  a  function  of  the  spatial  frequency  of  the 

bars . 

The  parameter  MRTD  was  originally  defined  for  a 
scanning  FLIR,  and  in  the  measurement  procedure  the  bars  of 
the  four-bar  pattern  were  oriented  perpendicular  to  the 
scan  direction  of  the  FLIR.  The  observer  is  allowed  to 
adjust  the  gain  and  level  settings  of  the  FLIR  as  well  as 
his  or  her  position  relative  to  the  system  display  during 
the  measurement  in  order  to  minimize  the  temperature 
difference  required  to  just  resolve  the  target. 
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MDTD  is  the  temperature  difference  between  a  square 
(or  circular)  target  and  a  uniform  background  when  viewed 
through  a  FLIR,  which  is  required  by  a  observer  to  just 

detect  the  target.  [Ref.  21] 

There  are  various  approaches  used  to  calculate  MRTD 
and  MDTD  values  in  the  infrared  community.  One  of  these 
formulas, and  maybe  the  most  common  one,  is  given  in  Ref. 22: 


MRTD(f)  = 


20 •  SNRt  ■  (FOVx  •  FOVjl2f-p 


1/2 

jc 


ra-D-D*  {jt  n- tjsc )1/2 •  {ocfi)V2 


MTF 

LYXXJ-  system 


(rrj 


in  cW 

ar 


(5.1) 


where : 

SNRt  =  perceived  signal-to-noise  ratio  set  according  to  the 
desired  probability  of  performing  the  given  task 
f  =  spatial  frequency  ( cycles /mrad) 

D  =  aperture  diameter  (m) 

D*  =  detectivity  (cm  Hz%  W'1) 

n  =  number  of  detectors 

T|sc  =  scan  efficiency 

T0  =  transmission  of  the  optics 

(X  =  in-scan  detector  angular  subtense  (mrad) 

(3  =  cross— scan  detector  angular  subtense  (mrad) 

MTF  system  =  the  system  modulation  transfer  function 


L  =  the  length- to-width  ratio  for  the  bar  chart  (7) 


Te  =  eye  integration  time  (0.2s) 

—  derivative  of  Planck’s  Law  (watt  cm  K  sr  ) 

dT 

FOVx  =  in-scan  field  of  view  (mrad) 

FOVy  =  cross-scan  field  of  view  (mrad) 

px  =  noise  filter  factor 

B .  MRTD  COMPUTATION 

In  this  project  a  hypothetical  second  generation  FLIR 
sensor  developed  in  Ref. 14  is  used  for  performance 
modeling.  The  parameters  of  the  sensor  are  presented  in 
Appendix  F  and  Figure  5 . 1  shows  the  plot  of  MRTD  values  of 
the  sensor  versus  spatial  frequency. 


MRTD  &  Spatial  Frequency 
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The  inclusion  of  a  polarized  filter  into  the  system  may 
alter  the  effect  of  some  parameters  in  MRTD  computation. 
Guimaraes  (Ref.  23)  has  reported  an  experimental  evaluation 
of  this  topic  and  concluded  that  the  polarizer 
transmissivity  (Tp) ,  polarizer  Modulation  Transfer  Function 
(MTFp)  ,  and  the  partial  derivative  of  Planck's  Equation  with 
respect  to  the  temperature  are  the  parameters  which  have 
the  greatest  impact  in  the  reformulation  of  MRTD  for  the 
polarized  case.  The  MRTD  values  of  the  sensor  which  is  used 
in  this  project,  were  computed  using  the  3D  Noise  concept. 
[Ref.  24]  Therefore  the  functional  relation  between  the 
MRTD  values  and  the  sensor  parameters  is  not  well  defined. 
To  identify  the  variation  of  the  MRTD  values  due  to  the 
polarization  effect,  the  algorithm  of  the  calculation  has 
to  be  known.  Bearing  in  mind  that  next  generation  FLIR 
sensor  performance  models  may  not  use  the  3D  noise  concept, 
and  that  the  impact  of  changing  the  other  parameters  in 
Equation  5.1  will  not  effect  the  overall  conclusion,  only 
Tp/  the  polarizer  transmissivity  value,  was  taken  into 
account  in  the  "polarization  filter  included  MRTD 
calculation" . 
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At  this  point,  in  agreement  with  the  conclusions  of 
Guimaraes  (Ref. 23)  the  MRTD  values  were  divided  by  0.85, 
which  is  manufacturer's  value  for  the  filter  used  in  the 
work.  It  must  be  recognized  that  the  process  used  here  and 
in  Ref. 23  i.e.,  dividing  MRTD  values  by  a  transmissivity, 

is  not  a  completely  appropriate  solution  to  the 
polarization-included  MRTD  calculation  issue,  but  that  it 
is  an  approximation  that  needs  further  discussion.  While 
the  transmittance  can  be  applied  to  radiant  flux  or 

irradiance,  it  should  not  be  applied  directly  to  the 

temperature  difference  term.  Appendix  G  shows  the  MRTD 

values  of  the  sensor  with  and  without  the  polarization 
filter.  In  the  rest  of  the  project  the  MRTD  values  and  MRTD 
curves  refer  to  the  case  in  which  the  polarizer  is  added  to 
the  system. 

C.  DETECTION  CRITERIA  AND  PROJECTED  AREA  OF  THE  TARGET 

In  order  to  compare  the  sensor's  maximum  range 

performance  with  and  without  the  polarization  filter,  we 

need  to  express  MRTD  values  in  terms  of  range.  So  far  all 

the  MRTD  calculations  have  been  made  as  a  function  of 

spatial  frequency.  In  performance  test  studies  if  the 

specific  task  is  given  and  target  dimensions  are  known,  the 

spatial  frequency  can  be  converted  into  distance  with  the 
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help  of  some  established  criteria,  (e.g*  Johnson  Criteria) 
which  relate  the  resolution  of  bar  charts  to  levels  of 
visual  discrimination. 

The  first  experimental  determination  of  these  criteria 
was  conducted  by  John  Johnson  in  1957.  Later  Moser  and 
O'Neil  made  separate  studies  relating  especially  to  ship 
targets.  Although  the  approaches  show  slight  differences, 
all  the  criteria  address  the  number  of  resolution  elements 
per  critical  dimension  required  by  a  typical  experienced 
operator  for  the  specified  task.  The  Johnson  criterion 
provides  a  relation  between  the  target  angular  subtense 
(TAS)  and  the  spatial  frequency  scale  on  the  MRTD  graph. 
The  number  of  cycles  that  appear  across  the  critical  target 
dimension,  which  can  just  be  resolved  by  the  sensor  at 
range  R,  can  be  found  by  the  relation 

N=f . TAS  (5.2) 

where  TAS=  Critical  Target  Dimension  (Dc)  /Range  and  f  is  the 
spatial  frequency.  Figure  5.2  shows  the  relation  between 
parameters  and  Table  5.1  shows  the  current  discrimination 
criteria  with  50%  probability  for  the  one  and  two 
dimensional  MRTDs .  The  same  approach  here  to  find  the  MRTD 
as  a  function  of  range  has  been  used  in  some  of  the 
previous  works  related  to  the  subject. 
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Dc 

◄ - ► 


For  2  resels  =1  cycle  on  Critical  Dimension 

Spatial  Frequency,  f  =  1/TAS  where  TAS=  Dc/R 

For  N  cycles  on  Critical  Dimension  f=  N  Dc/R  (cycles/mrad) 


Figure  5.2  Conversion  of  Spatial  Frequency  into  Range. 
After  Ref. [23]. 


Discrimination 

level 

Meaning 

One-D  cycles 
across 

critical 

dimension 

Two-D  cycles 
across  critical 
dimension,  Ns0 

Detection 

An  object  is 
present 

1 

0.75 

Recognition 

Class  to  which  an 
object  belongs 
( e .  g . ,  human , 
tank,  truck) 

4 

3 

Identification ! 

Object  is 
discerned  with 
enough  clarity  to 
specify  type  (MIA, 
T-62,  T-72  Tank) 

8 

6 

Table  5.1  Resolution  Cycles  on  Target  Required  for 
Specified  Levels  of  Recognition  for  Johnson  Criteria.  From 

Ref. [24]. 
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The  second  step  in  converting  spatial  frequency  into 

distance  is  to  calculate  the  "critical  dimension",  Dc.  This 

dimension  is  the  average  target  dimension  being  observed 
for  detection  or  recognition  and  simply  defined  as  the 
square  root  of  the  target-projected  area. 

In  this  project  a  ship  model  created  in  the  MuSES 

environment  is  used  as  a  target.  Its  dimensions  and  the 
orientation  of  the  target  when  the  radiance  measurements 

were  made  are  displayed  in  Appendix  C. 

The  same  orientation  is  used  for  the  calculation  of 
projected  area.  As  expected  the  aspect  angle  plays  a  major 
role  in  the  calculation.  The  projected  area  of  the  target 
as  a  function  of  its  dimensions  and  the  azimuth  and 

elevation  angle  is  found  by: 

(5.3) 

At  -  lh  cos  6  cos  <p  +  wh  cos  6  sin  <p  +  lw  sin  6 
where : 

1=  ship  length 
w=ship  width 
h- ship  height 
6=  elevation  angle 
(f)=  azimuth  angle 
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MRTD  Values 


Using  Equation  5.3,  for  the  scenario  conditions, 
projected  area  of  target  is  found  to  be  610.77m2  and  the 
critical  dimension,  Dc,  is  found  to  be  24.71m.  From  these 


parameters  and  using  Equation  5.2,  the  MRTD  values  of  the 
sensor  were  computed  as  a  function  of  range,  using  the 
values  N=0.75  for  detection,  N=3  for  recognition  and  N=6 
for  identification.  The  tabulated  values  are  presented  in 
Appendix  G.  Figure  5.3  shows  the  MRTD  for  detection, 
recognition,  and  identification  respectively. 


MRTD  (N=0.75,  N=3,  N=6)  &  RANGE 


— •—  Ns=6,  Identification  Task 
— a—  N=3,  Recognition  Task 
; — A —  N=0.75,  Detection  Task 


Figure  5.3  Minimum  Resolvable  Temperature  As  a  Function  Of 
Range,  For  Detection  (N=0.75),  Recognition  (N=3),  and 
Identification  (N=6) .  "N"  is  Number  of  Resolved  Cycles  on 
The  Critical  Dimension. 
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These  results  were  used  in  the  following  sections  to 
find  the  maximum  detection,  recognition  and  identification 
range  of  the  polarized  sensor  for  the  given  weather 
conditions . 

D.  APPARENT  TEMPERATURE  DIFFERENCE  COMPUTATION 

Apparent  differential  temperatures  are  equivalent 
blackbody  differential  temperatures  that  represent  a 
target- to-background  contrast  with  degradation  due  to 
atmospheric  effects.  Apparent  differential  temperature  or 
more  commonly  "apparent  delta  T"  depends  on  a  number  of 
parameters  including  source  emission  characteristics  and 
atmospheric  transmission  characteristics.  There  are  many 
different  techniques  for  calculating  apparent  delta  T. 

In  this  project  the  "radiance  computation  method”  is 
used  for  calculation  of  apparent  delta  T.  The  method  can  be 
summarized  by  the  following  seven  steps: 

•  Calculate  the  total  radiance  from  the  target  using 
the  MuSES  signature  computation  model  at  zero  range 
for  the  specific  weather  conditions  and  target 
orientation.  The  weather  file  and  the  result  file 
for  this  computation  are  presented  in  Appendix  D  and 
E  respectively.  The  target  radiance  calculations 
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were  made  as  explained  in  the  last  section  of 
Chapter  IV. 

•  Calculate  the  background  radiance,  including 
emission  and  reflection  components,  from  the  sea 
background,  using  the  SEARAD  Radiation  Model  for 
different  zenith  angles  that  correspond  to  different 
ranges . 

•  Correct  the  "zero  range  total  target  radiance" 

component  of  ATD  for  atmospheric  attenuation  for  a 
set  of  different  ranges  that  are  appropriate  to  the 
data  set.  At  this  point  the  wavelength  interval  is 
divided  into  four  sub-bands,  then  the  atmospheric 
transmittance  values  were  computed  corresponding  to 
each  sub-band  and  the  inband  radiance  values  were 
multiplied  by  those  transmittance  values.  By  doing 
so  we  have  reduced  the  error  that  comes  from 
averaging  transmittance  values  in  broadband. 

•  Correct  the  background  and  target  radiances  for  path 
radiance,  for  the  same  set  of  ranges  and  the 
transmission  coefficients,  computed  by  using  SEARAD 
propagation  code.  The  results  are  presented  in 
Appendix  I . 
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•  Calculate  corrected  total  apparent  target  radiance 
by  adding  the  path  radiance  and  target  radiance 
values  together.  This  step  can  be  expressed  in  an 
equation  as: 

Lapp. Target ( total ) =  ^Target  *  ^  (Range)  +  Lpath  Radiance  (5.4) 

•  Convert  these  radiances ,  i .  e  " total  ship  radiance 
plus  path  radiance"  and  "background  radiance"  values 
were  converted  into  temperature  values  by  using 


Planck's  formula  (Equation  2.6). 

Lapp. Target  (Range)  - ►  TApp. Target  (Range) 

Lapp  •  Background  ( Rang  e )  - ►  Tapp  .  Background  ( Range ) 

For  this  conversion  a  short  MATLAB  code  radconv.m  and 
a  look-up  table,  which  was  created  with  the  result  of  this 
code,  were  used.  The  code  and  the  look— up  table  are 


presented  in  Appendix  H. 

•  Subtract  the  temperature  values  found  from 
conversion  of  the  full  range  background  radiance 
from  the  temperature  values  found  from  conversion  of 
total  ship  plus  path  radiance,  to  find  the  apparent 
temperature  dif f erence (ATD)  for  the  specific  range. 

A  schematic  representation  of  the  ATD  calculation  process 


is  given  in  Figure  5.4. 


ATD  (R)—  TApp.Target  -  TApp.Background 


Sensor  Mounted  on  Aircraft 


TTarget»T  Background 

At  Ranee  “R” 


From  the  Target 

A  Ship  Target  on  Sea  Background 


Figure  5.4  A  Schematic  Representation  of  ATD  Calculation. 

Appendix  I  shows  the  computed  background  radiance, 
target  radiance  values  and  ATD  values  for  the  corresponding 
range  values.  The  ATD  values  were  then  plotted  versus  the 
range  values  as  seen  in  Figure  5.5.  From  the  plot  it  is 
observed  that  the  curve  fluctuates  in  the  short  range 
region  (1-5  km)  .  This  behavior  is  further  discussed  in 
Chapter  VI,  Conclusions  and  Recommendations. 

In  Appendix  K  the  radiance  difference  values  and 
corresponding  temperature  difference  values  are  tabulated. 
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By  using  a  simple  Microsoft  Excel  regression  model  a 
relation  which  is  almost  linear  and  independent  of  the 
range  was  found.  This  result  simply  states  that  both 
radiance  difference  and  temperature  difference  behave  m 
the  same  manner.  The  different  aspects  of  the  study  are 
considered  under  discussion  Chapter  VI. 


ATD&Range 


where  Sensor  Altitude  is  300m. 


E.  THE  EFFECT  OF  DEGREE  OF  POLARIZATION 

The  observed  degree  of  polarization  of  the  radiance 
can  be  defined  as: 


POL  = 


(w.)-(w,) 

<*.>+(*,> 


(5.5) 
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where : 


<NV>  =  mean  apparent  vertically  polarized  radiance. 

<Nh>  =  mean  apparent  horizontally  polarized  radiance. 
According  to  this  definition  a  degree  of  polarization 
greater  than  zero  implies  a  net  vertical  polarization  while 
the  degree  of  polarization  value  less  than  zero  implies 
horizontal  polarization. 

In  general  the  total  radiance  from  any  object  has  a 
vertical  and  horizontal  component  frequently  expressed  as  a 
function  of  the  total .  While  the  addition  of  these 
normalized  components  is  always  unity,  the  difference 
between  them  is  defined  as  the  degree  of  polarization. 
These  two  expressions  are  shown  in  Equations  5.6  and  5.7. 


_J^> _ (hl_=P0L 

(N,)+(Nh)  (N,)  +  (N„) 

(N  )  (Nh) 


(5.6) 

(5.7) 


Using  these  two  equations  we  can  express  the  fractional 
vertical  and  horizontal  polarization  radiance  in  terms  of 


the  degree  of  polarization  as: 

(Nv)  _  (1  +  POL) 

(Nv)  +  <Nh)  2 

K)  =  (1  -  POL) 

(Nv)  +  (Nh>  2 


(5.8) 


(5.9) 
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It  can  be  shown  from  the  above  equations  that  with  a 
zero  degree  of  polarization,  i.e. ,  in  an  unpolarized  case, 
half  of  the  total  radiance  shows  vertical  polarization  and 
the  other  half  shows  horizontal  polarization.  This  notation 
is  found  in  previous  publications. 

Defining  the  effect  of  the  degree  of  polarization  on 
radiance  components  in  this  way,  we  can  now  compute  the 
polarized  target  radiance  and  polarized  background  radiance 
that  enter  the  sensor  when  a  polarizer  is  added  to  the 
sensor.  As  mentioned  before,  in  performance  analysis 
knowledge  of  the  polarization  characteristics  of 
backgrounds  and  targets  is  very  important.  The  contrast 
difference  between  target  and  the  background  can  be 
improved  by  using  an  appropriate  polarizing  filter.  In  some 
previous  works  (e.g. ,  Ref.  20) ,  it  has  been  showed  that  the 
radiation  from  the  sea  surface  shows  vertical  polarization 
in  the  LWIR.  Knowing  this  fact  we  can  use  a  horizontally 
polarized  filter  to  eliminate  some  of  the  background 
radiation  and  obtain  better  contrast  difference  depending 
on  the  degree  of  polarization  from  a  target. 

In  this  project  the  first  goal  is  to  find  the  effect 
of  degree  of  polarization  on  the  apparent  temperature 
difference  and  investigate  whether  this  degree  of 


76 


polarization  makes  any  improvement  in  the  maximum  ranges  of 
a  sensor  when  a  polarized  filter  is  added.  In  order  to 
achieve  this  goal,  the  computations  were  made  for  the  cases 
where  sea  background  shows  vertical  polarization  with 
values  between  12%  -  14%.  The  exact  polarization  values 

were  computed  from  the  SEARAD  Radiance  Model  and  included 
into  the  calculations.  In  order  to  find  the  effect  of 
target  polarization  in  range  performance  the  cases,  in 

which  the  target  has  a  polarization  with  values  between  0%- 
8%,  both  vertically  and  horizontally,  were  studied.  As  an 
example,  when  a  perfect  horizontally  polarizing  filter  is 
used  and  the  target  has  a  2%  horizontal  degree  of 

polarization,  51%  of  the  target  radiance  would  pass  the 
filter  while,  for  an  average  15%  vertical  degree  of 

polarization  of  sea  surface  57.5%  of  the  background 

radiance  would  be  eliminated. 

Appendix  J  shows  the  radiance  values  from  target  and 
from  background  when  the  horizontal  polarizing  filter  is 
added  and  the  corresponding  ATD  values  for  this 
configuration.  Figure  5.6  shows  the  ATD  curves  for 
polarized  cases  versus  range  for  target  polarization  from  0 
to  8%  (horizontal) . 
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ATD &  Range 


-A7D  ( Target  0%  Dpol.,  with 
Polarized  Filter) 


-ATD  (Target  2%  Dpol,  with 
Polarized  Filter) 


— *r-A7D  (Target  8%  Dpol.,  with 
Polarized  Filter) 


-ATD  (No  Polarized  Filter;  with 
Unpolarized  Target,  Unpolarize 
Sea  Background ) 


Range  (km) 


Figure  5.6  Apparent  Temperature  Difference  Values  For 
Polarized  and  Unpolarized  Cases  versus  Range;  Target 
Polarization  as  Parameter. 


F.  RANGE  CALCULATIONS 

The  maximum  range  of  a  sensor  for  detection, 
recognition  or  identification  can  be  estimated  using  the 
Apparent  Temperature  Difference  (ATD)  and  the  MRTD-MDTD 
curves.  The  condition  for  detection  at  maximum  range  is 
that  the  ATD  will  match  the  required  MRTD  or  MDTD  at  the 
same  range. 

As  a  last  step,  the  ATD  and  the  MRTD  curves  for 
different  tasks  are  plotted  on  the  same  scale  as  a  function 
of  range  for  the  selected  target  with  a  critical  dimension 
of  24.71m.  In  each  case  the  intersection  of  corresponding 
ATD  and  MRTD  curves  gives  the  maximum  range  required  for 
that  specific  task.  Figures  5.7  to  5.12  show  the  different 
cases . 

G.  DISCUSSION  OF  THE  DATA 

Figure  5.6  shows  the  expected  decrease  in  ATD  with 
range,  all  curves  approaching  zero  in  the  long  range 
limit, as  opposed  to  Ref. 18.  As  expected,  increase  in 
horizontal  polarization  of  target  leads  to  increase  of  the 
ATD,  while  removing  the  filter  (transmitting  the  polarized 
sea  radiance)  leads  to  a  decrease.  These  results  are 
consistent  with  previous  experiences  in  Refs.  18,  23,  and 

27  . 
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5.7  to  5 . 12  are  shown  in  pairs  showing  tho 
detection,  recognition,  and  identification  ranges  for  each 
degree  of  polarization  0,  2%  and  8%  of  the  target.  In  each 

case  the  background  polarization  is  vertical  and  dependent 
on  range,  computed  using  the  SEARAD  Radiance  Model.  Each 
pair  consists  of  a  wide  range  plot  on  which  is  shown  the 
intersection  point  of  the  ATD  curves  for  no  polarization 
filter  and  for  horizontal  filter  with  the  same  target 
polarization  values.  As  expected  in  each  case  the  critical 
range  is  increased.  The  second  plot  of  each  pair  shows  an 
expanded  plot  of  ATD  versus  MRTD  to  display  the  limiting 
ranges  more  precisely.  On  this  scale  the  no-filter  case  for 
the  detection  task  cannot  be  included. 

In  Table  5.2  the  summary  of  the  results,  obtained 
using  the  model  is  given  for  the  specific  scenario 
conditions.  The  scenario  consists  of  a  ship  target  with  a 
critical  dimension  of  24.71m.  The  ship  target  model  defined 
in  the  MuSES  program  is  similar  to  the  Cyclone  class  ships 
(170  feet)  which  serve  as  patrol  coastal  boats-PC.  In  the 
scenario  the  sea  background  has  average  13%  vertical 
polarization  while  the  target  has  0%,  2%,  and  8%  horizontal 

polarization.  The  sensor  used  in  the  model  is  a  generic 
second  generation  FLIR  sensor  mounted  on  an  aircraft  that 


80 


elevation.  All  the 


approaches  the  target  at  300m 
calculations  and  predictions  in  the  model  are  made  within 

the  8-12jim  band. 


Range  (km) 

No  Polarization 

Polarization  Filter  Included 

Filter 

Target  Degree  of  Polarization 
(Horizontal) 

Specific  Task: 

0% 

2% 

8% 

Detection  (N=0.75) 

23 

30.70 

31.20 

32.15 

Recognition  (N=3) 

19.15 

21.60 

21.80 

22.25 

Identification  (N=6) 

13.8 

15.60 

15.75 

16.10 

Table  5.2  Summary  of  Maximum  Range  Prediction  Results 
Obtained  From  Application  of  The  Model  for  The  Described 
Scenario. 


The  model  was  applied  for  only  one  set  of 
meteorological  conditions  and  only  one  target .  Thus  there 
is  basically  only  one  data  point  for  each  criterion,  i.e., 
detection,  recognition,  and  identification. 


81 


Rang*  {km  ) 


Figure  5.9  ATD  &  MRTD  vs  Range  for  Recognition  Task 
(Cycles  Across  Critical  Dimension  (N)=3)  With  A 
Target  Polarization  0,+2%,+8%  and  Background. 
Polarization  Between  -12.7%  -14.6%. 


ATD  &  MRTD  (N  =  3) 


ATD  (C) 


ATD&MRTD  (N  = 


VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  main  objectives  of  this  project  were  to  determine 
the  effect  of  degree  of  polarization  on  maximum  range 
performance  analysis  using  an  advanced  signature  prediction 
model,  and  also  to  verify  previous  project  results  and 
models  which  have  studied  target-background  contrast 
improvement  using  horizontal  polarizers  on  sensors. 

After  obtaining  the  ATD  plot  for  the  given  conditions 
and  the  MRTD  plot  of  the  second-generation  sensor,  it 
becomes  easy  to  determine  the  effect  of  target  degree  of 
polarization  on  range  improvements.  As  defined  before  the 
intersection  of  the  ATD  with  the  MRTD  shows  the  maximum 
range  at  which  the  task  can  be  performed. 

The  main  conclusions  of  this  work  and  recommendations 
are  listed  below: 

•  Figure  5.7,  5.9  and  5 . 11  show  that  there  is 

certainly  an  improvement  in  maximum  range  of  the 
sensor  for  all  tasks  when  a  horizontal  filter  is 
included  provided  that  the  target  does  not  have  a 
negative  degree  of  polarization.  The  summary  of  the 
improvement  amounts  as  percentage  of  range  with 
polarizer  is  shown  in  Figure  6.1.  The  figure 
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Improvement  In  Max. Range 


indicates  that  when  the  target  has  0-6  degree  of 
polarization,  for  a  detection  task  the  improvement 
with  respect  to  an  unpolarized  target  is  770m 
(33.4%)  for  the  recognition  task  245m  (12.79%),  and 
for  the  identification  task  180m  (13.04%). 


Maximum  Range  Improvements 


45.00% 

40.00% 

35.00% 

30.00% 

25.00% 

20.00% 

15.00% 

10.00% 

5.00% 

0.00% 


39.78% 


2%  8% 

Target  Degree  of  Polarization 


□  Recognition  Task, 
Left  Bar 


□  Identification 
Task,  Middle  Bar 


□  Detection  T ask, 
Right  Bar 


Figure  6.1  Maximum  Range  Improvement  Values  For 
Specific  Tasks. 


•  The  model  used  in  this  study  involved  only  one 
target-sensor  scenario  under  one  meteorological 
condition,  so  no  variations  with  conditions  can  be 
concluded.  To  evaluate  the  effectiveness  of 

polarization  filtering  extensive  analysis  is  needed. 
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•  The  analysis  is  based  on  the  assumption  that  the 

target  has  a  non-negative  degree  of  polarization. 
More  study  and  a  polarized  target  model  are  needed 
to  understand  the  target  polarization 

characteristics  and  the  technigues  for  improving 
polarization  values  of  man  made  targets . 

•  In  the  analysis,  SEARAD  Radiance  Model  transmittance 

values  are  used  for  determination  of  apparent 
radiance  values  at  the  specified  wavelengths.  It  is 
known  that  SEARAD  calculates  the  average 

transmittance  value  over  the  selected  bandwidth,  and 
this  averaging  might  cause  a  suppression  of  strong 
emission  or  absorption  sub-bands  within  the 
specified  wavelength  range.  To  minimize  this  source 
of  possible  error,  the  8-12|im  region  over  which  the 
whole  analysis  is  done  was  divided  into  four  sub¬ 
bands  and  the  corresponding  spectral  transmittance 
values  computed.  In  the  analysis  the  spectral 
transmittance  values  were  multiplied  by  the  spectral 
radiance  values  of  the  source  to  obtain  apparent 
radiance  values .  The  results  showed  that  the 
difference  between  the  apparent  radiance  values 
within  8-12(jm  found  by  using  spectral  values  and 
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average  transmittance  values  is  not  very  significant 
for  the  scenario  in  the  model .  The  average 
difference  was  found  to  be  less  than  1  Watts/m2sr, 
spectral  values  0.442%  greater  than  average  values. 
The  difference  becomes  less  significant  with 
increasing  range.  Although  this  approach  appears  to 
have  small  effect  on  the  overall  calculations,  it 
may  be  an  important  factor  in  regions  other  than  8- 

12pm. 

•  In  the  whole  analysis  some  points  require  special 
care  such  as  path  radiance  characteristics.  It 
should  be  remembered  that  path  radiance  is  a 
component  of  both  the  sea  background  radiance  and 
target  radiance,  and  has  a  0  degree  of  polarization. 
Since  at  great  ranges  the  radiance  incident  on  the 
sensor  will  be  mostly  from  the  path  radiance,  the 
difference  between  target  apparent  temperature  and 
the  background  apparent  temperature  should  approach 
zero.  In  Figure  5.6  the  apparent  temperatures  of 
target  and  background  approach  each  other 
asymptotically.  This  result  can  be  seen  only  if  the 
effect  of  the  polarized  filter  on  each  radiance 
element  is  studied  separately.  In  this  analysis  the 
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above  procedure  is  applied  and  the  expected  results 
are  obtained.  In  previous  attempts  to  model  this 
phenomenon,  the  calculated  ATD  did  not  behave 
logically-i . e .  ,  tend  to  zero  at  long  range.  This  is 
believed  to  be  due  to  the  different  approaches  used 
to  determine  Trarget  (R)  an<^  Tsackground  (R)  •  use  of  MuSES 
and  better  modeling  here  give  more  acceptable 
behavior  for  ATD  and  as  a  result  for  maximum  range 
calculations .  This  point  is  one  the  most  important 
conclusions  of  this  project.  A  better  modeling  of 
ATD  calculations  comparing  the  transmittance 
calculation  method  with  Beer's  Law  modeling,  and 
integration  of  degree  of  polarization  into  the 
calculations  may  also  contribute  to  further  studies. 

•  Another  point  that  should  be  noted  is  the 
"fluctuation"  of  the  apparent  temperature  difference 
in  the  ranges  up  to  5km  (Figure  5.5) .  This  phenomena 
is  seen  both  in  unpolarized  and  in  polarized 
signatures  and  has  been  seen  and  interpreted  in  some 
previous  studies.  Although  no  explanation  has  been 
advanced  based  on  physical  events,  a  possible  reason 
for  this  "fluctuation"  is  the  variation  in  the 
spectral  radiance  of  the  sky  at  different  angles  of 
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Predicted  Radiance  Value  ( Watts  /  m2sr ) 


elevation  and  at  different  ranges.  The  parameters  of 
ATD  are  path  radiance  and  the  signatures  from  target 
and  background.  Target  signature,  and  the  sea 
surface  emission  component  of  background  signature 
follow  an  exponentially  decaying  form  with 
increasing  range,  while  path  radiance  follows  an 
exponentially  increasing  course  (in  the  form  of  (1- 
e’x)  )  .  As  seen  in  Figure  6.2  sky  reflection  shows  an 
unsteady  character  in  short  ranges. 

Variation  of  Background  Radiation  Components  with  Range 


—  S ky  Reflectant  Rad. 
—♦—Path  Radiance 
Sea  Emission 
- Total  Background 


0  y-r  - — - . - — - 1 - -  "* 

0.3  1.3  2.3  3.3  4.3  5.3 

Range  (km) 

Figure  6.2  Variation  of  Background  Radiance  Components  and 
Total  Background  Radiance  with  Slant  Range. 
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•  The  issue  of  conversion  of  source  temperature  to  and 
from  in-band  radiance  is  also  one  of  the  points  that 
requires  further  study.  In  this  project  conversion 
is  done  by  direct  use  of  Planck's  formula.  One 
surprising  consequence  was  seen  when  the  apparent 
temperature  differences  were  compared  with 
corresponding  radiance  differences  (Delta  W)  .  Having 
nearly  a  hundred  data  pairs,  a  linear  Microsoft 
Excel  regression  model  was  applied  and  an  almost 
linear  relation  found  between  ATD  and  Delta  W.  The 
result  simply  denotes  that  the  ATD  can  be 
approximated  for  a  given  range  by  multiplying  the 
apparent  radiance  difference  between  the  target  and 
the  background  by  a  constant  (in  this  particular 
scenario)  0.3516237  i.e., 

ATD ( R ) =  0.3516237  *  Delta  W  (6.1) 

The  approach  is  more  promising  at  the  ranges  greater 
than  1  km.  The  results  are  presented  in  Appendix  K. 

It  should  be  expressed  clearly  that  this  study  is  a 
successor  of  extensive  studies  previously  performed  in 
performance  measurement  modeling.  Although  this  study 
developed  an  improved  maximum  range  calculation  model 
compared  to  the  models  in  previous  studies  [Ref.  18,  23, 
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27]  it  is  not  yet  complete  due  to  the  lack  of  a  polarized 
target  model.  As  a  suggestion  for  future  work,  it  would  be 
appropriate  to  modify  the  whole  analysis  by  including  a 
polarized  target  model  and  a  complete  polarized  MRTD  model. 
The  results  of  this  work,  with  the  results  of  others  in  the 
same  area,  may  be  useful  to  demonstrate  the  effect  of 
target  polarization  on  performance  of  a  sensor  in  theory. 
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APPENDIX  A 


SEARAD  RADIATION  CODE 


Tljis  is  an  example  of  the  SEARAD  Radiance  Code  input  file. 

Ti?e  atmospheric  model  and  weather  data  are  obtained  from  the  weather  file. 


11000000000  288.15  0.00 

0000  0.000  5.000  6.000  .000  .000 

.000  90.680  0.000  .000  0.00  0  90.000  T 

1250  10  2 


This  is  an  example  of  the  SEARAD  output  file  for  above  input  file. 


*****  SEARAD,  A  MODIFICATION  OF  LOWTRAN7  ***** 


DATE:  11/27/2000  TIME:  11:24:34.38 

THERMAL  RADIANCE  MODE 

MULTIPLE  SCATTERING  USED 

SLANT  PATH  TO  SPACE 

HI  =  .300  KM 

HMIN  =  .000  KM 

ANGLE  =  90.680  DEG 

FREQUENCY  RANGE 

IV1  =  830  CM-1  (  12.05  MICROMETERS) 

XV2  =  1250  CM-1  (  8.00  MICROMETERS) 

IDV  =  10  CM-1 

IFWHM  =  2  CM-1 

IFILTER  =  0 

SUMMARY  OF  THE  GEOMETRY  CALCULATION 

HI  =  .300  KM 

H2  =  .000  KM 

ANGLE  =  90.680  DEG 

RANGE  =  30.360  KM 

BETA  =  .273  DEG 

PHI  =  89.547  DEG 

HMIN  =  -000  KM 

BENDING  =  .046  DEG 

LEN  =  0 


2  0 
00.300 
833 

0 
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SEA  AT  288.15  K  REPLACES  BLACK  BODY  BOUNDARY 


UPWIND  =  90.000  DEG  EAST  OF  LINE  OF  SIGHT 


ZERO  RANGE  UNPOLARIZED  VALUES 

SEA  EMISSION 
SKY  REFLECTION 
SUN  GLINT 

TOTAL  RADIANCE 
BLACK  BODY  TEMP. 


FULL  RANGE  UNPOLARIZED  VALUES 

SEA  EMISSION 
SKY  REFLECTION 
SUN  GLINT 
PATH  TO  FOOTPRINT 

TOTAL  RADIANCE 
BLACK  BODY  TEMP. 


ZERO  RANGE  POLARIZED  VALUES 

SEA  EMISSION 
SKY  REFLECTION 
SUN  GLINT 

'  TOTAL  RADIANCE 
BLACK  BODY  TEMP.  (C) 


18.01540  W  M-2  SR-1 
11.41543  W  M-2  SR-1 
.00000  W  M-2  SR-1 


(AV.  EMISS.  .5637) 


29.43082  W  M-2  SR-1 
10.6  C 


.01017  W  M-2  SR-1 
.00622  W  M-2  SR-1 
.00000  W  M-2  SR-1 
33.55154  W  M-2  SR-1 


(AV.  TRANS.  .0006) 


33.56792  W  M-2  SR-1 
18.0  C 


HORIZONTAL 
(W  M-2  SR-1) 


VERTICAL 
(W  M-2  SR-1) 


(H-V)  /  (H+V) 

(%) 


7.93007  10.08533  -12.0 

6.52375  4.89168  14.3 

.00000  .00000  .0 

14.45381  14.97701  -1.8 

-24.0  -22.5 
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FULL  RANGE  POLARIZED  VALUES 


SEA  EMISSION 

.00446 

.00571 

-12.3 

SKY  REFLECTION 

.00356 

.00265 

14.6 

SUN  GLINT 

.00000 

.00000 

.0 

PATH  TO  FOOTPRINT 

16.77577 

16.77577 

.0 

TOTAL  RADIANCE 

BLACK  BODY  TEMP.  (C) 

16.78379 

-17.5 

16.78413 

-17.5 

.0 
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APPENDIX  B 


IR  SIGNATURE  SIMULATION  SOFTWARES 

This  Appendix  contains  information  on  IR  signature  simulation, 
materials  and  background. modeling,  and  performance  prediction 
programs.  Some  of  the  programs  are  commercial  and  have  been 
widely  used.  Descriptions  of  the  programs  are  taken  from  Ref. 

13 .  Contact  points  and  more  information  about  the  programs  can 
be  found  in  the  same  reference. 

AEMAT  :Avlab  Electro-optical  Model  for  Aerial  Targeting 
This  model  also  known  as  EMAT  (Electro-optical  Model  for  Aerial 
■Pargeting) .  It  has  some  surface  target  analysis  capability  and 
includes  targets,  atmospheres,  backgrounds,  and  sensor  models  in 
all  wavebands.  Incorporates  SPIRITS  by  Aerodyne  Research.  Owner 
of  the  model  is  Horizons  Technology,  Inc. 

CHARM:  High  Altitude  Radiation  Model 

This  model  is  a  plume  flow  and  radiation  model  for  non-air 
breathing  vehicles  at  altitudes  above  40km,  the  model  SIRRM  is 
the  plume  radiation  model  applicable  for  both  air  breathing  and 
non-air  breathing  vehicles.  Owner  of  the  model  is  Photon 
Research  Associates. 

CREEP  :  Coating  Reflectance  Engineering  Evaluation  Code 

First  principle  method  for  computation  of  BRDF  of  paint 

coatings .  Performs  both  surface  and  volume  computations .  For 
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subsurface  material,  the  code  uses  computed  Mie  scattering  phase 
functions  with  multiple  scattering  computational  techniques  to 
determine  the  equivalent  BRDF .  Owner  of  the  model  is  ERIM. 

EMAT:  Electro-optical  Model  for  Aerial  Targeting 

GTSIG  (also  known  as  TCMII) :  Georgia  Tech  Thermal  Contrast 
Model.  Thermal  model  development  with  GTSIG  starts  with  a 
geometrical  model  and  allows  user  to  incorporate  thermal 
conditions  such  as  thermal  loading,  hardbody  heat  sources,  and 
so  forth.  It  models  both  targets  and  backgrounds.  GTSIG  is  also 
known  as  TCMII,  which  is  the  same  program  that  is  owned  by 
Georgia  Tech.  Research  Institute. 

IASPM  :  Infrared  Atmospherics  and  Signatures  prediction  Model 
IRMA:  Infrared  Modeling  and  Analysis 

It  models  passive  IR,  C02  laser  radar,  and  passive  millimeter 
wave  systems.  Output  is  high  resolution  imagery.  Models  ground 
mobile  and  high  value  fixed  targets  and  backgrounds.  Owner  of 
the  model  is  Nichols  Research  Corporation. 

PCNirATAM:  NATO  infrared  Air  Target  Model 

This  is  a  complete  package  to  compute  the  infrared  signature  of 
an  air  vehicle  and  the  associated  atmosphere.  The  package 
includes  executable  code  for  data  input,  screen  and  hardcopy 
plotting,  image  processing,  and  target  signature  calculations. 
Owner  of  the  model  is  Ontar  Corp. 
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PRXSM:  Physically  Reasonable  Infrared  Signature  Model 
This  model  predicts  the  surface  temperature  and  radiance  of 
vehicles,  soil,  and  vegetative  background  from  measured  or 
modeled  environment  conditions.  Radiometric  calculations  do  not 
include  atmospheric  or  sensor  effects.  The  results  are  sometimes 
referred  to  as  zero  distance  source  values .  Owner  of  the  model 
is  Keweenaw  Research  Center,  Michigan  Tech.  University. 

SIRIM  :  Simulated  Infrared  Image  Model 

This  model  is  used  to  generate  simulated  thermal  imagery  of 
targets  embedded  in  synthetic  or  collected  backgrounds.  SIRIM 
combines  CSG  (combinatorial  solid  geometry)  model,  an  automatic 
mesh  generator,  a  thermal  diffusion  model,  and  the  radiance 
calculation  into  one  integrated  thermal  signature  prediction 

code.  Owner  of  the  model  is  ERIM. 

SIRRM :  Standardized  Infrared  Radiation  Model 
SOSIS:  Ship  and  Ocean  Surface  Image  Simulation 
SPF:  Standard  Plume  Flow  Model 

SPIRITS:  Spectral  and  In-  Band  Radiometric  Imaging  of  Targets 
and  Scenes  Model.  SPIRITS  is  a  first-principles  target  signature 
model  developed  originally  for  aircraft,  but  now  applicable  to 
tactical  missiles  and  low  altitude  strategic  missiles.  It 
combines  a  variety  of  models  to  produce  a  flexible  and  high 
performing  target /signature  modeling  capability.  It  incorporates 
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LOWTRAN  or  MODTRAN  for  atmospherics.  Owner  of  the  model  is 


Aerodyne  Research,  Inc. 

SSGM:  Strategic  Scene  Generation  Model 
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APPENDIX  C 


TARGET  DIMENSIONS/  ORIENTATION  and  GEOMETRY  of  COMPUTATIONS 


<j)  '=  Elevation  Angle  (~0  degree) 

6  =  Azimuth  Angle  (~70  degrees) 
1  =  Length  (64.832  meters) 
h  =  Height  (7.934  meters) 
w  =  Width  (7.143  meters) 


J 

THIS  PAGE  INTENTIONALLY  LEFT  BLANK 
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APPENDIX  D 


SAMPLE  WEATHER  FILE  FOR  COMPUTATIONS 


The  Natural  environment (weather )  option  in  MuSES  uses  models  or 
measurements  of  the  weather  conditions  and  the  solar  position  to 
calculate  the  cooling  or  heating  that  occurs  from  the  influence  of 
weather . Currently  there  are  four  different  weather  files  present  in 
the  program.  The  parameters  are  as  follows: 

Global  Positioning: 

Latitude  (deg.): 47. 2  Longitude  (deg.):  88.5 

Time  Zone  (hrs.):4  Elevation  (m) :  332 

Date  of  The  Measurements:  07/19/1984. 

Available  Solution  Times:  From  07/19/1984  06:20  AM  to  07/20/1984 
08:59  AM 


Selected  Solution  Time:  07/19/1984  18:00  PM.  Late  afternoon  before 
dusk.  This  weather  file  is  selected  and  made  compatible  with  a 
standard  MODTRAN  atmosphere  (Midlatitude  Summer)  by  modifying  some 
of  the  default  surface  parameters . 


Time 

Air  Temp 
(C°) 

in 

Wind 

Speed 

(m/sec) 

Humidity  (%) 

■  ■ 
■ 

■nsm 

Wind  Direct. 
(Deg.) 

6:10  AM 

18.8 

0 

1.299 

78.833 

5 

155.808 

6:15  AM 

18.9 

2 

0.828 

78.25 

6 

150.912 

6:20  AM 

18.7 

4 

0.903 

77.667 

6 

151.92 

6:25  AM 

18.6 

6 

1.532 

77.083 

6 

180.216 

6:30  AM 

18.7 

10 

1.055 

76.5 

7 

169.704 

6:35  AM 

18.4 

20 

1.63 

75.917 

7 

162.72 

6:40  AM 

18.4 

24 

1.016 

75.333 

7 

171.936 

6:45  AM 

18.4 

34 

1.39 

74.75 

8 

159.408 

6:50  AM 

18.5 

50 

1.131 

74.167 

8 

149.184 

6:55  AM 

19.7 

116 

2.397 

73.583 

8 

137.088 

7:00  AM 

18.8 

94 

1.633 

73 

7 

163.512 

7:05  AM 

18.5 

34 

0.467 

72.583 

7 

182.232 

7:10  AM 

18.4 

26 

1.758 

72.167 

7 

147.456 

7:15  AM 

18 

14 

2.811 

71.75 

6 

161.28 

7:20  AM 

17.8 

16 

1.647 

71.333 

6 

142.056 

7:25  AM 

17.8 

14 

0.659 

70.917 

6 

170.928 

7:30  AM 

17.9 

12 

3.116 

70.5 

5 

157.248 

7:35  AM 

17.8 

10 

1.309 

70.083 

5 

214.848 

7:40  AM 

18 

10 

1.679 

69.667 

4 

148.464 

7:45  AM 

18.2 

28 

2.012 

69.25 

4 

166.32 

7:50  AM 

18.5 

70 

2.727 

68.833 

4 

152.496 

7:55  AM 

18.8 

160 

2.131 

68.417 

4 

168.48 

8:00  AM 

19.5 

282 

2.072 

68 

4 

156.528 

8:05  AM 

20 

306 

2.419 

67.917  | 

3 

153.792 
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8:10  AM 


8:15  AM 


8:20  AM 


8:25  AM 


8:30  AM 


8:35  AM 


8:40  AM 


8:45  AM 


8:50  AM 


8:55  AM 


9:00  AM 


9:05  AM 


9:10  AM 


9:15  AM 


9:20  AM 


9:25  AM 


9:30  AM 


9:35  AM 


9:40  AM 


9:45  AM 


9:50  AM 


9:55  AM 


10:00  AM 


1 0:05  AM 


10:10  AM 


10:15  AM 


10:20  AM 


10:25  AM 


10:30  AM 


10:35  AM 


10:40  AM 


10:45  AM 


10:50  AM 


10:55  AM 


11:00  AM 


11:05  AM 


11:10  AM 


11:15  AM 


11:20  AM 


11:25  AM 


11:30  AM 


11:35  AM 


11:40  AM 


11:45  AM 


11:50  AM 


11:55  AM 


12:00  PM 


12:05  PM 


12:10  PM 


1.402 

67.833 

3 

144.432 

1.209 

67.75 

3 

176.544 

2.516 

67.667 

3 

172.368 

1.515  i 

67.583 

3 

182.232 

1.483 

67.5  1 

3 

147.528 

2.031 

67.417 

3 

151.776 

1.293 

67.333 

2 

237.168 

1.631 

67.25 

2 

173.448 

1.807 

67.167 

2 

159.48 

1.446 

67.083 

2 

140.472 

1.967 

67 

2  I 

148.464 

1.504 

66.917 

2 

164.16 

1.03 

66.833 

2 

160.416 

2.862 

66.75 

2 

189.936 

1.435 

66.667 

2 

159.048 

1.752 

66.583 

1 

132.552 

2.231 

66.5 

1 

125.352 

0.707 

66.417 

1 

172.008 

2.117 

66.333 

1 

160.344 

1.52 

66.25 

1 

163.368 

2.134 

66.167 

1 

169.992 

1.522 

66.083 

1 

169.272 

3.271 

66 

1 

167.328 

1.01  i 

66.083 

1 

159.984 

1.016 

66.167 

1  1 

124.344 

1.779  ! 

66.25 

1 

170.208 

2.04 

66.333 

1 

154.944 

1.791 

66.417 

0 

168.264 

3.125 

66.5 

0 

169.92 

1.561 

66.583 

0 

175.824  [ 

2.837 

66.667 

0 

174.888 

0.753 

66.75 

0 

186.84 

0.916 

66.833 

0 

192.816 

1.224 

66.917 

0 

163.44 

1.726 

67 

0 

1 58.688 

1.956 

66.917 

0 

171.936 

1.387 

66.833 

0 

1  198.792 

0.884 

66.75 

0 

192.672 

1.577 

66.667 

0 

159.552 

1.431 

66.583 

0 

211.824 

2.232 

66.5 

0 

170.352 

2.599 

66.417 

0 

188.856 

0.714 

66.333 

0 

203.76 

2.549 

66.25 

0 

186.912 

0.992 

66.167 

0 

155.232 

1.933 

66.083 

0 

207.144 

1.344 

66 

0 

201.384 

2.328 

65.75 

0 

176.976 

1.727 

65.5 

0 

189.216 

104 


12:15  PM 


12:20  PM 


12:25  PM 


12:30  PM 


12:35  PM 


12:40  PM 


12:45  PM 


12:50  PM 


12:55  PM 


1:00  PM 


1:05  PM 


1:10  PM 


1:15  PM 


1:20  PM 


1:25  PM 


1 :30  PM 


1:35  PM 


1:40  PM 


1:45  PM 


1:50  PM 


1 :55  PM 


2:00  PM 


2:05  PM 


2:10  PM 


2:15  PM 


2:20  PM 


2:25  PM 


2:30  PM 


2:35  PM 


2:40  PM 


2:45  PM 


2:50  PM 


2:55  PM 


3:00  PM 


3:05  PM 


3:10  PM 


3:15  PM 


3:20  PM 


3:25  PM 


3:30  PM 


3:35  PM 


3:40  PM 


3:45  PM 


12:00  AM 


3:55  PM 


4:00  PM 


4:05  PM 


4:10  PM 


21.2 


21.3 


22.4 


21.5 


21.6 


21.9 


22.1 


21.8 


22.2 


22.2 


21.9 


21.8 


21.7 


21.8 


22.2 


21.9 


22.3 


22.4 


23.1 


22.8 


22.5 


22.7 


23.1 


22.8 


23.3 


22.9 


22.9 


23.9 


23.8 


23.6 


23.5 


23.4 


23 


23.3 


23.6 


23.8 


22.9 


23.3 


23.3 


23.7 


23.2 


23.3 


23.8 


23.7 


23.9 


23.7 


24 


23.8 


2.532 


1.106 


1.661 


2.452 


1.366 


1.851 


1.342 


2.102 


3.257 


1.933 


1.944 


2.694 


1.807 


1.839 


1.284 


2.331 


1.232 


1.541 


1.55 


1.375 


2.155 


1.243 


1.073 


1.734 


3.066 


2.306 


1.891 


0.649 


1.251 


1.179 


1.489 


1.439 


2.124 


1.72 


1.198 


1.814 


3.068 


0.935 


2.096 


0.808 


1.257 


1.646 


1.209 


0.746 


3.028 


2.494 


1.487 


1.95 


65.25 

0 

176.472 

65 

0 

152.496 

64.75 

0 

210.888 

64.5 

0 

190.584 

64.25 

0 

188.28 

64 

0 

219.6 

63.75  I 

0 

192.96 

63.5 

0 

153.576 

63.25 

0 

157.536 

63 

0 

161.136 

63.167 

0 

148.248 

63.333 

0 

184.176 

63.5 

0 

167.544 

63.667 

0 

172.584 

63.833 

0  ' 

163.656 

64 

0 

189.936 

64.167 

0 

168.408 

64.333 

0 

163.728 

64.5 

0 

198.936 

64.667  ] 

0  j 

170.928 

64.833 

0 

167.544 

65 

0 

171.072 

65.833 

0 

186.192 

66.667 

0 

153.792 

67.5 

0 

198.288 

68.333 

0 

147.168 

69.167 

0 

143.136 

70 

0 

139.464 

70.833 

0 

151.92 

71.667 

0 

201 .744 

72.5 

0 

166.104 

73.333 

0 

189.072 

74.167 

0 

191.52 

75 

0 

196.344 

73.833 

0 

157.032 

72.667 

0 

227.088 

71.5 

0 

163.8 

70.333 

0 

188.784 

69.167 

0 

144.936 

68 

0 

166.032 

66.833 

0 

175.896 

65.667 

0 

162.648 

64.5 

0 

199.008 

63.333 

0 

190.584 

62.167 

0 

190.368 

61 

0 

164.88 

59 

0 

186.984 

57 

0 

176.328 

105 


4:15  PM 

24.1 

758 

2.375 

55 

0 

176.472 

4:20  PM 

24.2 

750 

0.931 

53 

0 

156.672 

4:25  PM 

24.2 

744 

1.568  j 

51 

0 

187.56 

4:30  PM 

24.3 

738 

1.75 

49 

0 

218.592 

4:35  PM 

24.7 

714 

1.078 

47 

0 

154.944 

4:40  PM 

24.4 

708 

1.027 

45 

0 

158.256 

4:45  PM 

24.4 

700 

1.411 

43 

0 

151.488 

4:50  PM 

24.6 

692 

0.788 

41 

0  i 

192.24 

4:55  PM 

24.8 

684 

1.224 

39 

0  1 

220.176 

5:00  PM 

25 

674 

1.247 

37 

0 

156.744 

5:05  PM 

24.4 

662 

1.431 

37.083 

0 

193.032 

5:10  PM 

24 

650 

2.481 

37.167 

0 

183.456 

5:15  PM 

24 

640 

1.322 

37.25 

0 

163.08  | 

5:20  PM 

24.6 

626 

2.553 

37.333 

0 

164.232 

5:25  PM 

24.3 

618 

1.084 

37.417 

0 

177.48 

5:30  PM 

24.4 

604 

0.714 

37.5 

0 

173.592 

5:35  PM 

24.2 

594 

1.37 

37.583 

0 

215.928 

5:40  PM 

24.7 

582  1 

0.757 

37.667 

o  1 

139.176 

5:45  PM 

24.7 

568 

0.693 

37.75 

0 

184.608 

5:50  PM 

24.8 

552 

1.349 

37.833 

0 

180.36 

5:55  PM 

24.6 

542 

1.225 

37.917 

0 

194.112 

6:00  PM 

24.4 

530 

0.968 

38 

0 

196.344 

6:05  PM 

24.2 

516 

1.418 

38.417 

0 

194.256 

6:10  PM 

25.1 

502 

0.554 

38.833 

0 

179.784 

6:15  PM 

24.3 

488 

0.86 

39.25 

0 

176.976 

6:20  PM 

24.3" 

474 

0.981  1 

39.667 

0 

168.84 

6:25  PM 

24.7 

460 

0.941 

40.083 

0 

204.768 

6:30  PM 

24.4 

444 

0.824 

40.5 

0 

201.168 

6:35  PM 

24.1 

434 

1.095 

40.917 

0 

223.2 

6:40  PM 

23.9 

418 

0.87 

41.333 

0 

185.76 

6:45  PM 

23.8 

400 

0.973 

41.75 

0 

201.528 

6:50  PM 

23.7 

386 

1.339 

42.167 

0 

I  178.632 

6:55  PM 

23.9 

374 

1.111 

42.583 

0 

190.152 

7:00  PM 

23.9 

358 

1.211 

43 

0 

209.448 

7:05  PM 

23.9 

344 

1.674 

43 

0 

164.952 

7:10  PM 

23.9 

328 

0.624 

43 

0 

178.2 

7:15  PM 

23.7 

316 

1.223 

43 

0 

212.472 

7:20  PM 

23.3 

300 

0.753 

43 

0 

199.44 

7:25  PM 

23.6 

286 

0.791 

43 

0 

181.08 

7:30  PM 

23.4 

272 

0.482 

43 

0 

183.96 

7:35  PM 

23.3 

260 

0.622 

43 

0 

140.904 

7:40  PM 

23.7 

240 

0.626 

43 

0 

181.224 

7:45  PM 

23.4 

230 

0.505 

43 

0 

1  167.544 

7:50  PM 

23.4 

216 

0.303 

43 

0 

204.912 

7:55  PM 

23 

202 

0.783 

43 

0 

186.768 

8:00  PM 

22.4 

184 

0.2 

43 

0 

198.648 

8:05  PM 

22 

176 

0.337 

44.167 

0 

216.36 

8:10  PM 

21.9 

168 

0.2 

45.333 

0 

141.912 

8:15  PM 

21.8 

154 

0.202 

46.5 

0 

186.768 
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12:24  AM 


12:29  AM 


12:34  AM 


12:39  AM 


12:44  AM 


12:49  AM 


12:54  AM 


12:59  AM 


1:04  AM 


1 :09  AM 


1:14  AM 


1:19  AM 


1:24  AM 


1:29  AM 


1:34  AM 


1:39  AM 


1:44  AM 


1:49  AM 


1:54  AM 


1:59  AM 


2:04  AM 


2:09  AM 


2:14  AM 


2:19  AM 


2:24  AM 


2:29  AM 


2:34  AM 


2:39  AM 


2:44  AM 


2:49  AM 


2:54  AM 


2:59  AM 


3:04  AM 


3:09  AM 


3:14  AM 


3:19  AM 


3:24  AM 


3:29  AM 


3:34  AM 


3:39  AM 


3:44  AM 


3:49  AM 


3:54  AM 


3:59  AM 


4:04  AM 


4:09  AM 


4:14  AM 


4:19  AM 


4:24  AM 


13.7 


14.3 


14.4 


14.6 


14.6 


14.7 


14.6 


14.8 


14.5 


14.5 


14.2 


14.6 


14.5 


14.4 


13.9 


14.6 


14.6 


14.6 


14.8 


14.8 


14.9 


14.8 


14.8 


14.6 


14.2 


13.8 


13.3 


13.3 


13.2 


13.3 


13.6 


13.8 


13.9 


13.3 


13.3 


13.3 


13 


13.3 


13.3 


13.2 


13.3 


13.9 


13.9 


14.1 


13.8 


13.2 


12.6 


12.6 


12.4 


0.2 


0.2 


0.704 


0.2 


0.336 


0.227 


0.302 


0.218 


0.252 


0.205 


0.2 


0.421 


0.29 


0.369 


0.199 


0.266 


0.291 


0.531 


0.887 


0.415 


0.498 


0.304 


0.514 


0.855 


0.344 


0.747 


0.213 


0.2 


0.529 


0.938 


0.423 


0.572 


0.599 


0.28 


0.352 


0.64 


1.227 


0.761 


0.627 


0.689 


1.235 


0.551 


1.011 


0.877 


1.102 


0.341 


0.397 


0.595 


0.432 


67.578 


67.837 


68.096 


68.356 


68.615 


68.874 


69.133 


69.393 


69.652 


69.911 


70.17 


70.43 


70.689 


70.948 


71.207 


71.467 


71.726 


71.985 


72.244 


72.504 


72.763 


73.022 


73.281 


73.541 


73.8 


74.059 


74.319 


74.578 


74.837 


75.096 


75.356 


75.615 


75.874 


76.133 


76.393 


76.652 


76.911 


77.17 


77.43 


77.689 


77.948 


78.207 


78.467 


78.726 


78.985 


79.244 


79.504 


79.763 


80.022 


166.968 


172.584 


175.752 


180.648 


176.976 


184.248 


180.144 


183.6 


177.192 


170.208 


187.488 


191.664 


225.792 


214.848 


187.704 


188.208 


182.808 


182.88 


190.08 


206.424 


187.2 


206.784 


195.336 


193.248 


193.896 


159.768 


167.76 


178.128 


168.768 


173.16 


165.024 


175.752 


188.64 


201.168 


179.136 


179.784 


162.072 


175.248 


187.272 


169.632 


191.304 


193.392 


194.688 


207.648 


191.88 


182.304 


192.168 


194.04 


180.144 
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4:29  AM 

12.3 

4:34  AM 

12.2 

4:39  AM 

12.1 

4:44  AM 

12.1 

4:49  AM 

12.4 

4:54  AM 

12.3 

4:59  AM 

11.7 

5:04  AM 

11.7 

5:09  AM 

12 

5:14  AM 

12.4  ! 

5:19  AM 

12.3 

5:24  AM 

11.8 

5:29  AM 

11.9 

5:34  AM 

12.1 

5:39  AM 

12 

5:44  AM 

11.9 

5:49  AM 

11.9 

5:54  AM 

12.2 

5:59  AM 

12.3 

6:04  AM 

12.3 

6:09  AM 

12.3 

6:14  AM 

12.2 

6:19  AM 

12.3 

6:24  AM 

12.3 

6:29  AM 

12.2 

6:34  AM 

12.2 

6:39  AM 

12.1 

6:44  AM 

12.8 

6:49  AM 

13 

6:54  AM 

13.2 

6:59  AM 

13.6 

7:04  AM 

13.9 

7:09  AM 

14.3 

7:14  AM 

14.3 

7:19  AM 

14.7 

7:24  AM 

15 

7:29  AM 

15.1 

7:34  AM 

15.6 

7:39  AM 

15.8 

7:44  AM 

15.9 

7:49  AM 

16.3 

7:54  AM 

16.7 

7:59  AM 

16.9 

8:04  AM 

17.2 

8:09  AM 

17.4 

8:14  AM 

17.6 

8:19  AM 

17.7 

8:24  AM 

17.9 

8:29  AM 

18.2 

0.911 


0.573 


0.283 


0.259 


0.323 


0.904 


0.98 


0.464 


0.601 


0.474 


0.64 


0.254 


0.553 


0.984 


1.041 


0.271 


0.938 


0.423 


0.812 


0.556 


0.699 


0.492 


0.698 


0.627 


0.581 


0.207 


0.57 


0.245 


0.696 


0.254 


0.875 


0.345 


0.877 


0.6 


0.81 


0.305 


1.043 


0.539 


0.57 


1.42 


0.948 


1.574 


0.471 


0.997 


0.935 


0.85 


1.625 


0.243 


0.51 


80.281 


80.541 


80.8 


81.059 


81.318 


81.578 


81.837 


82.096 


82.356 


82.615 


82.874 


83.133 


83.393 


83.652 


83.911 


84.17 


84.43 


84.689 


84.948 


85.207 


85.467 


85.726 


85.985 


86.244 


86.504 


86.763 


87.022 


87.281 


87.541 


87.8 


88.059 


88.318 


88.578 


88.837 


89.096 


89.356 


89.615 


89.874 


90.133 


90.393 


90.652 


90.911 


91.17 


91.43 


91.689 


91.948 


92.207 


92.467 


92.726 


0 

181.728 

0 

170.928 

0 

190.728 

0 

191.016 

0 

189.36 

0 

183.384 

0 

187.56 

0 

183.816 

0 

184.464 

0 

179.352 

0 

190.8 

0 

179.856  ! 

0 

180.576 

0 

169.632 

0 

168.984 

0 

174.024 

0 

166.752 

0 

163.656 

0 

164.016 

1  o 

165.384 

0 

166.392 

0 

169.272 

0 

162.792 

0 

167.112 

0 

168.048 

0 

182.16 

0 

174.816 

0 

166.824 

0 

181.728 

0 

176.616 

0 

177.264 

0 

166.608 

0 

165.6 

0 

169.272 

0 

162.504 

0 

171.216 

0 

165.096 

0 

I  157.392 

0 

166.896 

0 

172.152 

0 

180.144 

0 

169.992 

0 

172.368 

0 

170.784 

1 

180.072 

1 

181.224 

1 

173.664 

1 

170.928 

1 

169.704 
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8:34  AM 

18.4 

324 

0.924 

92.985 

1 

202.752 

8:39  AM 

18.7 

336 

1.342 

93.244 

2 

189.792 

8:44  AM 

19.2 

350 

0.519 

93.504 

2 

190.368 

8:49  AM 

19.1 

366 

1.182 

93.763 

2 

191.664 

8:54  AM 

19.2 

376 

1.207 

94.022 

2 

205.704 

8:59  AM 

19.3 

390 

1.113 

94.281 

2 

162.072 
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APPENDIX  E 


TARGET  RADIANCE  COMPUTATION  RESULTS 

#  MuSES  ascii  Results  File 

#  Exported  from  model:  c:\Program  Files\MuSES\Work\Examples\tez1 

#  TueJun  13  00:22:47  2000 

#  Units  in  the  file:  SI  Units 
Radiance  Band:  8  to  9  micron 

Time  Step  10,  Time  =  Jul  19, 1984  18:00:00,  Lat.47.20,  Long., 88.5°,  Time  Zone  4  hrs 


Element 

(Number&Location  in 
the  Model) 

Radiance 

(W/mA2-sr) 

Area(m2) 

Eff.lntensity 

(W/sr) 

■ItkBWil 

7.958443 

1.2388 

9.8589 

IMBI 

7.992607 

4.1946 

33.5258 

101  (front) 

7.967958 

3.798 

30.2623 

8.007335 

1.512 

12.1071 

8.012812 

1.729 

13.8542 

6 

99  (front)  1 

7.965985 

1.373 

10.9373 

7 

mam 

7.993531 

3.537 

28.2731 

8 

kei«i 

7.561034 

2.776 

20.9894 

9 

Bk«l 

7.579693 

4.1102 

31.1541 

inai 

7.996742 

2.459 

19.6640 

fcJdKiSSBB 

8.006999 

4.3818 

35.0851 

8.00824 

4.4529 

35.6599 

97  (front) 

8.033536 

0.3109 

2.4976 

eisosssi 

7.587668 

3.928 

29.8044 

7.598871 

5.806 

44.1190 

8.010954 

6.6406 

53.1975 

EUQS&SI 

8.161586 

0.8872 

7.2410 

KEVffltBEI 

8.163546 

1.583 

12.9229 

mxmm 

8.01196 

4.2047 

33.6879 

mm 

8.014287 

3.1853 

25.5279 

8.156629 

2.3744 

19.3671 

8.15473 

3.167 

25.8260 

23 

UiltUI 

8.031284 

1.3385" 

10.7499 

8.020765 

1.7346 

13.9128 

25 

ws^hebsm 

8.157871 

3.558 

29.0257 

26 

E!BI 

8.178525 

4.2141 

34.4651 

5H 

8.039721 

0.4487 

3.6074 

28 

m 

8.194448 

4.8426 

39.6824 

29 

8.199583 

2.2895 

18.7729 

■n vnsfisi 

8.201283 

4.8426 

39.7155 

BBBB1 

9.210635 

2.0937 

19.2843 

9.238527 

2.0938 

19.3436 

8.20176 

4.8427 

39.7187 

8.201813 

4.8426 

39.7181 

una 

9.244595 

2.0938 

19.3563 

111 


1 09  (front) 


6 

8 

m 

9 

i*mrrm 


gywil 

aiiBl 


3  3 


TOTAL 


9.241716 

2.0938 

19.3503 

8.201728 

4.8426 

39.7177 

8.201805 

4.8426 

39.7181 

9.239729 

2.0938 

19.3461 

9.214361 

2.0938 

19.2930 

8.201061 

4.8426 

39.7145 

8.201627 

4.8426 

39.7172 

9.141351 

2.0949 

19.1502 

9.116871 

2.0937 

19.0880 

8.200602 

4.8426 

39.7122 

8.200427 

4.8426 

39.7114 

9.137383 

1.9801 

18.0929 

9.203058 

2.2074 

20.3148 

8.201587 

4.8426 

39.7170 

8.201874 

4.8426 

39.7184 

9.234715 

2.0937 

19.3347 

9.193591 

2.1922 

20.1542 

8.2017 

4.8426 

39.7176 

8.200892 

4.8426 

39.7136 

9.137054 

1.9953 

18.2312 

9.087182 

2.0937 

19.0258 

8.200092" 

4.8426 

39.7098 

8.199725 

4.8426 

39.7080 

9.069571 

2.0937 

18.9890 

9.063439 

2.0937 

18.9761 

8.199563 

4.8426 

39.7072 

8.199446 

4.8426 

39.7066 

9.05103 

2.0937 

18.9501 

8.982138 

1.563 

14.0391 

8.19859 

4.84261 

39.7025 

8.206901 

4.84261 

39.7427 

8.232753 

4.873 

40.1182 

8.0648711 

1.8129 

14.6208 

8.230124 

5.9936 

49.3281 

8.651789 

6.2584 

54.1464 

8.676771 

4.7116 

40.8815 

8.680046 

4.7116 

|  40.8969 

8.677021 

4.7116 

40.8827 

8.677439 

4.7116 

40.8846 

8.663925 

5.9623 

51.6569 

8.62009 

4.1013 

35.3536 

8.629618 

3.6995 

31.9253 

8.621908 

4.0475 

34.8972 

8.67183 

3.1649 

27.4455 

8.556823 

2.9062 

24.8678 

8.830137" 

2.5541 

22.5531 

8.73445 

2.0339 

17.7650 

9.08509 

5.7346 

52.0994 

287.2313 

2401.0082 
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#  MuSES  ascii  Results  File 

#  Exported  from  model:  c:\Program  Files\MuSES\Work\Examples\tez1 

#  Tue  Jun  13  00:22:47  2000 

#  Units  in  the  file:  SI  Units 
Radiance  Band:  9  to  10  micron 

Time  Step  10,  Time  =  Jul  19, 1984  18:00:00 


Element 

(Number&Location  in 
theModel) 


cUdMiMtTCI 


Radiance 

(W/mA2-sr) 

Area(m2) 

Eff.lntensity 

(W/sr) 

8.39591 1 

1.2388 

10.4009 

8.428118 

4.1946 

35.3526 

8.404932 

3.798 

31.9219 

8.454721 

1.512 

12.7835 

8.459952 

1.729" 

14.6273 

8.403075 

1.373" 

11.5374 

8.428987 

3.537 

29.8133 

7.980424 

2.776 

22.1537 

7.998057 

4.1102" 

32.8736 

8.432007 

2.459 

20.7343 

8.454469 

4.3818 

37.0458 

8.455703 

4.4529 

37.6524 

8.466532 

0.3109" 

2.6322 

8.005523 

3.928 

31.4457 

8.016232 

5.806 

46.5422 

8.4582 

6.6406 

56.1675 

8.582376 

0.8872" 

7.6143 

8.584213 

1.583 

13.5888 

8.459149 

4.2047 

35.5682 

8.461345 

3.1853 

26.9519 

8.57773 

2.3744] 

20.3670 

8.575952 

3.167 

27.1600 

8.477506 

1.3385 

11.3471 

8.467646 

1.7346 

14.6880 

8.578896 

3.558 

30.5237 

8.598026 

4.2141 

36.2329 

8.485527 

0.4487 

3.8075 

8.613157 

4.8426 

41.7101 

8.617968 

2.2895 

19.7308" 

8.61956 

4.8426 

41.7411 

9.562263 

2.0937" 

20.0205 

9.588289 

2.0938 

20.0760 

8.620008 

4.8427 

41 .7441 

8.620056 

4.8426 

41.7435 

9.593952 

2.0938 

20.0878 

9.591261 

2.0938 

20.0822 

8.619977 

4.8426 

41.7431 

8.620049 

4.8426 

41.7434 
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MMKMI 


312  (front) 


312  (front) 


HUM 


312  (front) 


312  (front) 


312  (front) 


312  (front) 


312  (front) 


312  (front) 


liiHHil 


78 


79  312  (front) 


Ham 


83  312  (front) 


TOTAL 


9.589403 

2.0938 

20.0783 

9.565712 

2.0938 

20.0287 

8.619352 

4.8426" 

41.7401 

8.619883 

4.8426 

41.7426 

9.497645 

2.0949 

19.8966 

9.474776 

2.0937 

19.8373 

8.618921 

4.8426 

41.7380 

8.618758 

4.8426 

41.7372 

9.494009 

1.9801 

18.7991 

9.555324 

2.2074 

21.0924 

8.619844 

4.8426 

41.7425 

8.620113 

4.8426 

41.7438 

9.584821 

2.0937 

20.0677 

9.546453 

2.1922 

20.9277 

8.6199 

4.8426 

41.7427 

8.619194 

4.8426 

41.7393 

9.493543 

1.9953 

18.9425 

9.446975 

2.0937 

19.7791 

8.618444 

4.8426 

41.7357 

8.6181 

4.8426 

41.7340 

9.430511 

2.0937 

19.7447 

9.424775 

2.0937 

19.7327 

8.617949 

4.8426 

41.7333 

8.617839 

4.8426 

41.7327 

9.413185 

2.0937 

19.7084 

9.348808 

1.563 

14.6122 

8.617038 

4.8426 

41.7289 

8.624819 

4.8426 

41.7665 

8.649023 

4.873 

42.1467 

8.501824 

1.8129 

15.4130 

8.646563 

5.9936 

51.8240 

9.041137 

6.2584 

56.5831 

9.064513 

4.71  "l  6 

42.7084 

9.067591 

4.7116 

42.7229 

9.064775 

4.7116 

42.7096 

9.065146 

4.7116 

42.7113 

9.052497 

5.9623 

53.9737 

9.009426 

4.1013 

36.9504 

9.018324 

3.6995 

33.3633 

9.011096 

4.0475 

36.4724 

9.059747 

3.1649 

28.6732 

8.951976 

2.9062 

26.0162 

9.207458 

2.5541 

23.5168 

9.118356 

2.0339 

18.5458 

9.4452 

5.7346 

54.1644 
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#  MuSES  ascii  Results  File 

#  Exported  from  model:  c:\Program  Files\MuSES\Work\Examples\tez1 

#  Tue  Jun  1 3  00:22:47  2000 

#  Units  in  the  file:  SI  Units 
Radiance  Band:  10  to  11  micron 

Time  Step  10,  Time  =  Jul  19, 1984  18:00:00 


Element 

(N  u  mber&Locati  on  i  n 
the  Model) 

Radiance 

(W/mA2-sr) 

Area(m2) 

Eff.  Intensity 
(W/sr) 

8.356647 

1.2388 

10.3522 

wmmssm 

8.385643 

4.1946 

Ml  llli  III 

101  (front) 

8.364804 

3.798 

31.7695 

BEKBM 

8.418959 

1.512 

12.7295 

I USE 

8.423718 

1.729 

14.5646 

% — 

99  (front) 

8.363135 

1.373 

11.4826 

lliltHBI 

8.386425 

3.537" 

29.6628 

EMlBHiBB 

7.95239 

2.776 

22.0758 

9 

7.968305 

4.1102 

32.7513 

10 

98  (front) 

8.389138 

2.459 

20.6289 

8.418777 

4.3818" 

36.8894 

8.419935 

4.4529 

37.4931 

97  (front)  I 

8.420107 

0.3109 

2.6178 

7.97499 

3.928 

31.3258 

7.984747" 

5.806 

46.3594 

8.42214 

6.6406 

55.9281 

8.520856 

0.8872 

7.5597 

8.522503 

1.583 

13.4911 

92  (front) 

8.422995 

4.2047 

35.4162 

KCIHB 

8.424973 

3.1853 

26.8361 

E9BI 

8.516693 

2.3744 

20.2220 

22 

8.515099 

3.167 

26.9673 

23 

5IJJEGM 

8.439624 

1.3385 

11.2964 

mmmn 

8.430789 

1.7346 

14.6240 

IIS 

8.517737 

3.558 

30.3061 

26 

KEinHSBB 

8.534717 

4.2141 

35.9662 

l^fi 

8.446896 

0.4487 

3.7901 

28 

8.54843 

4.8426 

41.3966 

29 

82  (front) 

8.55274 

2.2895 

19.5815 

8.554167 

4.8426 

41.4244 

KBB 

9.397755 

2.0937 

19.6761 

HLi&H/HTTCI 

9.421009 

2.0938 

19.7257 

Kiiinsissi 

8.554566 

4.8427 

41.4272 

8.55461 

4.8426 

41.4266 

9.42607 

2.0938 

19.7363 

36 

1 1 4  (front) 

9.423661 

2.0938 

19.7313 

\wm 

8.55454 

4.8426 

41.4262 

m 

8.554604 

4.8426 

41.4265 

115 


fefcEIMH 


feaiMi 


TOTAL 


9.421999 

2.0938 

19.7278 

9.400816 

2.0938 

19.6834 

8.55398 

4.8426 

41.4235 

8.554455 

4.8426 

41.4258 

9.340044 

2.0949 

19.5665 

9.319592 

2.0937 

19.5124 

8.553594 

4.8426 

41.4216" 

8.553448 

4.8426 

41.4209 

9.336844 

1.9801 

18.4879 

9.391642 

2.2074 

20.7311 

8.55442" 

4.8426 

41 .4256 

8.554662 

4.8426 

41.4268 

9.417978 

2.0937 

19.7184 

9.383682 

2.1922 

20.5709 

8.554516 

4.8426 

41.4261 

8.553838 

4.8426 

41.4228 

9.336308 

1.9953 

18.6287 

9.294678 

2.0937 

19.4603 

8.553166 

4.8426 

41.4196 

8.552858 

4.8426" 

41.4181 

9.279945 

2.0937 

19.4294 

9.274809 

2.0937 

19.4187 

8.552723 

4.8426 

41.4174 

8.552625 

4.8426 

41.4169 

9.264443 

2.0937 

19.3970 

9.206858 

1.563 

14.3903 

8.551907 

4.8426 

41.4135 

8.558875 

4.8426 

41.4472 

8.580549 

4.873 

41.8130 

8.456082 

1.8129] 

15.3300 

8.578346 

5.9936 

51.4152 

8.931857 

6.2584 

55.8991 

8.95279 

4.7116 

42.1820 

8.955557 

4.7116 

42.1950 

8.953047" 

4.7116 

42.1832 

8.953365 

4.7116 

42.1847 

8.942036 

5.9623 

53.3151 

8.901906 

4.1013 

36.5094 

8.909859 

3.6995 

32.9620 

8.903377 

4.0475 

36.0364 

8.948543 

3.1649 

28.3212 

8.851626 

2.9062 

25.7246 

9.080413 

2.5541 

23.1923 

9.000957" 

2.0339 

18.3070 

9.2932 

5.7346 

53.2928 
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#  MuSES  ascii  Results  File 

#  Exported  from  model:  c:\Program  Files\MuSES\Work\Examples\tez1 

#  TueJun  13  00:22:47  2000 

#  Units  in  the  file:  SI  Units 
Radiance  Band:  11  to  12  micron 

Time  Step  10,  Time  =  Jul  1 9, 1 984  1 8:00:00 


Element 

(Number&Location  in 
the  Model) 

Radiance 

(W/mA2-sr) 

Area(m2) 

Eff.lntensity 

(W/sr) 

8.003191 

1.2388 

9.9144 

■HUHiiqillil 

8.028575 

4.1946 

33.6767 

mil 

8.010359 

3.798 

30.4233 

8.064414 

1.512 

12.1934 

5 

mi 

8.068616 

1.729 

13.9506 

6 

99  (front)  7 

8.0089 

1.373 

10.9962 

nH 

IliIifilflSlflM 

8.029258 

3.537 

28.3995 

8 

EIEIZKSI 

7.62776 

2.776 

21.1747 

9 

7.641724 

4.1102 

31.4090 

10 

PiUnHISB 

8.03163 

2.459 

19.7498 

i^R3HJ5B 

8.064289 

4.3818 

35.3361 

12 

8.065336 

4.4529 

35.9141 

13 

97  (front) 

8.058656 

0.3109 

2.5054 

7.647551 

3.928 

30.0396 

15 

wmznrm 

7.656176 

5.806 

44.4518 

16 

8.067234 

6.6406 

53.5713 

EEHB 

8.144329 

0.8872 

7.2256 

18 

8.145765 

1.583 

12.8947 

19 

92  (front) 

8.067982 

4.2047 

33.9234 

EBI 

8.069715 

3.1853 

25.7045 

8.140698 

2.3744 

19.3293 

22 

8.139308 

3.167 

25.7772 

23 

90  (front) 

8.082613 

1.3385 

10.8186 

i 

8.074909 

1.7346 

14.0067 

25 

mmmsm 

8.141609 

3.558 

28.9678 

26 

KHU'ATiUB 

8.1563 

4.2141 

34.3715 

KiEiSBiBI 

8.089013 

0.4487 

3.6295 

28 

8.16837 

4.8426 

39.5561 

29 

im 

8.172127 

2.2895 

18.7101 

KflflHTOT 

8.17337 

4.8426 

39.5804 

H 

8.909106 

2.0937 

18.6530 

32 

116  (front) 

8.929345 

2.0938 

18.6963 

33 

80  (front) 

8.173718 

4.8427 

39.5829 

34 

79  (front) 

8.173758 

4.8426 

39.5822 

■fPVBSSEl 

8.933751 

2.0938 

18.7055 

mvosissi 

8.931651 

2.0938 

18.7011 

Em 

8.173696 

4.8426" 

39.5819 

8.173752 

4.8426 

39.5822 
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^SHBEKSBSI 

8.930202 

2.0938 

18.6981 

BaKKgjS! 

8.911756 

2.0938 

18.6594 

8.173208 

4.8426 

39.5796 

nii-lmB 

8.173621 

4.8426 

39.5816 

Esmoai 

8.8589 

2.0949 

18.5585 

^■■^obisbssi 

8.841086 

2.0937 

18.5106 

RiMElllHliW 

8.1728 

4.8426 

39.5776 

8.172744 

4.8426 

39.5773 

47 

109  (front) 

8.85615 

1.9801 

17.5361 

8.903842 

2.2074 

19.6543 

49 

72  (front)  j 

8.173592 

4.8426 

39.5814 

50 

71  (front)  ; 

8.173801 

4.8426 

39.5824 

51 

222  (front) 

8.926758 

2.0937 

18.6900 

52 

108  (front) 

8.896885 

2.1922 

19.5038 

l^fl^HRiVDSISSB 

8.173675 

4.8426 

39.5818 

8.173083 

4.8426 

39.5790 

55 

107  (front) 

8.855595 

1.9953 

17.6696 

56 

106  (front) 

8.81935 

2.0937 

18.4651 

57 

68  (front)  ; 

8.172499 

4.8426 

39.5761 

58 

67  (front) 

8.17223 

4.8426 

39.5748 

59 

105  (front) 

8.806512 

2.0937 

18.4382 

60 

104  (front) 

8.802034 

2.0937 

18.4288 

61 

66  (front) 

8.172112 

4.8426 

39.5743 

62 

65  (front) 

8.172026 

4.8426 

39.5739 

63 

225  (front) 

8.793005 

2.0937 

18.4099 

64 

226  (front) 

8.742841 

1.563 

13.6651 

65 

64  (front) 

8.1714 

4.8426 

39.5708 

66 

8.177474 

4.8426 

39.6002 

67 

62  (front) 

8.196365 

4.873 

39.9409 

68 

234  (front) 

8.093132 

1.8129 

14.6720 

69 

8.194445 

5.9936 

49.1142 

ga^KgiBj 

8.502962 

6.2584 

53.2149 

8.521213 

4.7116 

40.1485 

72 

121  (front) 

8.523636 

4.7116 

40.1600 

8.521456 

4.7116 

40.1497 

Bum 

8.521722 

4.7116 

40.1509 

tj— iiirai 

8.511845 

5.9623 

50.7502 

bhb 

8.475698 

4.1013 

34.7614 

W2KHKz££XSE5B% 

8.482619 

3.6995 

31.3814 

78 

126  (front) 

8.476961 

4.0475 

34.3105 

79 

25  (front) 

8.517634 

3.1649 

26.9575" 

80 

125  (front) 

8.432584 

2.9062 

24.5068 

81 

124  (front) 

8.632243 

2.5541 

22.0476 

82 

1 1  (front) 

8.563214 

2.0339 

17.4167 

83 

30(front) 

8.8182 

5.7346 

50.5688 
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2.  VARIATION  IN  RADIATION  FROM  SHIP  TARGET  SURFACE  USED  IN 

muses  program 

The  Radiance  values  on  the  scale  show  the  reference  value. 


140.0 


10.01 


_ 34.4  H/m*2-sr 

mW  HMTH0.0 


31.5  M/ii"Z-sr 


11:00  AM 


3:00  PM 


5:00  PM 


Wm~' 


7:00  PM 


29.2  H/nTZ-sr 


H0.0 


10.0| 


9:00  PM 


25.8  H/mAZ-sr 


23:00  PM 
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APPENDIX  F 


SECOND  GENERATION  FLIR  SENSOR  PARAMETERS 

U.S.  Army  CECOM  NVESD  FLIR92  Thu  Hay  18  23:24:55  2000 

output  file:  SADAII . 1  short  listing 

data  file:  SADAII 

command  line  arguments:  -d  SADAII  -o  SADAII  -p  MRT 
begin  data  file  listing  ...  . 

gen2:  sample  data  file  for  RISTA  II  2nd  generation  FLIR  with  SADA  II  FPA 


>spectral 


spectral_cut__on 

8.0 

microns 

spectral__cut_of  f 

12.0 

microns 

diffract! on_wave 1 eng th 

0.0 

microns 

>optics_l 

f__number 

0.0 

— 

ef f_f  ocal_length 

20.0 

cm 

ef  f__aperture__diameter 

10.0 

cm 

optics_blur_spot 

0.01 

mrad 

average_optical_trans 

0.7 

>optics__2 

HFOV: VFOV_aspect_ratio 

0.0 

magnification 

0.0 

— 

frame_rate 

30.0 

Hz 

fields_per_frame 

1.0 

—  — 

>detector 

hor  z_dimens  ion__ ( active ) 

25.4 

microns 

ver  t__dimens  ion_  ( active ) 

25.4 

microns 

peak_D_star 

1 . 5el0 

cm-sqrt (Hz) /W 

integration__time 

0.007 

microseconds 

1/  f_knee_f  requency 

3.0 

Hz 

>fpa_parallel 

#_detectors_in_TDI 

4.0 

— 

#_vert_detectors 

480.0 

— 

#_j3amples_j?er_HIFOV 

2.0 

— 

#_samples_per_VIFOV 

2.0 

■ 

3dB_response_frequency 

2032.0 

Hz 

scan_ef f iciency 

0.75 

— 

>electronics 

hi gh_pas s_3 db_cu t on 

1.0 

Hz 

high_pas  s_f il t er_order 

0.0 

■“  - 

low_pas  s__3  db_cuto  f  f 

100000.0 

Hz 

low_pass_f ilter_order 

0.0 

— — 

boos  t_amp 1 i t ude 

0.0 

— 

boost_frequency 

0.0 

Hz 

s  amp 1 e_and_ho 1 d 

HORZ  NO_ 

_HORZ__VERT 

>display 

display_brightness 

10.0 

milli-Lamberts 

display_height 

15.24 

cm 

d i s p 1 ay _v i e wi ng_di s  t anc e 

30.0 

cm 

>  c  r  t__di  sp  1  ay 

#_ac  t  i  ve_l  ine  s__on__CRT 

480.0 

hor  z_crt_spo  t_s  igma 

0.0 

mrad 

vert_crt_spot_sigma 

0.0 

mrad 

>eye 

t hr e s ho 1 d_SNR 

2.5 

“ — 

eye_integration_time 

0.1 

sec 

MTF 

EXP 

EXP_or_NL 

>  random_image_rno  t  i  on 

hor  z_rms_mo  t  i  on_amp  1  i  t  ude 

0.02 

mrad 

ver  t_rms_mo  t  i  on_amp  1  i  t  ude 

0.02 

mrad 

>sinusoidal_image_motion 

hor  z_rms„mo  t  i  on_ampl  i  tude 

0.0 

mrad 

ve  r  t_rms_mo  t  i  on_amp  1  i  tude 

0.0 

mrad 

>3d  noise  default 

noise_level 

MOD 

NO_LO_MOD_or_HI 

>spectral_detectivity 

#_points :  9  microns^ 

_ detectivity 

8.00 

0.666 

8.50 

0.708 

9.00 

0.750 

9.50 

0.792 

10.00 

0.833 

10.50 

0.875 

11.00 

0.917 

11.50 

0.958 

12.00 

1.00 

>end 

end  data  file  listing  .  .  . 

MESSAGES 

diagnostic () :  Using  default  3D  noise  components, 
diagnostic () :  Using  _MOD_  level  3D  noise  defaults. 

diagnostic () :  Diffraction  wavelength  set  to  spectral  band  midpoint, 
diagnostic () :  HFOVrVFOV  aspect  ratio  defaulted  to  1.33. 
diagnostic () :  Fields-of-view  calculated  by  model. 

diagnostic ( ) :  Electronics  high  pass  filter  defaulted  to  order  1 . 
diagnostic  ()  :  Electronics  low  pass  filter  defaulted  to  order  1. 


CALCULATED  SYSTEM  PARAMETERS 

f ield-of-view:  2.323h  x  1.746v  degrees 

40 . 54h  x  30 . 48v  mrad 

magnification:  16.323 

optics  blur  spot:  48.800  microns  (dif fraction- limited) 

0.244  mrad 

detector  IFOV:  0.127h  x  0.127v  mrad 

scan  velocity:  1621.29  mrad/second 

dwell  time:  7.833e-005  seconds 

TEMPERATURE  DEPENDENCE 

parameter  NETD  §  300  K  NETD  @  OK  noise  bandwidth 


white  NETD  0.185  deg  C  0.000  deg  C  1.003e+004  Hz 
classical  NETD  0.185  deg  C  0.000  deg  C  1.007e+004  Hz 

sigma_TVH  NETD  0.103  deg  C  0.000  deg  C  3.134e+003  Hz 
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sigma_TV  NETD  0.077  deg  C  0.000  deg  C 

sigma_V  NETD  0.077  deg  C  0.000  deg  C 

Planck  integral  1.978e-004  0.000e+000  W/ (cm*cm*K) 

w/D-star  2.439e+006  0.000e+000  sgrt (Hz) / (cm*K) 


TOTAL  HORIZONTAL  MTFs 


cy/mr 

H_SYS 

H_PRE 

H_TPF 

H_SPF 

0.000 

1.000 

1.000 

1.000 

1.000 

0.394 

0.875 

0.945 

0.954 

0.971 

0.787 

0.697 

0.880 

0.847 

0.936 

1.181 

0.527 

0.808 

0.728 

0.897 

1.575 

0.388 

0.731 

0.623 

0.853 

1.969 

0.282 

0.651 

0.537 

0.806 

'  2.362 

0.202 

0.570 

0.468 

0.757 

2.756 

0.143 

0.491 

0.414 

0.706 

3.150 

0.100 

0.416 

0.369 

0.653 

3.543 

0.069 

0.345 

0.333 

0.601 

3.937 

0.046 

0.280 

0.303 

0.549 

4.331 

0.031 

0.222 

0.277 

0.498 

4.724 

0.020 

0.171 

0.256 

0.449 

5.118 

0.012 

0.127 

0.237 

0.401 

5.512 

0.007 

0.091 

0.221 

0.356 

5.906 

0.004 

0.063 

0.207 

0.314 

6.299 

0.002 

0.040 

0.194 

0.275 

6.693 

0.001 

0.024 

0.183 

0.239 

7.087 

0.000 

0.012 

0.173 

0.206 

7.480 

0.000 

0.004 

0.164 

0.176 

7.874 

0.000 

0.000 

0.156 

0.149 

TOTAL  VERTICAL 
cy/mr 

MTFs 

H_SYS 

H_PRE 

H_SPF 

0.000 

1.000 

1.000 

1.000 

0.394 

0.918 

0.945 

0.972 

0.787 

0.827 

0.880 

0.940 

1.181 

0.731 

0.808 

0.905 

1.575 

0.634 

0.731 

0.867 

1.969 

0.538 

0.651 

0.827 

2.362 

0.448 

0.570 

0.785 

2.756 

0.365 

0.491 

0.742 

3.150 

0.290 

0.416 

0.698 

3.543 

0.225 

0.345 

0.654 

3.937 

0.171 

0.280 

0.610 

4.331 

0.125 

0.222 

0.566 

4.724 

0.089 

0.171 

0.523 

5.118 

0.061 

0.127 

0.481 

5.512 

0.040 

0.091 

0.440 

5.906 

0.025 

0.063 

0.401 

6.299 

0.015 

0.040 

0.363 

6.693 

0.008 

0.024 

0.328 

7.087 

0.004 

0.012 

0.295 

7.480 

0.001 

0.004 

0.263 

7.874 

0.000 

0.000 

0.235 
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PREFILTER  VALUES  AT  NYQUIST 

horz  H_PRE (7.87)  =  0.000  vert  H_PRE(7.87)  =  0.000 

SAMPLING  RATES 

horizontal  15.75  samples/mr 
vertical  15.75  samples/mr 
effective  15.75  samples/mr 

SENSOR  LIMITING  FREQUENCIES 

spatial  Nyquist 

horizontal  7.87  7.87 

vertical  7.87  7.87 

effective  7.87  7.87 

MRTD  3D  NOISE  CORRECTION  (AVERAGE) 

300  K  OK 

horizontal  1.000  0.000 

vertical  3.833  0.000 


MRTD  AT 

300  K  BACKGROUND 

TEMPERATURE 

2D 

cy/mr 

horz 

cy/mr 

vert 

cy/mr 

0 . 05 

0.394 

0.007 

0.05 

0.394 

0.065 

0.830 

0.065 

0 . 10 

0.787 

0.017 

0.10 

0.787 

0.100 

1.120 

0.085 

0 . 15 

1.181 

0.031 

0.15 

1.181 

0.134 

1.424 

0.110 

0.20 

1.575 

0.053 

0.20 

1.575 

0.176 

1.767 

0.144 

0.25 

1.969 

0.085 

0.25 

1.969 

0.228 

2.112 

0.188 

0.30 

2.362 

0.133 

0.30 

2.362 

0.295 

2 . 458 

0.245 

0.35 

2.756 

0.206 

0.35 

2.756 

0.385 

2.794 

0.319 

0.40 

3.150 

0.318 

0.40 

3.150 

0.508 

3.116 

0 . 416 

0.45 

3.543 

0.494 

0.45 

3.543 

0.680 

3.425 

0.542 

0 . 50 

3.937 

0.778 

0.50 

3.937 

0.929 

3.720 

0.706 

0 . 55 

4.331 

1.245 

0.55 

4.331 

1.301 

4.001 

0.921 

0.60 

4.724 

2.038 

0.60 

4.724 

1.872 

4.268 

1.200 

0.65 

5.118 

3.433 

0.65 

5.118 

2.785 

4.522 

1 . 564 

0.70 

5.512 

6.000 

0.70 

5.512 

4.314 

4.767 

2.039 

0.75 

5.906 

10.815 

0.75 

5.906 

7.021 

4.998 

2.658 

0.80 

6.299 

21.134 

0.80 

6.299 

12.206 

5.219 

3.464 

0.85 

6.693 

45.254 

0.85 

6.693 

23.284 

5.429 

4 . 515 

0.90 

7.087 

99.999 

0.90 

7.087 

51.641 

5.629 

5.886 

0.95 

7.480 

99.999 

0.95 

7.480 

99.999 

5.820 

7.672 

1.00 
FLIR92 . 

7.874  99.999  1-00 

.  .  SADAII.l:  end  of  listing 

7.874 

99.999 

6.003 

10.000 
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APPENDIX  G 


MRTD  VALUES  and  CONVERSION  TABLE 


Spatial 

Frequecy 

(cyc/mrad) 

Range  Values  For 
Different  Tasks  (km) 

MRTDJ2D 

MRTD(With 

pol.filter) 

N=0.75 

N=3 

N=6 

0.161 

5.30441 

1.3261 

0.663 

0.009 

0.011 

0.204 

6.72112 

1.6803 

0.84 

0.012 

0.014 

0.258 

8.50024 

2.1251 

1.063 

0.016 

0.019 

0.325 

10.7077 

2.6769 

1.338 

0.021 

0.025 

0.411 

13.5411 

3.3853 

1.693 

0.028 

0.033 

0.52 

17.1323 

4.2831 

2.142 

0.037 

0.044 

0.656 

21.613 

5.4033 

2.702 

0.049 

0.058 

0.83 

27.3457 

6.8364 

3.418 

0.065 

0.076 

1.12 

36.9003 

9.2251 

4.613 

0.085 

0.100 

1.424 

46.9161 

11.729 

5.865 

0.11 

0.129 

1.767 

58.2168 

14.554 

7.277 

0.144 

0.169 

2.112 

69.5834 

17.396 

8.698 

0.188 

0.221 

2.458 

80.9829 

20.246 

10.12 

0.245 

0.288 

2.794 

92.053 

23.013 

11.51 

0.319 

0.375 

3.116 

102.662 

25.665 

12.83 

0.416 

0.489 

3.425 

112.842 

28.211 

14.11 

0.542 

0.638 

3.72 

122.562 

30.64 

15.32 

0.706 

0.831 

4.001 

131.82 

32.955 

16.48 

0.921 

1.084 

4.268 

140.616 

35.154 

17.58 

1.2 

1.412 

4.522 

148.985 

37.246 

18.62 

1.564 

1.840 

4.767 

157.057 

39.264 

19.63 

2.039 

2.399 

4.998 

164.667 

41.167 

20.58 

2.658 

3.127 

5.219 

171.949 

42.987 

21.49 

3.464 

4.075  ! 

5.429 

178.867 

44.717 

22.36 

4.515 

5.312 

5.629 

185.457 

46.364 

23.18 

5.886 

6.925 

5.82 

191.75 

47.937 

23.97 

7.672 

9.026 

6.003 

197.779 

49.445 

24.72 

10 

11.765 
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APPENDIX  H 


CONVERSION  OF  INBAND  RADIANCE  INTO  TEMPERATURE 

1.  MATLAB  CODE  FOR  CONVERSION 

%*******prograin  Calculates  the  Exitance  of  Given 
Temperature ( s ) ****** 
clc ; 

Tmin  =  input ('Enter  the  initial  temperature  you  want  to 
calculate=  ' ) ; 

Tmax  =  input ('Enter  the  final  temperature  you  want  to 
calculate=  ' )  ; 

Emit  =  input ('Enter  the  emittance=  '); 

h=6 . 626e-34 ; 

c=2 . 9  9  8e8 ; 

k=1.38e-23; 

dt=0.001; 

%change  the  bandwidth  if  you  want 

low=8e-6;  %**  lower  limit  of  band** 
up=12e-6;  %**  upper  limit  of  band** 
dx=0 . 001e-6 ; 
m=0  ; 
n=0  ; 

for  T=Tmin : dt : Tmax 
m=m+l ; 
tot=0 . 0 ; 

for  x=low : dx : up-dx 
xm=x+dx/2 ; 

%planck  formula 
fl=(2*h*cA2)/(xmA5) ; 
f 2=1/ (exp( (h*c) / (xm*k*T) ) -1) ; 
fx=f l*f 2*Emit ; 
tot= (tot+fx*dx) ; 
end; 

fprintf ( 'Temp:%d  Radiance :%e\n' ,T, tot) 
end; 
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2.  LOOK  UP  TABLE  FOR  CONVERSION  OF  INBAND  RADIANCE  INTO  TEMPERATURE 


T(K°) 

Rad.(W/m2st) 

T(K°) 

Rad.(W/m2st) 

T(K°) 

Rad.(W/m2st) 

T  (K°) 

Rad.(W/m2st) 

245.00 

12.9028 

245.48 

13.0524 

245.97 

13.2064 

246.46 

1 3.3615 

245.01 

12.9059 

245.49 

13.0556 

245.98 

13.2095 

246.47 

13.3647 

245.02 

12.9090 

245.5 

13.0587 

245.99 

13.2127 

246.48 

13.3679 

245  03 

12.9122 

245.51 

13.0618 

246 

13.2158 

246.49 

13.3711 

245.04 

12.9153  i 

245.52 

13.0650 

246.01 

13.2190 

246.50 

13.3743 

245  05 

12.9184 

245.53 

13.0681 

246.02 

13.2222 

246.51 

13.3774 

245.06 

12.9215 

245.54 

13.0712 

246.03 

13.2253 

246.52 

13.3806 

245  07 

12.9246 

245.55 

13.0744 

246.04 

13-2285 

246.53 

13.3838  < 

245.08 

12.9277 

245.56 
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251.39 


251.40 


251.41 


251.42 


251.43 


251.44 


251.45 


251.46 


251.47 


251.48 


251.49 


253.01 


253.02 


253.03 


253.04 


253.05 


253.06 


253.07 


253.08 


253.09 


253.1 


253.11 


253.12 


253.13 


253.14 


253.15 


253.16 


253.17 


253.18 


253.19 


253.2 


253.21 


253.22 


253.23 


253.24 


253.25 


253.26 


253.27 


253.28 


253.29 


253.3 


253.31 


253.32 


253.33 


253.34 


253.35 


253.36 


253.37 


14.9914 


14.9949 


14.9983 


15.0017 


15.0052 


15.0086 


15.0120 


15.0155 


15.0189 


15.0224 


15.0258 


15.5547 


15.5583 


15.5618 


15.5653 


15.5688 


15.5724 


15.5759 


15.5794 


15.5829 


15.5864 


15.5900 


15.5935 


15.5970 


15.6006 


15.6041 


15.6076 


15.6111 


15.6147 


15.6182 


15.6217 


15.6253 


15.6288 


15.6323 


15.6359 


15.6394 


15.6429 


15.6465 


15.6500 


15.6535 


15.6571 


15.6606 


15.6641 


15.6677 


15.6712 


15.6747 


15.6783 


15.6818 


251.89 


251.9 


251.91 


251.92 


251.93 


251.94 


251.95 


251.96 


251.97 


251.98 


251.99 


253.52 


253.53 


253.54 


253.55 


253.56 


253.57 


253.58 


253.59 


253.6 


253.61 


253.62 


253.63 


253.64 


253.65 


253.66 


253.67 


253.68 


253.69 


253.7 


253.71 


253.72 


253.73 


253.74 


253.75 


253.76 


253.77 


253.78 


253.79 


253.8 


253.81 


253.82 


253.83 


253.84 


253.85 


253.86 


253.87 


253.88 


15.1638 


15.1673 


15.1707 


15.1742 


15.1777 


15.1811 


15.1846 


15.1880 


15.1915 


15.1950 


15.1984 


15.7350 


15.7385 


15.7421 


15.7456 


15.7492 


15.7527 


15.7563 


15.7598 


15.7634 


15.7669 


15.7705 


15.7740 


15.7776 


15.7811 


15.7847 


15.7882 


15.7918 


15.7954 


15.7989 


15.8025 


15.8060 


15.8096 


15.8131 


15.8167 


15.8203 


15.8238 


15.8274 


15.8309 


15.8345 


15.8381 


15.8416 


15.8452 


15.8488 


15.8523 


15.8559 


15.8595 


15.8630 


252.39 


254.10 


254.11 


254.12 


254.13 


254.14 


254.15 


254.16 


254.17 


254.18 


254.19 


254.20 


254.21 


254.22 


254.23 


254.24 


254.25 


254.26 


254.27 


254.28 


254.29 


254.30 


254.31 


254.32 


254.33 


254.34 


254.35 


254.36 


254.37 


254.38 


254.39 


15.3410 


15.3445 


15.3480 


15.3515 


15.3550 


15.3584 


15.3619 


15.3654 


15.3689 


15.3724 


15.9166 


15.9202 


15.9237 


15.9273 


15.9309 


15.9345 


15.9380 


15.9416 


15.9452 


15.9488 


15.9523 


15.9559 


15.9595 


15.9631 


15.9667 


15.9703 


15.9738 


15.9774 


15.9810 


15.9846 


15.9882 


15.9918 


15.9953 


15.9989 


16.0025 


16.0061 


16.0097 


16.0133 


16.0169 


16.0205 


16.0241 


16.0276 


16.0312 


16.0348 


16.0384 


16.0420 


16.0456 


252.90 


252.91 


252.92 


252.93 


252.94 


252.95 


252.96 


252.97 


252.98 


252.99 


254.53 


254.54 


254.55 


254.56 


254.57 


254.58 


254.59 


254.6 


254.61 


254.62 


254.63 


254.64 


254.65 


254.66 


254.67 


254.68 


254.69 


254.7 


254.71 


254.72 


254.73 


254.74 


254.75 


254.76 


254.77 


254.78 


254.79 


254.8 


254.81 


254.82 


254.83 


254.84 


254.85 


254.86 


254.87 


254.88 


254.89 


15.5161 


15.5196 


15.5231 


15.5266 


15.5301 


15.5336 


15.5372 


15.5407 


15.5442 


15.5477 


16.0960 


16.0996 


16.1032 


16.1068 


16.1104 


16.1140 


16.1176 


16.1212 


16.1248 


16.1284 


16.1320 


16.1356 


16.1392 


16.1428 


16.1464 


16.1501 


16.1537 


16.1573 


16.1609 


16.1645 


16.1681 


16.1717 


16.1753 


16.1789 


16.1826 


16.1862 


16.1898 


16.1934 


16.1970 


16.2006 


16.2043 


16.2079 


16.2115 


16.2151 


16.2187 


16.2224 


16.2260 


253.38 

15.6854 

253.89 

15.8666 

254.40 

16.0492 

254.9 

16.2296 

253.39 

15.6889  1 

253.9 

15.8702 

254.41 

16.0528 

254.91 

16.2332 

253.4 

15.6924 

253.91 

15.8737 

254.42 

16.0564 

254.92 

16.2368 

253.41 

15.6960 

253.92 

15.8773 

254.43 

16.0600 

254.93 

16.2405 

253.42 

15.6995 

253.93 

15.8809 

254.44 

16.0636 

254.94 

16.2441 

253.43 

15.7031 

253.94 

15.8844 

254.45 

16.0672 

254.95 

16.2477 

253.44 

15.7066 

253.95 

15.8880 

254.46 

16.0708 

16.2513 

253.45 

15.7102 

253.96 

15.8916 

254.47 

16.0744 

254.97 

16.2550 

253.46 

15.7137 

253.97 

15.8951 

254.48 

16.0780 

254.98 

16.2586 

253.47 

15.7172 

253.98 

15.8987 

254.49 

16.0816 

254.99 

16.2622 

253.48 

15.7208 

253.99 

15.9023 

254.50 

16.0852 

255 

16.2658 

253.49 

15.7243 

254 

15.9059 

254.51 

16.0888 

255 

16.2658 

253.5 

15.7279 

254.01 

15.9094 

254.52 

16.0924 

255.01 

16.2695 

253.51 

15.7314 

254.02 

15.9130 

254.53 

16.0960 

255.02 

16.2731 

255.03 

16.2767 

255.67 

16.5100 

256.18 

16.6975 

256.69 

16.8864 

255.04 

16.2803 

255.68 

16.5137 

256.19 

16.7012 

256.7 

16.8901  • 

255.05 

16.2840 

255.69 

16.5173 

256.2 

16.7049 

256.71 

16.8939 

255.06 

16.2876 

255.70 

16.5210 

256.21 

16.7086 

256.72 

16.8976 

255.07 

16.2912 

255.71 

16.5247 

256.22 

16.7123 

256.73 

16.9013 

255.08 

16.2949 

255.72 

16.5283 

256.23 

16.7160 

256.74 

.  16.9050 

255.09 

16.2985 

255.73 

16.5320 

256.24 

16.7197 

256.75 

16.9087 

255.1 

16.3021 

255.74 

16.5357 

256.25 

16.7234 

256.76 

16.9125 

255.11 

16.3058 

255.75 

16.5393 

256.26 

16.7271 

256.77 

16.9162 

255.12 

16.3094 

255.76 

16.5430 

256.27 

16.7307 

256.78 

16.9199 

255.13 

16.3130 

255.77 

16.5467 

256.28 

16.7344 

256.79 

16.9236 

255.14 

16.3167 

255.78 

16.5503 

256.29 

16.7381 

256.8 

16.9273 

255.15 

16.3203 

255.79 

16.5540 

256.3 

16.7418 

256.81 

16.9311 

255.16 

16.3239 

255.80 

16.5577 

256.31 

16.7455 

256.82 

16.9348 

255.17 

16.3276 

255.81 

16.5613 

256.32 

16.7492 

256.83 

16.9385 

255.18 

16.3312 

255.82 

16.5650 

256.33 

16.7529 

256.84 

16.9422 

255.19 

16.3348 

255.83 

16.5687 

256.34 

16.7566 

256.85 

16.9460 

255.2 

16.3385 

255.84 

16.5724 

256.35 

16.7603 

256.86 

16.9497 

255.21 

16.3421 

255.85 

16.5760 

256.36 

16.7640 

256.87 

16.9534 

255.22 

16.3457 

255.86 

16.5797 

256.37 

16.7677 

256.88 

16.9572 

255.23 

16.3494 

255.87 

16.5834 

256.38 

16.7714 

256.89 

16.9609 

255.24 

16.3530 

255.88 

16.5871 

256.39 

16.7751 

256.9 

16.9646 

255.25 

16.3567 

255.89 

16.5907 

256.4 

16.7788 

256.91 

16.9683 

255.26 

16.3603 

255.90 

16.5944 

256.41 

16.7825 

256.92 

16.9721 

255.27 

16.3639 

255.91 

16.5981 

256.42 

16.7862 

256.93 

16.9758 

255.28 

16.3676 

255.92 

16.6018 

256.43 

16.7899 

256.94 

16.9795 

255.29 

16.3712 

255.93 

16.6054 

256.44 

16.7936 

256.95 

16.9833 

255.3 

16.3749 

255.94 

16.6091 

256.45 

16.7973 

256.96 

16.9870 

255.31 

16.3785 

255.95 

16.6128 

256.46 

16.8011 

256.97 

16.9907 

255.32 

16.3822 

255.96 

16.6165 

256.47 

16.8048 

256.98 

16.9945 

255.33 

16.3858 

255.97 

16.6201 

256.48 

16.8085 

256.99 

16.9982 

255.34 

16.3894 

255.98 

16.6238 

256.49 

16.8122 

257 

17.0019 

255.35 

16.3931 

255.99 

16.6275 

256.5 

16.8159 

257.01 

17.0057 

255.36 

16.3967 

256.00 

16.6312 

256.51 

16.8196 

257.02 

17.0094 

255.37 

16.4004 

256.01 

16.6349 

256.52 

16.8233 

257.03 

17.0131 
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255.38 

16.4040 

256.02 

16.6385 

256.53 

16.8270 

257.04 

17.0169 

255.39 

16.4077 

256.03 

16.6422 

256.54 

16.8307 

257.05 

17.0206 

255.4 

16.4113 

256.04 

16.6459 

256.55 

16.8344 

257.06 

17.0243 

255.41 

16.4150 

256.05 

16.6496 

256.56 

16.8381 

257.07 

17.0281 

255.42 

16.4186 

256.06 

16.6533 

256.57 

16.8418 

257.08 

17.0318 

255.43 

16.4223 

256.07 

16.6570 

256.58 

16.8456 

257.09 

17.0356 

255.44 

16.4259 

256.08 

16.6606 

256.59 

16.8493 

257.10 

17.0393 

255.45 

16.4296 

256.09 

16.6643 

256.6 

16.8530 

257.11 

17.0430 

255.46 

16.4332 

256.1 

16.6680 

256.61 

16.8567 

257.12 

17.0468 

255.47 

16.4369 

256.11 

16.6717 

256.62 

16.8604 

257.13 

17.0505 

255.48 

16.4405 

256.12 

16.6754 

256.63 

16.8641 

257.14 

17.0543 

255.49 

16.4442 

256.13 

16.6791 

256.64 

16.8678 

257.15 

17.0580 

255.5 

16.4478 

256.14 

16.6828 

256.65 

16.8716 

257.16 

17.0618 

255.51 

16.4515 

256.15 

16.6864 

256.66 

16.8753 

257.17 

17.0655 

255.65 

16.5027 

256.16 

16.6901 

256.67 

16.8790 

257.18 

17.0692 

255.66 

16.5064 

256.17 

16.6938 

256.68 

16.8827 

257.19 

17.0730 

257.20 

17.0767 

257.7 

17.2647 

258.2 

17.4540 

258.70 

17.6447 

257.21 

17.0805 

257.71 

17.2685 

258.21 

17.4578 

258.71 

17.6485 

257.22 

17.0842 

257.72 

17.2722 

258.22 

17.4616 

258.72 

17.6524 

257.23 

17.0880 

257.73 

17.2760 

258.23 

17.4654 

258.73 

17.6562 

257.24 

17.0917 

257.74 

17.2798 

258.24 

17.4692 

258.74 

17.6600 

257.25 

17.0955 

257.75 

17.2836 

258.25 

17.4730 

258.75 

17.6638 

257.26 

17.0992 

257.76 

17.2873 

258.26 

17.4768 

258.76 

17.6677 

257.27 

17.1030 

257.77 

17.2911 

258.27 

17.4806 

258.77 

17.6715 

257.28 

17.1067 

257.78 

17.2949 

258.28 

17.4844 

258.78 

17.6753 

257.29 

17.1105 

257.79 

17.2987 

258.29 

17.4882 

258.79 

17.6792 

257.30 

17.1142 

257.8 

17.3024 

258.3 

17.4920 

258.80 

17.6830 

257.31 

17.1180 

257.81 

17.3062 

258.31 

17.4958 

258.81 

17.6868 

257.32 

17.1217 

257.82 

17.3100 

258.32 

17.4997 

258.82 

17.6907 

257.33 

17.1255 

257.83 

17.3138 

258.33 

17.5035 

258.83 

17.6945 

257.34 

17.1292 

257.84 

17.3176 

258.34 

17.5073 

258.84 

17.6983 

257.35 

17.1330 

257.85 

17.3213 

258.35 

17.5111 

258.85 

17.7022 

257.36 

17.1367 

257.86 

17.3251 

258.36 

17.5149 

258.86 

17.7060 

257.37 

17.1405 

257.87 

17.3289 

258.37 

17.5187 

258.87 

17.7098 

257.38 

17.1442 

257.88 

17.3327 

258.38 

17.5225 

258.88 

17.7137 

257.39 

17.1480 

257.89 

17.3365 

258.39 

17.5263 

258.89 

17.7175 

257.40 

17.1518 

257.9 

17.3403 

258.4 

17.5301 

258.90 

17.7214 

257.41 

17.1555 

257.91 

17.3440 

258.41 

17.5339 

258.91 

17.7252 

257.42 

17.1593 

257.92 

17.3478 

258.42 

17.5377 

258.92 

17.7290 

257.43 

17.1630 

257.93 

17.3516 

258.43 

17.5416 

258.93 

17.7329 

257.44 

17.1668" 

257.94 

17.3554 

258.44 

17.5454" 

258.94 

17.7367 

257.45 

17.1705 

257.95 

17.3592 

258.45 

17.5492 

258.95 

17.7406 

257.46 

17.1743 

257.96 

17.3630 

258.46 

17.5530 

258.96 

17.7444 

257.47 

17.1781 

257.97 

17.3668 

258.47 

17.5568 

258.97 

17.7482 

257.48 

17.1818 

257.98 

17.3705 

258.48 

17.5606 

258.98 

17.7521 

257.49 

17.1856 

257.99 

17.3743 

258.49 

17.5644 

258.99 

17.7559 

257.50 

17.1893 

258 

17.3781 

258.5 

17.5683 

259.00 

17.7598 

257.51 

17.1931 

258.01 

17.3819 

258.51 

17.5721 

259.01 

17.7636 

257.52 

17.1969 

258.02 

17.3857 

258.52 

17.5759 

259.02 

17.7675 

APPENDIX  I 


BACKGROUND-TARGET  RADIANCE  and  ATD  COMPUTATIONS  FOR  UNPOLARIZED  CASE 

The  method  to  calculate  the  field  values  of  this  table  is  given  in 
Chapter  5  Section  C. 


Zenith 

Angle 

(Degree) 

Range 

(km) 

Path 

Rad. 

Radiance  From 
Background* 
(W/m2-sr) 

Bandwidth 

(micrometer) 

Trans. 

Radiance  From 
Target* 
(W/m2-sr) 

ATD 

(C°) 

8-9 

0.8391 

7.014 

9-10 

0.9205 

8.070 

165 

0.31 

3.631 

31.9615 

10-11 

0.9225 

8.012 

3.79 

11-12 

0.8933 

7.403 

Total 

30.500 

8-9 

0.7618 

6.368 

9-10 

0.8693 

7.621 

120 

0.6 

5.842 

31.874 

10-11 

0.8678 

7.537 

4.00 

11-12 

0.8203 

6.798 

Total 

28.325 

8-9 

0.6513 

5.444 

9-10 

0.7855 

6.886 

105 

1.159 

9.353 

31.47 

10—11 

0.7759 

6.739 

5.86 

11-12 

0.7014 

5.813 

Total 

24.883 

8-9 

0.6378 

5.331 

9-10 

0.7745 

6.790 

104 

1.24 

9.805 

31.5145 

10-11 

0.7637 

6.633 

4.78 

11-12 

0.6861 

5.686 

Total 

24.441 

8-9 

0.6229 

5.207 

9-10 

0.7621 

6.681 

103 

1.333 

10.31 

31.4058 

10-11 

0.7499 

6.513 

5.00 

11-12 

0.6688 

5.543 

Total 

23.944 

8-9 

0.6061 

5.066 

9-10 

0.7479 

6.557 

102 

1.443 

10.89 

31.15807 

10-11 

0.7341 

6.376 

5.47 

11-12 

0.6492 

5.380 

Total 

23.379 

8-9 

0.6007 

5.021 

9-10 

0.7433 

6.516 

101.7 

1.479 

11.08 

31.1859 

10-11 

0.729 

6.332 

5.42 

11-12 

0.6429 

5.328 

Total 

23.197 

8-9 

0.5989 

5.006  [ 

9-10 

0.7417 

6.502 

101.6 

1.492 

11.14 

31.1959 

10-11 

0.7272 

6.316 

5.41 

11-12 

0.6407 

5.310 

Total 

23.135 
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101.5  1.505  11.2  31.2058 


101.4  1.518  11.27  31.2169 


101.3  1.531  11.34  31.0277 


101.2  1.545  11.41  31.0391 


101.1  1.558  31.051 


100.5  1.647  11.91  31.1369 


100.2  1.694  12.14  31.1893 


100  1.728  12.3  31.2284 


99.8  1.763  12.47  31.271 


8-9 


9-10 


10-11 


11-12 


Total 


8-9 


9-10 


10-11 


11-12 


Total 


8-9 


9-10 


10-11 


11-12 


Total 


8-9 


9-10 


10-11 


11-12 


Total 


8-9 


9-10 


10-11 


11-12 


Total 


8-9 


9-10 


10-11 


11-12 


Total 


8-9 


9-10 


10-11 


11-12 


Total 


8-9 


9-10 


10-11 


11-12 


Total 


8-9 


9-10 


10-11 


11-12 


Total 


0.597 


0.7401 


0.7254 


0.6385 


0.5951 


0.7384 


0.7236 


0.6362 


0.5932 


0.7368 


0.7218 


0.634 


0.5912 


0.7351 


0.7199 


0.6317 


0.5893 


0.7334 


0.718 


0.6293 


0.5769 


0.7226 


0.7059 


0.6146 


0.5703 


0.7167 


0.6995 


0.6068 


0.5658 


0.7127 


0.695 


0.6013 


0.5612 


0.7086 


0.6904 


0.5957 


4.990 


6.488 


6.300 


5.292 


23.071 


4.975 


6.473 


6.285 


5.273 


23.005 


4.959 


6.459 


6.269 


5.254 


22.942 


4.942 


6.444 


6.253 


5.235 


22.874 


4.926 


6.429 


6.236 


5.216 


22.807 


4.822 


6.335 


6.131 


5.094 


22.382 


4.767 


6.283 


6.076 


5.029 


22.155 


4.730 


6.248 


6.036 


4.983 


21.997 


4.691 


6.212 


5.996 


4.937 


21.837 
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1  8-9 _ 0.5126  4.285 

9-10  0.6638  5.819 

98  2.157  14.23  31.611  10-11 _ 0.6406 _ 5.564  4.77 

11-12  0.5365  4.446 

_ Total _ 20.114 _ 

8- 9 _ 0.4964 _ 4.149 

9- 10  0.6485 _ 5.685 

97.5  2.3  14.82  31.777  10-11  0.6236  5.416  4.51 

11-12  0.5167  4.282 

Total  19.533 
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97.4 


2.331 


14.95 


31.812 


4.44 


8- 9 _ 0.493 _ 4.121 

9- 10  0.6452  5.656 

10- 11  0.6199  5.384 

11- 12  0.5125  4.247 

Total  19.409 


8- 9 _ 0.4472  3.738 

9- 10  0.5999  5.259 

96.2  2.782  16.68  31.94008  10-11  0.5698 _ 4.949  4.28 

11-12  0.4555  3.775 

Total  17.721 
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96.1 


2.827 


16.84 


31.38594 


5.28 


8- 9 _ 0.4429 _ 3.702 

9- 10  0,5956  5,221 

10- 11 _ 0.565 _ 4.907 

11- 12  0.4501 _ 3.730 

Total  17.561 


8- 9  0.2987  2.497 

9- 10  0,4376 _ 3.836 

93.6  4.801  22.55  32.662  10-11  0,3933 _ 3.416  3.25 

11-12  0.2713  2.248 

_  Total  11.998 
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APPENDIX  J 


COMPUTATIONS  WITH  POLARIZED  FILTER 


1.  BREAKDOWN  of  BACKGROUND  RADIANCE 

When  the  horizontal  polarizer  is  added  to  the  sensor 
the  fraction  of  radiance  that  enters  the  sensor  will  vary 
depending  on  the  degree  of  polarization  value  (Dpol) . 


Range 

(km) 

Sea 

emission 

Dpol 

Value* 

Incident 

Rad.Value 

Sky 

Reflection 

Dpol 

Value** 

Incident 

Rad. 

Value 

Path  Rad. 

DPol(%0) 

Total 

0.31 

28.2065 

0.1 

14.089147 

0.1235 

11.6 

0.06891 

3.6314 

1.8157 

15.97376" 

0.319 

28.1378 

0.2 

14.040762 

0.1268" 

21 

0.07671 

3.70386 

1.85193 

15.96941 

0.599 

25.3283 

KB 

12.195576 

0.7033 

76.6 

0.62101 

5.8429 

2.92145" 

15.73804 

1.728 

14.9885 

6.7523193 

3.9385 

26.7 

2.49504 

12.3013 

6.15065 

15.39801 

2.157 

12.8075 

10.7 

5.7185488 

4.5758 

22.2 

2.79581 

14.2279 

7.11395 

15.62831 

2.874 

10.2973 

11.4 

4.5617039 

4.1221 

20.4 

2.4815 

17.011 

8.5055 

15.54871 

3.45 

8.7899 

11.8 

3.8763459 

4.1881 

18.3 

2.47726 

18.9227 

9.46135 

15.81496 

4.317 

7.1319 

12.2 

3.1309041 

3.8133 

17 

2.23078 

21 .369 

10.6845 

16.04618 

4.801 

6.4201 

12.4 

2.8120038 

2.1471 

16.2 

1.24747 

22.5454 

11.2727 

15.33217 

5.408 

5.6673 

12.6 

2.4766101 

2.0202" 

15.4 

1.16566 

23.8613 

11.93065 

15.57292 

5.773 

5.2531 

12.7 

2.2929782 

1.9306 

15 

1.1101 

24.5774 

12.2887 

15.69177 

6.677 

4.3621 

12.8 

1.9018756 

1.6999 

14.2 

0.97064 

26.137 

13.0685 

15.94102 

7.921 

3.4151 

13 

1.4855685 

1 .3839 

13.8 

0.78744 

27.875 

13.9375 

16.21051 

8.738 

2.9214" 

13.1 

1.2693483 

1.2167 

13.4 

0.68987 

28.8093 

14.40465 

16.36387 

9.748 

2.443 

13.3 

1.0589538 

1.0357 

13.1 

0.58569 

29.7811 

14.89055 

16.53519 

1 1 .03 

1.9244 

13.4 

0.8332652 

0.8429 

12.9" 

0.47582 

30.7803 

15.39015 

16.69923" 

12.71 

1.4354 

13.5" 

0.6208105 

0.643 

12.6 

0.36201 

31.7865 

15.89325 

16.87607 

15.03 

0.9709 

13.7 

0.4189434 

0.4436 

12.4 

0.2493 

32.7657 

16.38285 

17.0511 

18.47 

0.554 

13.8 

0.238774 

0.2574 

12.3 

0.14453 

33.6628 

16.8314 

17.2147 

20.92 

0.3752 

14 

0.161336 

0.1755 

12.31 

0.09854 

34.054 

17.027 

17.28688 

24.24 

0.224 

14.1 

0.096208 

0.1054 

12.3 

0.059181 

34.3916 

17.1958 

17.35119 

29.09 

0.1072 

14.3 

0.0459352 

0.0506 

12.3 

0.02841 

34.6583 

17.32915 

17.4035 

29.7T 

0.09771 

14.3 

0.0418687 

0.0462 

12.3 

0.02594 

34.68063 

17.340321 

17.40813 

30.36 

0.0886 

14.3 

0.037965l1 

0.0419 

12.3 

0.02353 

34.70212 

17.35106 

17.41255 

31.05 

0.07991 

14.3 

0.0342414 

0.0378 

12.3 

0.02122 

34.72279 

17.3614 

17.41686 

31.78 

0.07166 

14.4 

0.0306705 

0.0339 

12.3 

0.01903 

34.74259 

17.3713 

17.421 

32.56 

0.06383 

14.4 

0.0273192 

0.0302 

12.3 

0.01696 

34.76155 

17.38078 

17.42505 

33.39 

0.05644 

14.4 

0.0241563 

0.0267 

12.3 

0.01499 

34.77967 

17.38984 

17.42898 

34.28 

0.04949 

14.4 

0.0211817 

0.0234 

12.4 

0.01315 

34.79701 

17.39851 

17.43284 

35.23 

0.04299 

14.5 

0.0183782 

0.0203 

12.4 

0.01141 

34.81353 

17.40677 

17.43655 

36.26 

0.03694 

14.5 

0.0157919 

0.0174 

12.4 

0.00978 

34.82927 

17.41464 

17.44021 

37.38 

0.03134 

14.5 

0.0133979 

0.0148 

12.4 

0.00832 

34.84424 

17.42212 

17.44384 

38.6 

0.02618 

14.6 

0.0111789 

0.0124 

12.4 

0.00697 

34.85853 

17.42927 

17.44741 

*  All  the  Dpol  values  for  sea  surface  emission  are  vertical  polarization 

**  All  the  Dpol  values  for  Sky  reflection  component  are  horizontal  polarization 
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APPENDIX  K 


ATD  -  DELTA  W  COMPARISON 


Range 

(km) 

Target+path.ra 
d.value  (0%) 

Target+path.re 
d.value  (2%) 

Target+path.re 
d.value  (8%) 

Bckground 
Total  Radiance 

Delta  W  (0% 
Dpol) 

Delta  W  (2% 
Dpol.) 

Delta  W  (%8 
Dpol.) 

0.310 

17.0653 

17.3703 

18.2853 

15.9738 

1.0916 

1.3966 

2.3115 

0.599 

17.0833 

17.3665 

18.0339 

15.7380 

1 .3452 

1.6285 

2.2959 

1.728 

17.1494 

17.3693 

17.8419 

15.3980 

1.7513 

1.9713 

2.4439 

2.157 

17.1707 

17.3719 

17.9527 

15.6283 

1.5424 

1.7436 

2.3244 

2.874 

17.1862 

17.3752 

17.8599 

15.5487 

1.6375 

1.8264 

2.3112 

3.450 

17.2274 

17.3828 

17.8487 

15.8150 

1.4125 

1.5678 

2.0338 

4.317 

17.2574 

17.3889 

17.7832 

16.0462 

1.2112 

1.3427 

1.7371 

4.801 

17.2713 

17.3913 

17.7512 

15.3322 

1.9391 

2.0591 

2.4190 

5.408 

17.2884 

17.3955 

17.7170 

15.5729 

1.7154 

1.8226 

2.1441 

5.773 

17.2967 

17.3969 

17.6974 

15.6918 

1.6050 

1.7051 

2.0056 

6.677 

17.3177 

17.4027 

17.6577 

15.9410 

1.3767 

1.4617 

1.7166 

7.921 

17.3403 

17.4084 

17.6125 

16.2105 

1.1298 

1.1979 

1.4020 

8.738 

17.3526 

17.4115 

17.5884 

16.3639 

0.9887 

1.0477 

1.2245 

9.748 

17.3655 

17.4150 

17.5635 

16.5352 

0.8303 

0.8799 

1.0284 

11.029 

17.3792 

17.4190 

17.5384 

16.6992 

0.6800 

0.7198 

0.8391 

12.711 

17.3939 

17.4239 

17.5139 

16.8761 

0.5178 

0.5478 

0.6379 

15.030 

17.4082 

17.4287 

17.4902 

17.0511 

0.3571 

0.3776 

0.4391 

18.469 

17.4226 

17.4344 

17.4699 

17.2147 

0.2079 

6.2197 

0.2552 

20.922 

17.4297 

17.4377 

17.4619 

17.2869 

0.1428 

0.1508 

0.1750 

24.240 

17.4305 

17.4352 

17.4493 

17.3512 

0.0793 

0.0840 

0.0981 

29.087 

17.4457 

17.4480 

17.4550 

17.4035 

0.0422 

0.0445 

0.0515 

29.707 

17.5365 

17.5386 

17.5450 

17.4081 

0.1284 

0.1305 

0.1369 

30.360 

17.5285 

17.5304 

17.5362 

17.4126 

0.1160 

0.1179 

0.1236 

31.050 

17.5214 

17.5231 

17.5283 

17.4169 

0.1045 

0.1062 

0.1114 

31.782 

17.5146 

17.5162 

17.5208 

17.4210 

0.0936 

0.0952 

0.0998 

32.559 

17.5091 

17.5105 

17.5147 

17.4251 

0.0841 

0.0855 

0.0896 

33.388 

17.5032 

17.5044 

17.5081 

17.4290 

0.0742 

0.0754 

0.0791 

34.276 

17.4981 

17.4992 

17.5024 

17.4328 

0.0652 

0.0663 

0.0696 

35.229 

17.4935 

17.4944 

17.4972 

17.4366 

0.0569 

0.0578 

0.0607 

36.259 

17.4893 

17.4901 

17.4925 

17.4402 

0.0491 

0.0499 

0.0523 

37.378 

17.4851 

17.4858 

17.4879 

17.4438 

0.0413 

0.0420 

0.0440 

38.603 

17.4819 

17.4825 

17.4842 

17.4474 

0.0345 

0.0351 

0.0368 
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Radiance  Diffrence  (Watts/mA2  str) 


Temperature  Difference  (ATD)  &  Radiance  Difference  (Delta  W) 


144 


LIST  OF  REFERENCES 


1.  Cooper  A.W. ,  "Polarization  Issues  in  Infrared  Imaging  of 
Ships  Against  Background, "  Proceedings  of  the  DDR&E 
Electromagnetic /Electro-Optic  Performance  Prediction  and 
Products  Symposium,  Monterey,  1997 

2.  Cooper  A.W.  and  Crittenden  E.C.,  Electro-Optic  Sensors  and 
Systems,  Manuscript,  Naval  Postgraduate  School,  September 

1998 

3.  Kaplan  H. ,  Practical  Applications  of  Infrared  Thermal 
Sensing  and  Imaging  Equipment,  SPIE  Optical  Engineering 
Press,  Washington,  1993 

4.  Wolfe  L.  W.  and  Zissis  J.  G.,  The  Infrared  Handbook,  Office 
of  Naval  Research,  Department  of  Navy,  Washington,  DC,  1978 

5.  Khalil  S.  and  Hovanessian  S.A.,  Introduction  to  Electro- 
Optical  Imaging  and  Tracking  Systems,  Artech  House,  Inc., 
MA,  1993 

6.  Hackforth  L.H.,  Infrared  Radiation,  McGraw-Hill  Book 
Company,  Inc.,  New  York,  1960 

7.  Kruse  P.W.,  McGlauchlin  L.,and  McQuistan  R.B.,  Elements  of 
Infrared  Technology,  John  Willey  &  Sons,  Inc.  New  York, 
1962 


145 


Jamieson  J.A. ,  McFee  R.H.,  Plass  G.N.,  Grube  R.H.,  and 
Richards  R.G. ,  Infrared  Physics  and  Engineering,  McGraw- 
Hill  Book  Company,  1963 

Wollenweber  G.  F.,  "Thermal  Background  Modeling  and  Its  Use 
in  Science  and  Technology",  SPIE  Proceedings,  Vol.1687, 


1992 

10.  Zeisse  C.R.,  SeaRad,  A  Sea  Radiance  Prediction  Code, 
Technical  Report  1702,  Naval  Command,  Control,  and  Ocean 
Surveillance  Center,  RDT  and  E  Division,  November  1995 

11.  Accetta  J.S.,  Shumaker  D.L.,  The  Infrared  and  Electro- 
Optical  Systems  Handbook,  Vol.6,  SPIE  Optical  Engineering 
Press  and  Infrared  Information  Analysis  Center,  1996 

12.  Jordan  D.L.,  Lewis  G-,  "Infrared  Polarization  Signatures," 
AGARD  Meeting  on  Atmospheric  Propagation  Effects  Through 
Natural  and  Man-made  Obscurants  for  Visible  to  MM-Wave 


Radiation,  May  1993 

Morey  B.,  Ellis  K.  ,  Perry  D.  ,  and  Gleichman  K.  ,  Infrared 
Signature  Simulation  of  Military  Targets,  Infrared 
Information  Analysis  Center  Environmental  Research 


Institute  of  Michigan,  September  1994 
14.  Goksin  C.,  Evaluation  of  Tactical  Decision  Programs  for 
Predictions  of  Field  Performance  of  IR  Sensors,  Master  s 
Thesis,  Naval  Postgraduate  School,  September  2000 


146 


15.  Cooper  A.W.,  Lentz  W.J.,  Walker  P.L.,  and  Chan  P.M., 

"Infrared  Polarization  Measurements  of  Ship  Signatures  and 
Background  Contrast,"  SPIE  Proceedings,  Vol .  2223,  1994 

16.  Gregoris  D.J.,  Yu  S.,  Cooper  A.W. ,  and  Milne  E.A., 

"Polarization  Measurements  of  Sun  Glint  from  the  Sea 

Surface",  SPIE  Procedings ,  Vol. 1687,  1992 

17.  Nee  S.M.F.,  Cole  T.  ,  Yoo  C.,  and  Burge  D.  , 

Characterization  of  Infrared  Polarizers",  SPIE  Procedings, 
Vol.  3121,  1997 

18.  Lagaras  S.  E.,  Modeled  Detection  And  Recognition  Range  For 
A  Polarization  Filtered  FLIR  Sensor,  Master's  Thesis,  Naval 
Postgraduate  School,  June  1999 

19.  Maxwell  J.R.,  Beard  J.L.,  and  Due  C.,  "Polarization  In  The 
Thermal  Infrared" ,  paper  presented  at  the  Workshop  on 
Detection,  Discrimination,  and  Classification  of  Targets  in 
Clutter,  13-15  November  1990 

20.  Walker  P.L.,  Lentz  W. J. ,  Cooper  A.W. ,  "Atmospheric  and  Sea 
State  Dependence  of  Polarized  Infrared  Contrast  ,  SPIE 
Proceedings,  Vol.  2469,  April  1995 

21.  Accetta  J.S.,  Shumaker  D.L.,  The  Infrared  and  Electro- 
Optical  Systems  Handbook,  Vol. 4,  SPIE  Optical  Engineering 
Press  and  Infrared  Information  Analysis  Center,  1996 


147 


22.  Shumaker  D.L.,  Wood  J.T.  and  Thacker  C.R.,  Infrared. 
Imaging  Systems  Analysis ,  The  Environmental  Research 
Institute  of  Michigan,  1993 

23.  Guimares  F.C.E.,  Investigation  of  Minimum  Resolvable 

Temperature  Difference  Formulation  For  Polarized  Thermal 
Imaging  Range  Prediction,  Master's  Thesis,  Naval 

Postgraduate  School,  September  1999 

24.  Driggers  R. ,  Cox  P. ,  and  Edwards  T. ,  Introduction  to 
Infrared  and  Electro-Optical  Systems,  Artech  House  Inc., 
MA,  1999 

25.  Thermo Analytics  Inc.,  "MuSES  5.0  Manual", 94X  Airport  Road 
P.O.  Box  66,  Calumet,  Mi  49913,  May  2000 

26.  Cox  C.,  and  Munk  W.  ,  "Bulletin  of  The  Scripps  Institution 
of  Oceanography",  Vol .  6,  1956 

27.  Yu  C.,  Estimate  of  Maximum  Detection  Range  For  FLIR  From 

E0MET95  Measurement  Data,  Master's  Thesis,  Naval 

Postgraduate  School,  December  1997 


148 


INITIAL  DISTRIBUTION  LIST 


No.  Copies 

1.  Defense  Technical  Information  Center . 2 

87725  John  J.  Kingman  Rd. ,  STE  0944 
Ft.  Belvoir,  VA  22060-6218 

2  .  Dudley  Knox  Library . 2 

Naval  Postgraduate  School 
411  Dyer  Rd. 

Monterey,  CA  93943-5101 

3  .  Engineering  and  Technology  Curricular  Office . 1 

Code  34 

Naval  Postgraduate  School 
Monterey,  CA  93943-5107 

4.  Professor  Alfred  W.  Cooper,  Code  PH/Cr  . 3 

Department  of  Physics 

Naval  Postgraduate  School 
Monterey,  CA  93943-5117 

5.  Professor  Ron  J.  Pieper,  Code  EC/Pr  . 1 

Department  of  Physics 

Naval  Postgraduate  School 
Monterey,  CA  93943-5117 

6  .  Dr .  Andreas  K .  Goroch . 1 

Naval  Research  Laboratory 
Monterey,  CA,  93943-5101 

7 .  Kara  Kuwetleri  Komutanligi . 1 

Personel  Daire  Baskanligi 

06200  Bakanliklar,  Ankara 
Turkey 

8 .  Kara  Kuwetleri  Komutanligi  Kutuphanesi . 1 

06200  Bakanliklar,  Ankara 

Turkey 

9.  Kara  Harp  Okulu  Komutanligi  Kutuphanesi  .  1 

06654  Bakanliklar,  Ankara 
Turkey 


149 


10.  Kara  Kuwetleri  Egitim  ve  Doktrin  Komutanligi-EDOK- 

Kutuphanesi  . 

Balgat ,  Ankara 
Turkey- 

11.  Mehmet  Yildirim,  1st  Lieutenant  Turkish  Army  . 

Kara  Kuwetleri  Komutanligi 

Teknik  ve  Pr j .  Yon.  D.  Bsk.ligi 
Tank  Sb.  Pl.Sb. 

Bakanliklar,  Ankara 
Turkey 


150 


